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Abstract 

We review lattice results relevant for pion and kaon physics with the aim of making them 
easily accessible to the particle physics community. Specifically, we review the determina- 
tion of the light-quark masses, the form factor /+(0), relevant for the semileptonic K — s- tt 
transition at zero momentum transfer as well as the ratio fx/ fn of decay constants and 
discuss the consequences for the elements V us and V u d of the CKM matrix. Furthermore, 
we describe the results obtained on the lattice for some of the low-energy constants of 
SU(2) L xSU(2) fl and SU(3) L xSU(3) fl Chiral Perturbation Theory and review the deter- 
mination of the Bk parameter of neutral kaon mixing. We introduce quality criteria and 
use these when forming averages. Although subjective and imperfect, these criteria may 
help the reader to judge different aspects of current lattice computations. Our main re- 
sults are summarized in section 1.2, but we stress the importance of the detailed discussion 
that underlies these results and constitutes the bulk of the present review. 
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1 Introduction 



Flavour physics provides an important opportunity for exploring the limits of the Standard 
Model of particle physics and in constraining possible extensions of theories "Beyond the 
Standard Model". As the LHC explores a new energy frontier, the importance of flavour 
physics will grow still further, in searches for signatures of new physics through precision 
measurements and/or in helping to unravel the theoretical framework behind direct discoveries 
of new particles. The major theoretical limitation arises from the precision with which strong 
interaction effects can be quantified and large-scale numerical simulations of lattice QCD 
provide the opportunity of computing these effects from first principles. In this paper we 
review the current status of lattice results for a variety of physical quantities in low-energy 
physics; our aim is to provide the answer to the frequently posed question "What is currently 
the best lattice value for a particular quantity?" in a way which is readily accessible to 
non-lattice-experts. This is generally not an easy question to answer; different collaborations 
use different lattice actions (discretizations of QCD) with a variety of lattice spacings and 
volumes, and with a range of masses for the u and d quarks. Not only are the systematic 
errors therefore different, but also the methodology used to estimate these uncertainties vary 
between collaborations. Below we summarize the main features of each of the calculations and 
provide a framework for judging and combining the different results. Sometimes it is a single 
result which provides the "best" value; more often it is a combination of results from different 
collaborations. Indeed, the consistency of values obtained using different formulations adds 
significantly to our confidence in the results. 

1.1 Scope of the present review 

The Flavianet Lattice Averaging Group (FLAG) was constituted in November 2007, at a 
general meeting of the European Network on Flavour Physics (Flavianet), with the remit 
of providing the current lattice results to the network's working groups and to the wider 
community. In this paper we focus on physically important quantities in kaon and pion 
physics, specifically the masses of the light quarks, the CKM matrix element V us determined 
from leptonic and semileptonic kaon decays, the Bk parameter of K°-K° mixing and the 
low-energy constants of the SU(2)lx SU(2)r and SU(3)lx SU(3)r chiral Lagrangians. In the 
future, we plan to extend the range of topics covered and to update the material regularly. 
Below, we review papers which appeared before the closing date of the present edition, 28 
February 2011. The averages are based on the results quoted in the tables, which exclusively 
list articles that appeared before this date. Most of the data concern simulations for which 
the masses of the two lightest quarks are set equal. This is indicated by the notation: Nf = 
2 + 1 + 1, for instance, denotes a lattice calculation with four dynamical quark flavours and 
m u = rrid ^ rn s ^ m c . Our review is also available on the FLAG webpage [1], which will 
be updated as new lattice results appear. A compilation of some of the results discussed 
below can also be found on the web page of Laiho, Lunghi and Van de Water [2, 3]. These 
authors restrict themselves to simulations with Nf = 2 + 1 dynamical flavours, but in addition 
cover lattice results for bound states containing charmed or beauty quarks, which are beyond 
the scope of the present review. Our analysis does include simulations with Nf = 2. The 
significance of these data for the physics conclusions to be drawn from the work done on the 
lattice is discussed in section 2.2. 

It is by now generally accepted that, taken together, QCD and QED provide a very 
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accurate approximation for the laws governing the low energy properties of nature - other 
degrees of freedom, such as those of the W and Z bosons, appear to be frozen at low energies. 
Moreover, since QED is infrared stable, the electromagnetic interaction can be dealt with 
perturbatively at low energy, but for QCD, this is not the case. The present review deals 
with the information obtained about the low energy structure of that theory by simulating it 
on a lattice. 

Lattice QCD is a mature field today as a result of a 30-year history of theoretical and 
computational developments, many of which were tested on quenched simulations. Indeed 
the remarkable recent progress in the precision of lattice calculations is due to improved 
algorithms, better computing resources and, last but not least, conceptual developments, such 
as improved actions that reduce lattice artifacts, regularizations which preserve (remnants of) 
chiral symmetry, understanding the finite size effects, non- perturb at ive renormalization, etc. 
For recent discussions of these developments, we refer to [4, 5]. A concise characterization 
of the various discretizations that underlie the results reported in the present review is given 
in appendix A. 1.2 (cf. [6] for a more extended overview). The recent developments in the 
domain of computer algorithms were triggered by [7-10]. For an outline of the progress in 
computer architecture, we refer to [11], which also contains an overview of recent developments 
in simulation-algorithms. 

Each discretization has its merits, but also its shortcomings. One of the main points 
emerging from this review is that the large variety of discretizations employed by the various 
collaborations lead to consistent results, confirming universality within the accuracy reached. 
In our opinion, only the eventual convergence of the various lattice results, obtained with 
different discretizations and methods, can lead to a reliable determination, for which all 
systematic errors can be claimed to be fully under control. 

The lattice spacings reached in recent simulations go down to 0.05 fm or even smaller. In 
that region, growing autocorrelation times may slow down the sampling of the configurations 
[12-16]. Many groups check for autocorrelations in a number of observables, including the 
topological charge, for which a rapid growth of the autocorrelation time is observed if the 
lattice spacing becomes small. In the following, we assume that the continuum limit can be 
reached by extrapolating the existing simulations. 

Lattice simulations of QCD currently involve at most four dynamical quark flavours. To 
connect these with the QCD sector of the Standard Model, which is characterized by Nf = 6, 
one considers a sequence of effective theories, obtained from full QCD by sending one of 
the quark masses after the other to infinity: QCD5, QCD4, ...The series terminates with 
gluodynamics, QCDo- The physical value of the coupling constant is conventionally specified 
in terms of the value in QCD5, in the MS scheme at scale [i = Mz- We characterize the 
masses of the light quarks by their running MS values in QCD4 at scale [i = 2 GeV. 2 

In this language, lattice calculations with Nf = 3 concern QCD3. This theory deviates 
from full QCD by Zweig-suppressed corrections of 0(l/m 2 ). Note that calculations with 
Nf = 2 dynamical flavours often include strange valence quarks interacting with gluons, so 
that bound states with the quantum numbers of the kaons can be studied - what these 
simulations share with QCD2 is that strange sea quark fluctuations are neglected. Likewise, 
calculations done in the quenched approximation, in which sea quark effects are treated as a 

2 Note that this convention is not universally adopted - the scale at which QCD3 is matched to QCD4 is 
not always taken below 2 GeV. The difference between the numerical values in QCD3 and QCD4 is small, but 
with the increase of precision on the lattice, it does become necessary to specify which one of these effective 
theories the quoted running masses refer to (see section 3). 
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mean field, can be made to involve light valence quarks and mesons as well as baryons, while 
in QCDo, the spectrum exclusively contains glueballs. 

For Nf < 3, the connection between the observables and the running coupling constant 
and quark masses of the effective theory is now understood at the non-perturbative level. In 
the case of QCD3, for instance, which is obtained from the QCD sector of the Standard Model 
by sending m t ,mi, and m c to infinity, the running quark masses m u (fi), m,d{(j), m s (fi) and 
the running coupling constant g((j,) can unambiguously be related to the masses M^+, M K +, 
M K o and the mass of one of the baryons that occur in the spectrum, so that the simulation 
yields parameter free results for all of the observables of the effective theory in terms of these 
masses. The lattice results obtained so far, which concern the isospin limit of QCD3, agree 
with experiment, within the accuracy of the calculation. The results reported for the masses 
of the long-lived baryons in BMW 08 [17], for instance, have a precision of better than 4%. In 
particular, if the scale is fixed with Ms, the mass of the f2 can be calculated to an accuracy 
of 1.5%. The fact that the result agrees with experiment indicates that the residual effects 
generated by the neglected heavy degrees of freedom are small. 

We concentrate on pion and kaon observables and assume that, at the precision indicated, 
simulations which do not account for the heavy degrees of freedom represent a sufficiently 
accurate approximation to full QCD. In particular, for fx/ fw and /+(0), where flavour sym- 
metry suppresses the dynamics, we expect the effects generated by the presence of heavy 
sea quarks to be too small to be seen, despite the fact that the level of accuracy reached 
is impressive indeed. We emphasize that the estimates quoted below and labeled Nf = 2 
or Nf = 2 + 1 concern the values of the relevant quantities in full QCD - the label merely 
indicates the method by means of which these estimates are obtained. 

At low energies, the fact that QCD has an approximate, spontaneously broken chiral 
symmetry plays a crucial role. In the limit where the masses of the three lightest quarks are 
set to zero, this symmetry becomes exact, so that the spectrum contains a massless octet 
of Nambu-Goldsone bosons. There are lattice formulations of the theory that preserve this 
symmetry and we refer to appendix A. 1.2 for an overview. We rely on the understanding 
of the dependence on the lattice spacing reached with Symanzik's effective theory [18, 19]. 
This framework can be combined with Chiral Perturbation Theory (%PT), to cope with the 
Nambu-Goldstone nature of the lowest excitations that occur in the presence of light quarks 
[20-40]. The form of the effective theory depends on the discretization used - see appendix 
A. 4 for a brief characterization of the different variants in use. 

As the calculations become more precise, it is important to emphasize the need for contin- 
ued close collaboration between the lattice and xPT communities. Lattice computations are 
generally performed at u and d quark masses which are heavier than the physical ones and 
the results are then extrapolated to the physical point. Although these masses are decreasing 
very significantly with time (see, in particular, the most recent calculations by BMW [41] and 
PACS-CS [42]), it remains true that this extrapolation is one of the most significant sources 
of systematic error. Indeed several lattice collaborations have found that the expansion of 
Mj, M^-, f w , fx to first non-leading order in the quark masses m M ,m^,m s [SU(3) xPT to 
NLO] does not yield an adequate parametrization of the data, and so cannot reliably be used 
for guiding the chiral extrapolation. In part, the problem may be due to the fact that in the 
simulations, the mean mass of u and d is usually significantly heavier than in nature, while 
m s is taken in the vicinity of the physical mass of the strange quark, so that the kaon mass 
is then too heavy for ;\PT to be of use, but it may also indicate that contributions of NNLO 
need to be accounted for to have sufficient accuracy (an analytic approximation for the latter 
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that goes beyond the double-log approximation [43] is proposed in [44, 45]). One method to 
circumvent the impasse is to (a) work at a smaller value of m s , so that the data taken concern 
the region where the kaon mass does not exceed its physical value [46, 47] and (b) use NNLO 
xPT for the extrapolation to the quark masses of physical interest. Alternatively, one may 
try to apply brute force, lower the mean mass of u and d to the physical value and avoid using 
an extrapolation altogether [41, 42]. Which one of these approaches yields better results yet 
remains to be seen, but all of them can teach us something about the low-energy regime of 
QCD and can contribute to a better determination of the low-energy constants of xPT. 

1.2 Summary of the main results 

We end this introduction with a summary of our main results, but we stress that they should 
be considered together with the many comments and observations which accompany their 
determination. The lattice results for the running masses of the three lightest quarks are 
discussed in detail in section 3. Using the MS scheme and fixing the running scale at 2 GeV, 
the lattice data with iVV = 2 + 1 dynamical flavours lead to 

Tfi 

m ud = 3.43 ±0.11 MeV, m s = 94±3MeV, — - = 27.4 ±0.4, (1) 

with m ud = ^(m u + m d ). We emphasize that these results represent a very significant im- 
provement of our knowledge: as compared to the estimates for the light quark masses given 
by the Particle Data Group [48], the uncertainties in the above estimates for m s , m ud and 
m s /m u d are reduced by about an order of magnitude, in all three cases ! Note also that the 
estimates (1) agree with the only Nf = 2 determination that satisfies our quality criteria [49]: 

Tfi 

m ud = 3.6±0.2MeV, m s = 95±6MeV, — - = 27.4 ±0.4. (2) 

m ud 

The lattice simulations are performed with m u = m d , so some additional input is required to 
obtain m u and m d separately (see the discussion in section 3.4). We quote 

Tfi 

m u = 2.19 ±0.15 MeV, m d = 4.67 ± 0.20 MeV , — = 0.47 ±0.04. (3) 

m d 

In the determination of the CKM matrix elements V us and V ud , the term /+(0) (form factor 
relevant for the semileptonic transition K° — > ir~ , evaluated at zero momentum transfer) and 
the ratio fx/ f-w of the kaon and pion decay constants play a central role. Both of these 
quantities can now be determined rather precisely on the lattice. Adding statistical and 
systematic errors in quadrature, the results for /+(0) read 

/+(0) = 0.959 ± 0.005 , (direct, N f = 2 + 1), (4) 

/+(0) = 0.956 ± 0.008 , (direct, N f = 2). 

On the basis of the detailed discussion in section 4.3, we conclude that the value 

/+(0) = 0.956 ± 0.008 (direct) (5) 

represents a conservative estimate for the range permitted by the presently available direct 
determinations of /+(0) in lattice QCD. 
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For fx/ fn, the results are 

JkIU = 1-193 ± 0.005 , (direct , Nf = 2 + 1) , (6) 

fx/fn = 1-210 ± 0.018 , (direct , N f = 2), 

where the first number is the average of three calculations whereas the second stems from 
a single one. We refer to these results, which do not rely on the Standard Model, as direct 
determinations . 

Two precise experimental results constrain the determination of V us and V ud - the observed 
rate for the semileptonic K° — > ir~ decay determines the product | V us \f + (0), while the ratios 
of the kaon's and pion's leptonic widths determine the ratio \V us \fK/\V ud \f n . In section 4.4, 
these results are combined with the lattice determinations of fx/ fn and /+(0) to test the 
unitarity of the CKM matrix. We find that the first row obeys this property within errors: 
\V ud \ 2 + \V US \ 2 + \V ub \ 2 = 1.002 ± 0.015 and 1.036 ± 0.037 for N f = 2 + 1 and N f = 2, 
respectively. 3 In obtaining these results, we have not made use of the recent determinations 
of \V ud \ from super- allowed nuclear /3-decays [50], but these are perfectly consistent with the 
outcome of the lattice calculations. In fact, the unitarity test sharpens considerably if the 
results obtained on the lattice are combined with those found in /3-decay: the interval allowed 
for the sum |^^| 2 + |Kjs| 2 + |^wft| 2 then shrinks by more than an order of magnitude. Since 
this interval includes unity, CKM unitarity passes the test. 

Within the Standard Model of particle physics, the CKM matrix is unitary. As explained 
in section 4.5, this condition and the experimental results for |V^ s |/+(0), \V us \fK/\V ud \f n 
imply three equations for the four quantities \V ud \, \V US \, fx/ fin /+(0)- Thus it is sufficient 
to add one additional constraint, such as the lattice result for fx/fir or /+(0), to determine 
all four of these quantities. Averaging over the numbers obtained for fx/ fir an d /+(0) with 
Nf = 2 + 1 dynamical quark flavours, we find 



\V ud \ = 0.97427 ± 0.00021 , 
/+(0) = 0.9597 ±0.0038, 

The results obtained from the data with Nf = 2 
\V ud \ = 0.97433 ± 0.0042 , 
/ + (0) = 0.9604 ± 0.0075 , 



\V US \ = 0.2254 ±0.0009, (7) 
^f- = 1.1925 ±0.0050. 

J TV 

are the same within errors: 
\V US \ = 0.2251 ±0.0018, (8) 
Ij- = 1.194 ±0.010. 

Jtv 



Since the lattice results are consistent with CKM unitarity, it was to be expected that the 
above values for /+(0) and fx/fn are consistent with those derived from the direct determi- 
nations, listed in equation (6). The main effect of the unitarity constraint is a reduction of 
the uncertainties. 

As far as the low-energy constants are concerned, there are many results coming from 
lattice calculations which are of interest. A coherent picture does not emerge yet for all cases, 
but for two quantities related to the SU(2) chiral expansion, the lattice provides results which 
are competitive with what has been obtained from phenomenology. These are 

F 

£ 3 = 3.2 ± 0.8 and = 1.073 ± 0.015 (9) 

F 

The experimental information on \ V u b\ shows that the contribution from this term is negligibly small. 
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which are estimates based on Nf = 2, Nf = 2+1 and even Nf = 2+1 + 1 calculations. We refer 
the interested reader to Section 5 for further results and all details about the calculations. 

For the Bk parameter of K-K mixing we quote {Bk refers to the MS-scheme at scale 2 
GeV, while Bk denotes the corresponding renormalization group invariant parameter): 

B K = 0.536 + 0.017, 5^ = 0.738 + 0.020 {Nf = 2 + 1), (10) 
B K = 0.516 ± 0.022 , B K = 0.729 ± 0.030 (Nf = 2) . 

1.3 Plan of the paper 

The plan for the remainder of this paper is as follows. In the next section we suggest criteria 
by which the quality of lattice calculations can be judged and compared. We recognize of 
course that these are necessarily subjective and generally give an incomplete picture of the 
quality of a given simulation. Nevertheless we feel that some framework of quality standards 
is necessary. There then follow five sections with physics results: quark masses (section 3), 
V us as determined from leptonic and semileptonic Kaon decays (section 4), the low energy 
constants of SU(2)lxSU(2)r and SU(3)lxSU(3)r, chiral perturbation theory (sections 5.1 
and 5.2 respectively) and the Bk parameter which contains the non-perturbative strong 
interaction effects in K°-K° mixing (section 6). 

The final point we wish to raise in this introduction is a particularly delicate one. As 
stated above, our aim is to make lattice QCD results easily accessible to non-lattice-experts 
and we are well aware that it is likely that some readers will only consult the present paper 
and not the original lattice literature. We consider it very important that this paper is not 
the only one which gets cited when the lattice results which are discussed and analyzed here 
are quoted. Readers who find the review and compilations offered in this paper to be useful 
are therefore kindly requested also to cite the original sources - the bibliography at the end 
of this paper should make this task easier. Indeed we hope that the bibliography will be one 
of the most widely used elements of the whole paper. 

2 Quality criteria 

In this article we present a collection of results and references concerning lattice results of 
physical quantities in low energy particle physics. In addition however, we aim to help the 
reader in assessing the reliability of particular lattice results without necessarily studying the 
original article in depth. We understand that this is a delicate issue and are well aware of 
the risk of making things "simpler than they are". On the other hand, it is rather common 
nowadays that theorists and experimentalists use the results of lattice calculations in order 
to draw the physics conclusions coming from a new measurement or a new analysis of a set of 
measurements, without a critical assessment of the quality of the various calculations. This 
may lead to a substantial underestimate of systematic errors or in the worst case to a distorted 
physics conclusion. We believe that despite the risks mentioned, it is important to provide 
some compact information about the quality of a calculation. 

2.1 Colour code 

In the past, at lattice conferences, some speakers have used a "colour code" to provide such 
information, for example when showing plots [51]. We find that this way of summarizing the 
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quality of a lattice calculation serves well its purpose and adopt it in what follows. We identify 
a number of sources of systematic errors for which a systematic improvement is possible and 
assign to each calculation a colour with respect to each of these: 

when the systematic error has been estimated in a satisfactory manner and convincingly 

shown to be under control; 

when a reasonable attempt at estimating the systematic error has been made, although 
this could be improved; 
■ when no or a clearly unsatisfactory attempt at estimating the systematic error has been 
made. 

The precise criteria used in determining the colour coding is unavoidably time-dependent 
because as lattice calculations become more accurate the standards against which they are 
measured become tighter. Today, our definition of the colour code related to the systematic 
error coming from i) the chiral extrapolation, ii) the continuum extrapolation, Hi) the finite- 
volume effects, and (where applicable) iv) the renormalization and v) the running of operators 
and their matrix elements are as follows: 

• Chiral extrapolation: 

* M nMn < 250 MeV 

250 MeV < M nMn < 400 MeV 

■ M^min > 400 MeV 

It is assumed that the chiral extrapolation is done with at least a three-point analysis 
- otherwise this will be explicitly mentioned in a footnote. In case of nondegeneracies 
among the different pion states M n ,imn stands for a root-mean-squared (RMS) pion 
mass. 

• Continuum extrapolation: 

* 3 or more lattice spacings, at least 2 points below 0.1 fm 
2 or more lattice spacings, at least 1 point below 0.1 fm 

■ otherwise 

It is assumed that the action is 0(a)-improved, i.e. the discretization errors vanish 
quadratically with the lattice spacing - otherwise this will be explicitly mentioned in 
a footnote. Moreover the colour coding criteria for non-improved actions change as 
follows: one lattice spacing more needed. 

• Finite- volume effects: 

* M^minL > 4 or at least 3 volumes 
-<W7r,mm-k > 3 and at least 2 volumes 

■ otherwise 

These criteria apply to calculations in the p-regime, and it is assumed that L m \ n > 2 
fm, otherwise this will be explicitly mentioned in a footnote and a red square will be 
assigned. In case of nondegeneracies among the different pion states M Wj mm stands for 
a root-mean-squared (RMS) pion mass. 

• Renormalization (where applicable): 

* non-perturbative 

2-loop perturbation theory 

■ otherwise 
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• Running (where applicable): 

For scale-dependent quantities, such as quark masses or Bk, it is essential that contact 
with continuum perturbation theory can be established. Various different methods 
are used for this purpose (cf. Appendix A. 3): Regularization-independent Momentum 
Subtraction (RI/MOM), Schrodinger functional, direct comparison with (resummed) 
perturbation theory. In the case of the quark masses, a further approach has been 
proposed recently: determination of m s via the ratio m c /m s . Quite irrespective of 
the particular method used, the uncertainty associated with the choice of intermediate 
renormalization scales in the construction of physical observables must be brought under 
control. This is best achieved by performing comparisons between non-perturbative and 
perturbative running over a reasonably large range of scales. These comparisons were 
initially only made in the Schrodinger functional (SF) approach, but are now also being 
performed in RI/MOM schemes. 

In the framework of the Schrodinger functional, the comparison of the lattice results 
for the relevant renormalization factors with perturbation theory has thoroughly been 
explored. Among the calculations relying on the RI/MOM framework, the most recent 
ones are aiming for a level of control over running and matching which is of comparable 
quality. However, since these approaches are new, we postpone the formulation of 
quantitative criteria until the systematics associated with their use is better understood. 
We mark those data for which information about non-perturbative running checks is 
available and give some details, but do not attempt to translate this into a colour-code. 

The mass of the pion plays an important role in our colour coding, but there are fermion 
action discretizations in which pion states are non-degenerate. This is the case in particu- 
lar for the twisted-mass formulation of lattice QCD. In that discretization, isospin symme- 
try and parity are traded for a chiral symmetry which simplifies the renormalization of the 
theory [52] and, at "maximal twist", eliminates leading 0(a) discretization errors from all 
observables [53]. As a consequence of the loss of isospin symmetry, the charged and neutral 
pions have different masses, even when the u and d quarks are mass degenerate. This mass 
difference is an 0(a 2 ) discretization effect, whose actual size depends on the specific choice 
of the lattice regularization (in both the gauge and fermionic part of the action). In the 
Nf = 2 simulation of [54] it is found 1 — M 7r o/M 7r + ~ 0.2 for a ~ 0.09 fm. Since the charged 
pion mass turns out also to be larger than the mass of the neutral pion, when applying the 
criterion adopted for assessing the quality of the chiral extrapolation to twisted mass fermion 
simulations we will conservatively choose as M njia i n the mass of the charged pion. 

In staggered fermion calculations, 4 the Nambu-Goldstone pion has "taste" partners whose 
masses are those of the Nambu-Goldstone boson plus a discretization error of order a 2 , which 
can be significant at larger lattice spacings. For instance, in MILC's simulations, at a = 
0.15 fm and for a Nambu-Goldstone mass of 241 MeV, the largest taste partner (i.e. the taste 
singlet) mass is 673 MeV and the RMS taste-partner mass is 542 MeV [59]. The situation 
improves significantly as the lattice spacing is reduced: the singlet and RMS masses become 
341 MeV and 334 MeV at a = 0.045 fm for a Nambu-Goldstone of mass 324 MeV. While it 
is possible to pick out the Nambu-Goldstone pion in the valence sector, the contributions of 
the sea correspond to some complicated admixture of the 16 taste partners. MILC reduces to 
one the contribution of the four quark tastes associated with a single quark flavor by taking 

4 We refer the interested reader to a number of good reviews written about the subject [6, 55-58]. 
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the fourth root of the quark determinant, and hence the name rooted staggered fermions. 
In any event, the effective pion mass in the sea can be significantly larger than that of 
the valence Nambu-Goldstone and, in comparing staggered fermion calculations with those 
performed using other discretizations, it makes sense to quote as an M^ jm i n a value which is 
characteristic of the sea pion mass contributions. For most quantities, the RMS taste-partner 
pion mass is a reasonable choice. 

Of course any colour coding has to be treated with caution. First of all we repeat that the 
criteria are subjective and evolving. Sometimes a single source of systematic error dominates 
the systematic uncertainty and it is more important to reduce this uncertainty than to aim 
for green stars for other sources of error. In spite of these caveats we hope that our attempt 
to introduce quality measures for lattice results will prove to be a useful guide. In addition we 
would like to underline that the agreement of lattice results obtained using different actions 
and procedures evident in many of the tables presented below provides further reassurance. 

Finally it is important to stress that the colour coding does not imply any credit for the 
development of theoretical or technical ideas used in the simulations. Frequently new ideas 
and techniques are first tested on smaller lattices, perhaps on quenched configurations or at 
heavier quark masses, to demonstrate their feasibility. They may then be applied on more 
modern lattices by other collaborations to obtain physical results, but since our primary aim 
is to review the latest results concerning quantities of physical interest, we do not try to 
attribute credit to the underlying ideas or techniques. 

For a coherent assessment of the present situation, the quality of the data plays a key 
role, but the colour code cannot be made visible in the figures. Simply showing all data 
on equal footing would give the misleading impression that the overall consistency of the 
information available on the lattice is questionable. As a way out, the figures do indicate the 
quality in a rudimentary way: data for which the colour code is free of red tags are shown as 
green symbols, otherwise they are depicted in red. Note that the pictures do not distinguish 
updates from earlier data obtained by the same collaboration - the reader needs to compare 
the legends to identify closely correlated data. 

2.2 Averages and estimates 

For some observables there may be enough independent lattice calculations of good quality 
that it makes sense to average them and propose such an average as the best current lattice 
number. In order to decide if this situation is realized for a certain observable we rely on 
the colour coding. Unless special reasons are given for making an exception, we restrict the 
averages to data for which the colour code does not contain any red tags. In some cases, 
the averaging procedure nevertheless leads to a result which in our opinion does not cover 
all uncertainties - this happens, in particular, if some of the data have comparatively small 
systematic errors. In such cases, we may provide our estimate as the best current lattice 
number. This estimate is not obtained with a prescribed mathematical procedure, but is 
based on our own assessment of the information collected on the lattice. 

There are two other important criteria which also play a role in this respect, but which 
cannot be represented with the colour coding, because a systematic improvement is not possi- 
ble. These are: i) the publication status, and %%) the number of flavours. As far as the former 
is concerned we adopt the following policy: we consider in our averages only calculations 
which have been published, i.e. which have been endorsed by a referee. The only exception 
to this rule is an obvious update of numbers appearing in a previous publication. Other cases 
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will be listed and identified as such also by a (coloured) symbol, but not used in averages: 

• Publication status: 

A published or plain update of published results 
P preprint 

C conference contribution 

Several active lattice collaborations do their calculations with only two dynamical flavours 
(Nf = 2). We are not aware of visible differences between results obtained in this framework 
and QCD with Nf = 3. At the present level of accuracy, it is thus perfectly meaningful to 
compare results obtained for Nf = 2 with experiment. On the other hand, the two theories 
are intrinsically different; the manner in which quantities like the coupling constant or the 
quark masses depend on the running scale depends on Nf. 

Since we are not aware of an a priori way to quantitatively estimate the difference between 
Nf = 3 and Nf = 2 calculations we will present separate averages (where applicable) for the 
two sets of calculations. Averages of Nf = 3 and Nf = 2 calculations will not be provided. A 
first lattice calculation with Nf = 4 dynamical flavours has recently appeared [60]. While in 
principle it would be cleaner to keep Nf = 4 separate from the rest, we believe that, for the 
quantities under discussion in this review, the effects of the c-quark are below the precision 
of current lattice calculations: if appropriate, we will average Nf = 4 and Nf = 3 results. 

We stress that Nf = 2 calculations have in several better control over systematic 

effects than Nf = 3 calculations, as will be clear in colour coding tables below, and therefore 
play a very important role today in the comparison with experiment. Several calculations 
discussed below provide examples of the relevance of Nf = 2 calculations today, and in 
particular they enter in some of our estimates. In the future, as the accuracy and the control 
over systematic effects in lattice calculations will increase, it will hopefully be possible to 
see a difference between Nf = 2 and Nf = 3 calculations and so determine the size of the 
Zweig-rule violations related to strange quark loops. This is a very interesting issue per se, 
and one which can be quantitatively addressed only with lattice calculations. 

One of the problems that arises when forming averages is that not all of the data sets 
are independent - in fact, some rely on the same ensembles. For this reason, we do not 
average updates with earlier results. Note that, in this respect, the figures give a somewhat 
distorted picture: they include all of the data listed in the tables and, moreover, indicate their 
quality only in a very rudimentary way. Specific problems encountered in connection with 
correlations between different data sets are mentioned in the text. In our opinion, the variety 
of algorithms and configurations used and the fact that we attach conservative errors to our 
estimates ensure that these correlations do not distort our results in a significant manner. 

We take the mean values from the standard procedure, where x 2 1S evaluated by adding 
the statistical and systematic errors in quadrature (data sets for which the error estimate 
only accounts for the statistical uncertainties are not taken into consideration). The standard 
procedure also offers an estimate for the net uncertainty 5 to be attached to the mean: if the 
fit is of good quality (x^ in /dof < 1), calculate 5 from \ 2 = Xmin + 1> otherwise, stretch the 
result obtained in this way by the factor S = y x^ nin /doi. The problem with this recipe is 
that the systematic errors are not stochastic. In particular, applying it to the lattice data on 
nT-ud, m s, ms/iriud and fx/ fn, one arrives at a total error that is smaller than the smallest 
systematic error of the individual data sets. In our opinion, the latter represents a lower 
limit for the total error to be attached to the average. It is not a trivial matter, however, 
to split the total error of the average into a statistical and a systematic part. We are not 
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aware of a generally accepted prescription for doing this and refrain from proposing one in 
this review. In application to the lattice data to be discussed below, where the statistical 
errors are small compared to the systematic ones and thus barely affect the result, replacing 
the quantity 5 by the smallest systematic error of the individual data sets does not lead to 
a significant underestimate of the total error. In the following, we adopt this prescription, 
which has the advantage of being simple and transparent. For a more sophisticated discussion 
of the problem, we refer to [61]. At the precision used in the present review, the numerical 
results obtained for the total errors with the method proposed there are the same. 

3 Quark masses 

Quark masses are fundamental parameters of the Standard Model. An accurate determination 
of these parameters is important for both phenomenological and theoretical applications. The 
charm and bottom masses, for instance, enter the theoretical expressions of several cross 
sections and decay rates, in heavy quark expansions. The up, down and strange quark masses 
govern the amount of explicit chiral symmetry breaking in QCD. From a theoretical point of 
view, the values of quark masses provide information about the flavour structure of physics 
beyond the Standard Model. The Review of Particle Physics of the Particle Data Group 
contains a review of quark masses [48], which covers light as well as heavy flavours. The 
present summary only deals with the light quark masses (those of the up, down and strange 
quarks), but discusses the lattice results for these in more detail. 

Quark masses cannot be measured directly with experiment, because quarks cannot be 
isolated, as they are confined inside hadrons. On the other hand, quark masses are free 
parameters of the theory and, as such, cannot be obtained on the basis of purely theoretical 
considerations. Their values can only be determined by comparing the theoretical prediction 
for an observable, which depends on the quark mass of interest, with the corresponding 
experimental value. What makes light quark masses particularly difficult to determine is the 
fact that they are very small (for the up and down) or small (for the strange) compared to 
typical hadronic scales. Thus, their impact on typical hadronic observables is minute and it 
is difficult to isolate their contribution accurately. 

Fortunately, the spontaneous breaking of SU(3)l<8>SU(3)# chiral symmetry provides ob- 
servables which are particularly sensitive to the light quark masses: the masses of the resulting 
Nambu-Goldstone bosons (NGB), i.e. pions, kaons and etas. Indeed, the Gell-Mann-Oakes- 
Renner relation [62] predicts that the squared mass of a NGB is directly proportional to the 
sum of the masses of the quark and antiquark which compose it, up to higher order mass 
corrections. Moreover, because these NGBs are light and are composed of only two valence 
particles, their masses have a particularly clean statistical signal in lattice QCD calculations. 
In addition to which, the experimental uncertainties on these meson masses are negligible. 

Three flavour QCD has four free parameters: the strong coupling, a s (alternatively Aqcd) 
and the up, down and strange quark masses, m u , and m s . However, present day lattice 
calculations are performed in the isospin limit, and the up and down quark masses get replaced 
by a single parameter: the isospin averaged up and down quark mass, m u d = \{m u + m^). A 
lattice determination of these parameters requires two steps: 

1. Calculations of three experimentally measurable quantities are used to fix the three bare 
parameters. As already discussed, NGB masses are particularly appropriate for fixing 
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the light quark masses. Another observable, such as the mass of a member of the baryon 
octet, can be used to fix the overall scale. It is important to note that most of present 
day calculations are performed at values of m u d which are still substantially larger than 
its physical value, typically four times as large. Thus, reaching the physical up and down 
quark mass point still requires a significant extrapolation. The only exceptions are the 
PACS-CS 08, 09, 10 [42, 63, 64] and BMW 08, 10A [17, 65] calculations discussed below, 
where masses very close to the physical value are reached. The situation is radically 
different for the strange quark. Indeed, modern simulations can include strange quarks 
whose masses bracket its physical value, and only interpolations are needed. 

2. Renormalizations of these bare parameters must be performed to relate them to the 
corresponding cutoff-independent, renormalized parameters. 5 These are short distance 
calculations, which may be performed perturbatively. Experience shows that one-loop 
calculations are particularly unreliable for the renormalization of quark masses, in all 
discretizations. Thus, at least two loops are required to have trustworthy results. Nev- 
ertheless, it is best to perform the renormalizations non-perturbatively, to avoid poten- 
tially large perturbative uncertainties due to neglected higher order terms. 

Of course, in quark mass ratios the renormalization factor cancels, so that this second step is 
no longer relevant. 

3.1 Contributions from the electromagnetic interaction 

As mentioned in section 2, the present review relies on the hypothesis that, at low energies, 
the Lagrangian £ QC d + £ QED describes nature to a high degree of precision. Moreover, we 
assume that, at the accuracy reached by now and for the quantities discussed here, the 
difference between the results obtained from simulations with three dynamical flavours and full 
QCD is small in comparison with the quoted systematic uncertainties. The electromagnetic 
interaction, on the other hand, cannot be ignored. Quite generally, when comparing QCD 
calculations with experiment, radiative corrections need to be applied. In lattice simulations, 
where the QCD parameters are fixed in terms of the masses of some of the hadrons, the 
electromagnetic contributions to these masses must be accounted for. 6 

The electromagnetic interaction plays a crucial role in determinations of the ratio m u /md, 
because the isospin breaking effects generated by this interaction are comparable to those from 
m u rrid (see subsection 3.4). In determinations of the ratio m s /m u d, the electromagnetic 
interaction is less important, but at the accuracy reached, it cannot be neglected. The reason 
is that, in the determination of this ratio, the pion mass enters as an input parameter. Because 

5 Throughout this review, the quark masses m u , nid and m s refer to the MS scheme at running scale 
H — 2 GeV and the numerical values are given in MeV units. 

6 Since the decomposition of the sum £qcd + £qed into two parts is not unique, specifying the QCD part 
requires a convention. In order to give results for the quark masses in the Standard Model at scale [i = 2 GeV, 
on the basis of a calculation done within QCD, it is convenient to match the two theories at that scale. We 
use this convention throughout the present review. Note that a different convention is used in the analysis of 
the precision measurements carried out in low energy pion physics. The result ao M n + = 0.220(5) (2) (6) of the 
NA48/2 collaboration [66] for the 1 — S-wave -kit scattering length, for instance, also concerns QCD in the 
isospin limit, but refers to the convention where M n is identified with M n +. When comparing lattice results 
with experiment, it is important to fix the QCD parameters in accordance with the convention used in the 
analysis of the experimental data (for a more detailed discussion, see [67-70]). 



15 



M n represents a small symmetry breaking effect, it is rather sensitive to the perturbations 
generated by QED. 

We distinguish the physical mass Mp, P £ {vr + , 7r°, K + , K }, from the mass Mp within 
QCD alone. The e.m. self-energy is the difference between the two, Mj, = Mp — Mp. Since 
the self-energy of the Nambu-Goldstone bosons diverges in the chiral limit, it is convenient 
to replace it by the contribution of the e.m. interaction to the square of the mass, 

A p = Mp — Mp = 2 MpMl + C»(e 4 ) . (11) 

The main effect of the e.m. interaction is an increase in the mass of the charged particles, 
generated by the photon cloud that surrounds them. The self-energies of the neutral ones 
are comparatively small, particularly for the Nambu-Goldstone bosons, which do not have a 
magnetic moment. Dashen's theorem [71] confirms this picture, as it states that, to leading 
order (LO) of the chiral expansion, the self-energies of the neutral NGBs vanish, while the 
charged ones obey A^ K+ = A^ + . It is convenient to express the self-energies of the neutral 
particles as well as the mass difference between the charged and neutral pions within QCD 
in units of the observed mass difference, A w = M 2 + — M 2 n : 

A^o = e w o , A^ ee e K o An , Ml + - M 2 = e m A n . (12) 

In this notation, the self-energies of the charged particles are given by 

Al + = (1 + e^o - e m ) A,, , A^ + = (1 + e + e K o - e m ) A n , (13) 

where the dimensionless coefficient e parametrizes the violation of Dashen's theorem, 

A^-A^-A^ + A^EEeA,. (14) 

Any determination of the light quark masses based on a calculation of the masses of ir + , K + 
and K° within QCD requires an estimate for the coefficients e, e n o, e K o and e m . 

The first determination of the self-energies on the lattice was carried out by Duncan, 
Eichten and Thacker [72]. Using the quenched approximation, they arrived at Mj^ + — M^. = 
1.9 MeV. Actually, the parametrization of the masses given in that paper yields an estimate for 
all but one of the coefficients introduced above (since the mass splitting between the charged 
and neutral pions in QCD is neglected, the parametrization amounts to setting e m = ab 
initio). Evaluating the differences between the masses obtained at the physical value of the 
electromagnetic coupling constant and at e = 0, we obtain e = 0.50(8), e n o = 0.034(5) and 
e^o = 0.23(3). The errors quoted are statistical only: an estimate of lattice systematic errors 
is not possible from the limited results of [72] . The result for e indicates that the violations of 
Dashen's theorem are sizable: according to this calculation, the non-leading contributions to 
the self-energy difference of the kaons amount to 50% of the leading term. The result for the 
self-energy of the neutral pion cannot be taken at face value, because it is small, comparable 
with the neglected mass difference M w + — M^o. To illustrate this, we note that the numbers 
quoted above are obtained by matching the parametrization with the physical masses for 7r°, 
K + and K°. This gives a mass for the charged pion that is too high by 0.32 MeV. Tuning 
the parameters instead such that M n + comes out correctly, the result for the self-energy of 
the neutral pion becomes larger: e n o = 0.10(7). Also, the uncertainties due to the systematic 
errors of the lattice calculation need to be taken into account. 
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In an update of this calculation by the RBC collaboration [73] (RBC 07), the electromag- 
netic interaction is still treated in the quenched approximation, but the strong interaction is 
simulated with Nt = 2 dynamical quark flavours. The quark masses are fixed with the physical 
masses of ir°, K + and K°. The outcome for the difference in the electromagnetic self-energy 
of the kaons reads Mj, + — M^. = 1.443(55) MeV. This corresponds to a remarkably small vi- 
olation of Dashen's theorem. Indeed, a recent extension of this work to Nf = 2 + 1 dynamical 
flavours [74] leads to a significantly larger self-energy difference: M 1 K+ — MX = 1.87(10) MeV, 
in good agreement with the estimate of Eichten et al. Expressed in terms of the coefficient e 
that measures the size of the violation of Dashen's theorem, it corresponds to e = 0.5(1). 

More recently, the BMW collaboration has reported preliminary results for the electromag- 
netic self-energies, based on simulations in which a U(l) degree of freedom is superimposed 
on their Nt = 2 + 1 QCD configurations [75]. The result for the kaon self-energy difference 
reads — M^. = 2.2(2) MeV, that is e = 0.66(14), where the error indicates the statistical 
uncertainties. 

The input for the electromagnetic corrections used by the MILC collaboration is specified 
in [76]. In their analysis of the lattice data, e^o, e K o and e m are set equal to zero. For the 
remaining coefficient, which plays a crucial role in determinations of the ratio m u /mci, the very 
conservative range e = 1±1 was used in MILC 04 [77], while in more recent work, in particular 
in MILC 09 [6] and MILC 09A [59], this input is replaced by e = 1.2 ± 0.5, as suggested by 
phenomenological estimates for the corrections to Dashen's theorem [78, 79]. First results 
of an evaluation of the electromagnetic self-energies based on Nf = 2 + 1 dynamical quarks 
in the QCD sector and on the quenched approximation in the QED sector are also reported 
by the MILC collaboration [80]. Further calculations using staggered quarks and quenched 
photons are underway [81]. 

The effective Lagrangian that governs the self-energies to next-to-leading order (NLO) of 
the chiral expansion was set up in [82]. The estimates in [78, 79] are obtained by replacing 
QCD with a model, matching this model with the effective theory and assuming that the 
effective coupling constants obtained in this way represent a decent approximation to those 
of QCD. For alternative model estimates and a detailed discussion of the problems encoun- 
tered in models based on saturation by resonances, see [83-85]. In the present review of the 
information obtained on the lattice, we avoid the use of models altogether. 

There is an indirect phenomenological determination of e, which is based on the decay 
rj — > 3ir and does not rely on models. The result for the quark mass ratio Q obtained from a 
dispersive analysis of this decay implies e = 0.70(28) (see section 3.4). While the values found 
on the lattice [72-75] are lower, the phenomenological estimate used by MILC 09A [59] is on 
the high side. In the following, we take the central value from 7y-decay, but stretch the error, 
so that the numbers given in [59, 72-75] are all within one standard deviation: e = 0.7(5). 

We add a few comments concerning the physics of the self-energies and then specify the 
estimates used as an input in our analysis of the data. The Cottingham formula [86] represents 
the self-energy of a particle as an integral over electron scattering cross sections; elastic as 
well as inelastic reactions contribute. For the charged pion, the term due to elastic scattering, 
which involves the square of the e.m. form factor, makes a substantial contribution. In 
the case of the 7r°, this term is absent, because the form factor vanishes on account of charge 
conjugation invariance. Indeed, the contribution from the form factor to the self-energy of the 
7r + roughly reproduces the observed mass difference between the two particles. Furthermore, 
the numbers given in [87-89] indicate that the inelastic contributions are significantly smaller 
than the elastic contributions to the self-energy of the tt + . The low energy theorem of Das, 
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Guralnik, Mathur, Low and Young [90] ensures that, in the limit m u ,m^ — > 0, the e.m. self- 
energy of the 7T° vanishes, while the one of the ir + is given by an integral over the difference 
between the vector and axial spectral functions. The estimates for e w o obtained in [72] are 
consistent with the suppression of the self-energy of the tt° implied by chiral SU(2)xSU(2). 
In our opinion, e n o = 0.07(7) is a conservative estimate for this coefficient. The self-energy 
of the K° is suppressed less strongly, because it remains different from zero if m u and 
are taken massless and only disappears if m s is turned off as well. Note also that, since 
the e.m. form factor of the K° is different from zero, the self-energy of the K° does pick up 
an elastic contribution. The lattice result for e^-o indicates that the violation of Dashen's 
theorem is smaller than in the case of e. In the following, we use e^-o = 0.3(3). 

Finally, we consider the mass splitting between the charged and neutral pions in QCD. 
This effect is known to be very small, because it is of second order in m u — m^. There 
is a parameter-free prediction, which expresses the difference — M^ in terms of the 
physical masses of the pseudoscalar octet and is valid to NLO of the chiral perturbation 
series. Numerically, the relation yields e m = 0.04 [91], indicating that this contribution does 
not play a significant role at the present level of accuracy. We attach a conservative error also 
to this coefficient: e m = 0.04(2). The lattice result for the self-energy difference of the pions, 
reported in [74], M 7 + — M 7 = 4.50(23) MeV, agrees with this estimate: expressed in terms 
of the coefficient e m that measures the pion mass splitting in QCD, the result corresponds 
to e m = 0.04(5). The corrections of next-to-next-to- leading order (NNLO) have been worked 
out [92] , but the numerical evaluation of the formulae again meets with the problem that the 
relevant effective coupling constants are not reliably known. 

In summary, we use the following estimates for the e.m. corrections: 

e = 0.7(5) , € w0 = 0.07(7) , e K o = 0.3(3) , e m = 0.04(2) . (15) 

While the range used for the coefficient e affects our analysis in a significant way, the numerical 
values of the other coefficients only serve to set the scale of these contributions. The range 
given for e n o and e K o may be overly generous, but because of the exploratory nature of the 
lattice determinations, we consider it advisable to use a conservative estimate. 

Treating the uncertainties in the four coefficients as statistically independent and adding 
errors in quadrature, the numbers in equation (15) yield the following estimates for the 
e.m. self-energies, 

Mj + =4.7(3) MeV, Mj = 0.3(3) MeV , Mj + - Mj = 4.4(1) MeV , (16) 
= 2.5(7) MeV, = 0.4(4) MeV , M^. + — M^ = 2.1(6) MeV , 

and for the pion and kaon masses occurring in the QCD sector of the Standard Model, 

M w + = 134.8(3) MeV , M^o = 134.6(3) MeV , M n + - M^o = 0.2(1) MeV , (17) 
M K + = 491.2(7) MeV , M K o = 497.2(4) MeV , M K + - M K o = -6.1(6) MeV . 

The self-energy difference between the charged and neutral pion involves the same coefficient 
e m that describes the mass difference in QCD - this is why the estimate for Mj + — Mj is so 
sharp. 

3.2 Pion and kaon masses in the isospin limit 

As mentioned above, most of the lattice calculations concerning the properties of the light 
mesons are performed in the isospin limit of QCD (m u — — > at fixed m u + mj. We 
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denote the pion and kaon masses in that limit by M n and Mk, respectively. Their numerical 
values can be estimated as follows. Since the operation u f-> d interchanges tt + with 7r~ and 
K + with K°, the expansion of the quantities and + M^ ) in powers of m u — irid 

only contains even powers. As shown in [93] , the effects generated by m u — rrid in the mass 
of the charged pion are strongly suppressed: the difference — M n represents a quantity 
of 0[(m u — md) 2 {m u + md)] and is therefore small compared to the difference M 2 + — M 2 , 

for which an estimate was given above. In the case of \{M^+ + — M^, the expansion 

does contain a contribution at NLO, determined by the combination 2Lg — L5 of low energy 
constants, but the lattice results for that combination show that this contribution is very 
small, too. Numerically, the effects generated by m u — md in and in \ +M^ ) are 

negligible compared to the uncertainties in the electromagnetic self-energies. The estimates 
for these given in equation (15) thus imply 

M w = 134.8(3) MeV , M K = 494.2(5) MeV . (18) 

This shows that, for the convention used above to specify the QCD sector of the Standard 
Model, and within the accuracy to which this convention can currently be implemented, 
the mass of the pion in the isospin limit agrees with the physical mass of the neutral pion: 
M w - M n0 = -0.2(3) MeV. 

3.3 Lattice determination of m s and m u d 

We now turn to a review of the lattice calculations of the light quark masses and begin with 
m s , the isospin averaged up and down quark mass, m u d, and their ratio. Most groups quote 
only m u d, not the individual up and down quark masses. We then discuss the ratio m u /md 
and the individual determination of m u and m^. 

Quark masses have been calculated on the lattice since the mid-nineties. However early 
calculations were performed in the quenched approximation, leading to unquantifiable sys- 
tematics. Thus in the following, we only review modern, unquenched calculations, which 
include the effects of light sea quarks. 

Tables 1 and 2 list the results of Nf = 2 and Nf = 2 + 1 lattice calculations of m s and m u d- 
These results are given in the MS scheme at 2 GeV, which is standard nowadays. The tables 
also show the colour coding of the calculations leading to these results. The corresponding 
results for m s /m u d are given in Table 3. As indicated in the Introduction, we treat Nf = 2 
and Nf = 2 + 1 calculations separately. The latter include the effects of a strange sea quark, 
but the former do not. 

3.3.1 Nf = 2 lattice calculations 

We begin with Nf = 2 calculations. A quick inspection of Table 1 indicates that only the 
most recent calculation, ETM 10B [94], controls all systematic effects. Only ETM 10B [94] 
and ETM 07 [49] really enter the chiral regime, with pion masses down to about 300 MeV. 
All the others have significantly more massive pions, the lightest being about 430 MeV, in 
the calculation by CP-PACS 01 [101]. Moreover, the latter calculation is performed on very 
coarse lattices, with lattice spacings a > 0.11 fm and only one loop perturbation theory is 
used to renormalize the results. 

One calculation of m s which performs quite well in terms of colour coding is that of AL- 
PHA 05 [98] . This is due to the fact that both non-perturbative running and non-perturbative 
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CP-PACS 01 


noii 


A i 


■ 


★ 



* 
* 
* 
* 


a 3.6(1)(2) 

- 4.452(81) (38) (_+° 7 ) 

4.25(23)(26) 

3.85(12)(40) 


95(2)(6) 

119.5(5.6)(7.4) 
105(3) (9) 


* 


4.08(23) (19) (23) 


111(6) (4) (6) 


* 
* 


- 4.3(4)(tS:£) 
b 


101(8)(±g 5 ) 
97(4)(18) s 


* 


- 4.7(2)(3) 


119(5)(8) 


■ 
■ 


- 3.223(t£) 

- 3.45(10)(±iJ) 


84.5(+l 2 7 °) 
89(2)(1 2 6 )* 



' The calculation includes quenched e.m. effects. 

' The data used to obtain the bare value of m s are from RBC/UKQCD 04 [99] 

* This value of m s was obtained using the kaon mass as input. If the <j> meson mass is used instead, the 
authors find m s — 90l n . 

a The masses are renormalized non-perturbatively at scales 1/a ~ 2 -f- 3GeV in the Nf — 2 RI/MOM 

scheme. In this scheme, non-perturbative and N 3 LO running for the quark masses are shown to agree 

from 4 GeV down 2 GeV to better than 3% [102]. 
b The masses are renormalized and run non-perturbatively up to a scale of 100 GeV in the Nf = 2 SF 

scheme. In this scheme, non-perturbative and NLO running for the quark masses are shown to agree 

well from 100 GeV all the way down to 2 GeV [98]. 



Table 1: Nf = 2 lattice results for the masses m u d and m s , together with the colour coding 
of the calculation used to obtain these. If information about non-perturbative running is 
available, this is indicated in the column "running", with details given at the bottom of the 
table. 

renormalization are performed in a controlled fashion, using Schrodinger functional methods. 
However, the lightest sea pion in the calculation is 600 MeV, significantly above our thresholds 
for a controlled chiral extrapolation. In that sense, this calculation should be viewed more 
as a means of giving a first estimate of possible sea quark corrections to the equally precise 
quenched calculation ALPHA 99 [112], than a full fledged dynamical fermion calculation. 

The conclusion of our analysis of Nf = 2 calculations is that the results of ETM 10B [94] 
(which updates and extends ETM 07 [49]), is the only one to date which satisfies our selection 
criteria. Thus we do not offer an average of Nf = 2 results for the strange and isospin averaged 
quark masses and their ratios, but simply quote the results of ETM 10B [94]: 

772 

Nf = 2 : m s = 95(2) (6) MeV, m u d = 3.6(1)(2) MeV , — - = 27.3(5) (7) . (19) 

The combined statistical and systematic errors on these results are 7%, 6% and 3% respec- 
tively. The errors are smaller in the ratio because many systematics cancel, most notably 
those associated with renormalization and the setting of the scale. Moreover, because of 
statistical correlations between m s and m u d, those errors are also reduced in the ratio. 
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PACS-CS 10 
MILC 10A 
HPQCD 10 
BMW 10A, 10B+ 
RBC/UKQCD 10 A 
Blum 10 t 
PACS-CS 09 
HPQCD 09 
MILC 09A 
MILC 09 
PACS-CS 08 
RBC/UKQCD 08 
CP-PACS/ 
JLQCD 07 
HPQCD 05 
MILC 04, HPQCD/ 
MILC/UKQCD 04 



[«.-> 



[64] 
[103] 
[104] 
105] 
[106] 
[74] 
[42] 
[107] 
[59] 

[6] 
[63] 
[108] 

[109] 

[110] 

[77, 1111 



P 
C 
A 
P 
P 
P 
A 
A 
C 
A 
A 
A 

A 

A 

A 



★ 
★ 
★ 

★ 



★ 
★ 
★ 



★ 
★ 
★ 



* 
* 
* 
★ 



★ 
* 
★ 
■ 

* 

★ 



★ 

★ 
★ 
★ 
★ 
★ 
★ 



2.78(27) 
3.19(4) (5) (16) 
3.39(6)* 
3.469(47) (48) 
3.59(13) (14) (8) 
3.44(12) (22) 
2.97(28)(3) 
3.40(7) 

3.25 (1)(7)(16)(0) 
3.2(0) (1) (2) (0) 
2.527(47) 
3.72(16) (33) (18) 

3.55(19)(±f ) 



86.7(2.3) 

92.2(1.3) 
95.5(1. 1)(1.5) 
96.2(1.6) (0.2) (2.1) 
97.6(2.9)(5.5) 
92.75(58) (95) 
92.4(1.5) 

89.0(0.2) (1.6) (4.5) (0.1) 
88(0) (3) (4) (0) 
72.72(78) 

107.3(4.4) (9.7) (4.9) 
90.1(4.3)(ti 6 3 7 ) 



3.2(0) (2) (2) (0)* 87(0) (4) (4) (0)* 



2.8(0) (1) (3) (0) 



76(0) (3) (7) (0) 



* Value obtained by combining the HPQCD 10 result for m 3 with the MILC 09 result for m 3 /m u d. 
+ The fermion action used is tree-level improved. 

* The calculation includes quenched e.m. effects. 

* The bare numbers are those of MILC 04. The masses are simply rescaled, using the ratio of the 2-loop to 
1-loop renormalization factors. 

a The masses are renormalized and run non-perturbatively up to a scale of 40 GeV in the Nf — 3 SF 

scheme. In this scheme, non-perturbative and NLO running for the quark masses are shown to agree well 

from 40 GeV all the way down to 3 GeV [64]. 
b The masses are renormalized and run non-perturbatively up to a scale of 4 GeV in the Nf = 3 RI/MOM 

scheme. In this scheme, non-perturbative and N 3 LO running for the quark masses are shown to agree 

from 6 GeV down to 3 GeV to better than 1% [105]. 
c The masses are renormalized non-perturbatively at a scale of 2 GeV in a couple of Nf = 3 RI/SMOM 

schemes. A careful study of perturbative matching uncertainties has been performed by comparing results 

in the two schemes in the region of 2 GeV to 3 GeV [106]. 

Table 2: jV^ = 2 + 1 lattice results for the masses m u d and m s , together with the colour 
coding of the calculation used to obtain these. If information about non-perturbative running 
is available, this is indicated in the column "running", with details given at the bottom of 
the table. 



3.3.2 Nf = 2 + 1 lattice calculations 

We turn now to Nf = 2 + 1 calculations. These and the corresponding results are summarized 
in Table 2 and 3. Somewhat paradoxically, these calculations are more mature than those with 
Nf = 2. This is thanks, in large part, to the head start and sustained effort of MILC, which 
has been performing Nf = 2 + 1 rooted staggered fermion calculations for the past ten or so 
years. They have covered an impressive range of parameter space, with lattice spacings which, 
today, go down to 0.045 fm and valence pion masses down to approximately 180 MeV [59]. 
The most recent updates, MILC 10A [103] and MILC 09A [59], include significantly more 
data and use two loop renormalization. Since these data sets subsume those of their previous 
calculations, these latest results are the only ones that must be kept in any world average. 

We recall that valence and sea quarks are treated differently in the rooted staggered 
approach, as discussed in section 2, and that the lightest effective sea-pion mass of all of 
their ensembles is about 260 MeV rather than 180 MeV. Regarding the specific subject of 
light quark masses, a weak point of present day staggered fermion calculations is the use 
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BMW 10A, 10B+ 


[65, 105] 


P 


RBC/UKQCD 10A 


[106] 


P 


Blum 10* 


[74] 


P 


PACS-CS 09 


[42] 


A 


MILC 09A 


[59] 


C 


MILC 09 


[6] 


A 


PACS-CS 08 


[63] 


A 


RBC/UKQCD 08 


[108] 


A 


MILC 04, HPQCD/ 


[77, 111] 


A 


MILC/UKQCD 04 



★ 


★ 


* 


27.53(20) (8) 


• 


• 


* 


26.8(0.8) (1.1) 


• 


■ 


• 


28.31(0.29)(1.77) 


★ 


■ 


■ 


31.2(2.7) 


• 


★ 


* 


27.41(5) (22) (0) (4) 


• 


★ 


★ 


27.2(1) (3) (0)(0) 


★ 


■ 


■ 


28.8(4) 


• 


■ 


★ 


28.8(0.4)(1.6) 


• 


• 


• 


27.4(1) (4) (0)(1) 



ETM 10B 


[94] 


A 


• 




• 


27.3(5)(7) 


RBC 07* 


[73] 


A 


■ 


■ 


★ 


28.10(38) 


ETM 07 


[49] 


A 


• 


■ 


• 


27.3(0.3) (1.2) 


QCDSF/UKQCD 06 


[96] 


A 


■ 


★ 


■ 


27.2(3.2) 



+ The fermion action used is tree-level improved. 
' The calculation includes quenched e.m. effects. 

Table 3: Lattice results for m s /m u ^, together with the colour coding of the calculation used 
to obtain them. The upper half of the table corresponds to Nf = 2 + 1 calculations, the lower 
half to Nf = 2 results. 

of perturbative renormalization. As already noted, it is preferable, a priori, to renormal- 
ize non-perturbatively, to avoid introducing uncontrolled systematic errors from neglected, 
higher order contributions. That this is a problem, in practice, for the staggered calculations 
can be seen from the following observation. In MILC 04 [77] and HPQCD/MILC/UKQCD 
04 [111], only one loop perturbation theory was used. 0{a 2 s ) corrections were included in the 
HPQCD 05 [110] update. The only difference between this update and the earlier calculation 
is the inclusion of these corrections: the bare lattice results are the same in both. Yet, the 
0{a 2 s ) correction increased the renormalized value of m s by 11 MeV, i.e. 14%. In light of this, 
one may question whether the perturbative expansion is under control. Moreover, prelimi- 
nary results including non-perturbative renormalization suggest that the omitted terms may 
contribute an additional increase of ~ 20% in m s [113]. 

Of all the entries in Tables 2 and 3, Blum 10 [74] is the only one which provides an 
improved determination of the electromagnetic self-energies on the lattice and points the way 
for future promising developments: still, at present, these data do not pass the selection 
criteria formulated in section 2.2 as the authors did not investigate all sources of systematic 
error. There are four other calculations with iVy = 2 + 1, however, that do and that we discuss 
now. 

The most recent one is BMW 10A, 10B [65, 105] which, as shown by the color coding in 
Tables 2 and 3 has addressed all sources of systematic effects. One of the most impressive 
aspects of this calculation is that they have reached the physical up and down quark mass 
by interpolation and did not need any chiral extrapolation. Moreover, their calculation was 
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performed at five lattice spacings ranging from 0.054 to 0.116 fm, with full non-perturbative 
renormalization and running and in volumes of up to (6 fm) 3 guaranteeing that all sources of 
systematic error are controlled. However, according to the rules set forth in section 2.2, these 
results will have to be published before they can be included in our averages. 

The second calculation that we discuss is PACS-CS 10 [64], which updates 09 [42]: these 
calculations have been performed with 0(o)-improved Wilson fermions on lattices with a ~ 
0.09 fm. The important aspect of this calculation is that it also probes pion masses down 
to ~ 135 MeV, i.e. down to the physical mass point. This is achieved by reweighting 
the simulations performed in PACS-CS 08 [63] at M n ~ 160 MeV. If adequately controlled, 
this procedure eliminates the need to extrapolate to the physical mass point and, hence, the 
corresponding systematic error. Moreover, renormalization of quark masses is implemented 
non-perturbatively, through the Schrodinger functional method [114]. As it stands, the main 
drawback of the calculation, which makes the inclusion of its results in a world average of 
lattice results inappropriate at this stage, is that for the lightest quark mass the volume is 
very small, corresponding to LAl n ~ 2.0, a value for which finite volume effects will be difficult 
to control. Another problem is that the calculation was performed at a single lattice spacing, 
forbidding a continuum extrapolation. Furthermore, it is unclear at this point what might be 
the systematic errors associated with the reweighting procedure. The PACS-CS collaboration 
plans to perform physical point simulations with larger and finer lattices. 

The third calculation of particular interest is RBC/UKQCD 10A [106] (which updates 
RBC/UKQCD 08 [108]). This calculation uses the theoretically appealing domain wall 
fermions. It reaches down to 330 MeV pion masses in the sea and to 242 MeV valence pions. 
Moreover, it makes use of non-perturbative renormalization a la Rome-Southampton [115]. 
The latest update removes the main drawback of the previous calculation by adding a second, 
finer lattice spacing of about 0.09 fm to the first rather coarse lattice spacing of 0.11 fm, and 
allows a first estimate of the systematic error from the continuum extrapolation. 

Finally we come to the fourth calculation which satisfies our selection criteria, HPQCD 10 [104] 
(which updates HPQCD 09 [107]). The strange quark mass is computed using a precise de- 
termination of the charm quark mass, m c (m c ) = 1.273(6) GeV [104, 116], whose accuracy is 
better than 0.5%, and a calculation of the quark mass ratio m c /m s = 11.85(16) [107], which 
achieves a precision slightly above 1%. The determination of m s via the ratio m c /m s circum- 
vents the problem of running by calculating a scale-invariant ratio non-perturbatively. The 
MILC 09 determination of the quark mass ratio m s /m u d [6] is also used in HPQCD 10 [104] 
to calculate m U( i- 

The high precision quoted by HPQCD 10 on the strange quark mass relies in large part 
on the precision reached in the determination of the charm quark mass [104, 116]. This 
calculation uses a novel approach based on the lattice determination of moments of charm- 
quark pseudoscalar, vector and axial- vector correlators. These moments are then combined 
with four-loop results from continuum perturbation theory to obtain a determination of the 
charm quark mass in the MS scheme. In the preferred case, in which pseudoscalar correlators 
are used for the analysis, there are no renormalization factors required, since the corresponding 
axial-vector current is partially conserved in the staggered lattice formalism. Overall, the 
0.5% accuracy claimed in the calculation of m c is certainly impressive, since it implies that 
all sources of uncertainty, coming from the scale determination, discretization effects (the 
bare charm mass in lattice units is as large as am c ~ 0.85), perturbation theory as well as 
non-perturbative effects in the estimate of moments, are all controlled at the per mille level. 

This discussion leaves us with three results to average for m s , those of MILC 09A [59], 
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Figure 1: Mass of the strange quark (MS scheme, running scale 2 GeV). The upper part shows 
the lattice results listed in Tables 1 and 2. Full and empty symbols correspond to simulations 
with Nf = 2 + 1 and Nj = 2, respectively. Diamonds represent results based on perturbative 
renormalization, while squares indicate that, in the relation between the lattice-regularized 
and renormalized MS masses, non-perturbative effects are accounted for. The vertical bands 
indicate our estimates (19) and (20). The lower part shows recent determinations obtained 
from the evaluation of sum rules, together with the earliest result based on this method, as 
well as the most recent estimate of the Particle Data Group. 



RBC/UKQCD 10A [106] and HPQCD 10 [104]. Treating them as independent measurements 
and simply adding the statistical and systematic errors in quadrature, we obtain a good fit, 
with x 2 = 2.5 for 3 data points and 1 free parameter. With the prescription specified in 
section 2.2, the average becomes m s = 92.8(1.3) MeV. When repeating the exercise for m u d, 
we replace MILC 09A by the more recent analysis reported in MILC 10A [103] and again 
obtain a good fit: \ 2 = 2.3 for 2 degrees of freedom and m u d = 3.38(6) MeV for the average. 
The outcome of the averaging procedure thus amounts to a determination of m s and m uc i of 
1.4%. and 1.8%, respectively. 

The results in BMW 10A, 10B [65, 105], which are not included in the averaging procedure 
on account of our criteria concerning the publication status, offer an independent, remarkably 
sharp determination of m u d and m s . These data are perfectly consistent with the averages 
quoted above. 

The heavy sea quarks affect the determination of the light quark masses only through 
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Figure 2: Mean mass of the two lightest quarks, m u d 
various symbols is explained in the caption of Fig. 1. 



+ rrid). The meaning of the 



contributions of order 1/mf., which moreover are suppressed by the Zweig-rule. We expect 
these contributions to be small, but do not know of a reliable quantitative evaluation. The 
problem originates in the fact that the relation between the parameters of QCD3 and those of 
full QCD can currently be analyzed only in the framework of perturbation theory. The f3- and 
7-functions, which control the renormalization of the coupling constants and quark masses, 
respectively, are known to four loops [124-127]. The precision achieved in this framework for 
the decoupling of the t- and 6-quarks is excellent, but the c-quark is not heavy enough: at the 
percent level, the corrections of order 1/mj? cannot be neglected and the decoupling formulae of 
perturbation theory do not provide a reliable evaluation, because the scale m c (m c ) ~ 1.27 GeV 
is too low for these formulae to be taken at face value. Consequently, the accuracy to which 
it is possible to identify the running masses of the light quarks of full QCD in terms of those 
occurring in QCD3 is limited. For this reason, it is preferable to characterize the masses m u , 
rrid, m s in terms of QCD4, where the connection with full QCD is under good control. The 
role of the c-quarks in the determination of the light quark masses will soon be studied in 
detail - some simulations with 2+1+1 dynamical quarks have already been carried out [60]. 

A crude upper bound on the size of the effects due to the neglected heavy quarks can 
be established within the Nt = 2 + 1 simulations themselves, without invoking perturbation 
theory. In [17] it is found that when the scale is set by M=, the result for Ma agrees well 
with experiment within the 2.3% accuracy of the calculation. Because of the very strong 
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correlations between the statistical and systematic errors of these two masses, we expect the 
uncertainty in the difference Mg — M\ to also be of order 2%. To leading order in the chiral 
expansion, this mass difference is proportional to m s — m u d- We conclude that the agreement 
of Nf = 2 + 1 calculations with experiment yields an upper bound on the sensitivity of m s to 
heavy sea quarks of order 2%. 

In order to stay on the conservative side in this rapidly developing field, our final estimates 
for m u d and m s come with sizable uncertainties: 

N f = 2 + 1: m ud = 3.43(11) MeV , m s = 94(3) MeV . (20) 

As discussed in section 2.2, these numbers are not obtained from a mathematical prescription, 
but are chosen so as to cover all of the high quality data discussed in the present section (note 
also that the central values are slightly higher than those of the averages discussed above) . 

The estimates (20) represent the conclusions we draw from the information gathered on 
the lattice until now. They are shown as vertical bands in Figures 1 and 2, together with the 
Nf = 2 results (19) and a few estimates from other sources. The values quoted in the 2010 
edition of the Review of Particle Properties [122] are m u d = 3.9(9), m s = 101(jJi) MeV. The 
errors attached to our estimates are smaller by an order of magnitude. In our opinion, the 
remarkable recent progress achieved with the simulation of light dynamical quarks and the 
fact that all of the data that pass our quality criteria are consistent with one another fully 
justify this reduction of the uncertainties. 

The lower half of Figure 1 shows that the sum rule results for m s agree with the Nf = 2+1 
data, but Figure 2 indicates that in the case of m u d, the sum rule estimates are on the high 
side. For the most recent sum rule evaluation [121], the origin of this apparent discrepancy 
is readily identified. In that work, the sum rules are used to estimate the value of m s + m U( i- 
In order to disentangle the two terms, the authors take the ratio m s /m u d from the early 
literature on xPT. If the numerical value of this ratio is updated to the range given in (24), 
the value of m s practically stays put, but the outcome for m u d is lowered substantially, so 
that the difference shrinks to about 1 a. In other words, there is no discrepancy to speak of 
between the lattice and sum rule results for the masses of the light quarks - the size of m s 
and m u d is well understood. 

In the ratio m s /m u d, one of the sources of systematic error - the uncertainties in the 
renormalization factors - drops out. Also, we can compare the lattice results with the leading 
order formula of xPT, 

lo M| + +M| -M 2 + ^ 
m ud M 2 + 

which relates the quantity m s /m u d to a ratio of meson masses in QCD. Expressing these in 
terms of the physical masses and the four coefficients introduced in (12)-(14), linearizing the 
result with respect to the corrections and inserting the observed mass values, we obtain 

— = 25.9 - 0.1 e + 1.9 e^o - 0.1 e K o - 1.8 e m . (22) 
m u d 

If the coefficients e, e n o , e K o and e m are set equal to zero, the right hand side reduces to the 
value m s /m u d = 25.9 that follows from Weinberg's leading order formulae for m u /md and 
m s /md [128], in accordance with the fact that these do account for the e.m. interaction at 
leading order, but neglect the mass difference between the charged and neutral pions in QCD. 
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Figure 3: Results for the ratio m s /m u d, in which the renormalization factors cancel. The 
upper part of the figure shows the lattice results listed in Table 3 (full and empty symbols 
correspond to Nf = 2 + 1 and Nf = 2). The vertical bands indicate our estimates (19) and 
(24). The lower part shows results obtained on the basis of xPT or from QCD sum rules, as 
well as the most recent estimate of the Particle Data Group. 



Inserting the estimates (15), the LO prediction becomes 



^ = 25.9(1) . (23) 
m ud 

The quoted uncertainty does not include an estimate for the higher order contributions, 
but only accounts for the error bars in the coefficients, which is dominated by the one in the 
estimate given for e n o . The result shows that, at the accuracy reached in lattice determinations 
of the ratio m s /m u d, the uncertainties due to the electromagnetic corrections is smaller than 
the systematic errors from other sources. 

The lattice results in Table 3 indicate that the corrections generated by the non-leading 
terms of the chiral perturbation series are remarkably small, in the range 5-11%. Despite 
the fact that the SU(3) flavour symmetry breaking effects in the Nambu-Goldstone boson 
masses are very large [M\ ~ 13M|), the mass spectrum of the pseudoscalar octet obeys the 
SU(3)xSU(3) formula (21) very well. 

Our average for m s /m U( i is based on the results of MILC 09A and RBC/UKQCD 10A - 
the value quoted by HPQCD 10 does not represent independent information as it relies on the 
result for m s /m u d obtained by the MILC collaboration. Adding the statistical and systematic 
errors in quadrature, we obtain a good fit with m s /m uc i = 27.39(23) and x 2 = 0-2 for 2 data 
points and one free parameter. The fit is dominated by MILC 09A. Since the errors associated 
with renormalization drop out in the ratio, the uncertainties are even smaller than in the case 
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of the quark masses themselves: the above number for m s /m uc i amounts to an accuracy of 
0.8%. Again, the results of the BMW collaboration are not included in the fit because they are 
not yet published. The result reported in BMW 10A, 10B, m s /m uc i = 27.53(22), is perfectly 
consistent with the above average and reaches comparable accuracy. 

At this level of precision, the uncertainties in the electromagnetic and strong isospin 
breaking corrections are not completely negligible. The error estimate in the LO result (23) 
indicates the expected order of magnitude. The uncertainties in m s and m u d associated with 
the heavy sea quarks cancel at least partly. In view of this, we are convinced that our final 
estimate, 
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N f = 2 + 1: —2- = 27.4(4) , (24) 
m ud 

is on the conservative side - it corresponds to an accuracy of 1.5%. 

The lattice results show that the LO prediction of in (23) receives only small cor- 
rections from higher orders of the chiral expansion: according to (24), these generate a shift 
of 5.8 ± 1.5%. Our estimate does therefore not represent a very sharp determination of the 
higher order contributions. Note however, that in comparison with the value m s /m U( i = 26(4) 
quoted by the Particle Data Group [122], the uncertainty is reduced by a factor of 10. 

The ratio m s /m U( i can also be extracted from the masses of the neutral Nambu-Goldstone 
bosons: neglecting effects of order (m u — m^) 2 also here, the leading order formula reads 

m s /m u d = ^M^/M^ — 5. Numerically, this gives m s /m uc [ = 24.2. The relation has the 
advantage that the e.m. corrections are expected to be much smaller here, but it is more 
difficult to calculate the 77-mass on the lattice. The comparison with (24) shows that, in this 
case, the contributions of NLO are somewhat larger: 13 ± 2%. 



3.4 Lattice determination of m u and rrid 

The determination of m u and rrid separately requires additional input. MILC 09 A [59] uses the 
mass difference between K° and K + , from which they subtract electromagnetic effects, using 
Dashen's theorem with corrections, as discussed in section 3.1. The up and down sea quarks 
remain degenerate in their calculation, fixed to the value of m u d obtained from M n o . Instead 
of subtracting electromagnetic effects using phenomenology, RBC 07 [73] and Blum 10 [74] 
actually include a quenched electromagnetic field in their calculation. This means that their 
results include corrections to Dashen's theorem, albeit only in the presence of quenched 
electromagnetism. Since the up and down quarks in the sea are treated as degenerate, very 
small isospin corrections are neglected, as in MILC's calculation. The results are summarized 
in Table 4. 

In order to discuss these results, we consider the LO formula 



m„ lo 



Ml. - ML + M 2 + 

- A + A" 7r+ (OK\ 

m d ~ M| - M 2 R+ + M 2 + ' 

Using equations (12)— (14) to express the meson masses in QCD in terms of the physical ones 
and linearizing in the corrections, this relation takes the form 

TYL lo 

— H = 0.558 - 0.084 e - 0.02 e n o + 0.11 e m . (26) 
Inserting the estimates (15) and adding errors in quadrature, the LO prediction becomes 

^ = 0.50(4) . (27) 
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Collaboration 



Ref. 



<? / cf & £ f 



m d 



m u /m d 



HPQCD 10* 
BMW 10A, 10B+ 
Blum 10* 
MILC 09A 
MILC 09 

MILC 04, HPQCD/ 
MILC/UKQCD 04 



[104] 
[65, 105] 
[74] 
[59] 
[6] 



★ 
* 
■ 

★ 
★ 



* 
★ 

* 
* 



★ - 2.01(14) 

★ a 2.15(03)(10) 

★ - 2.24(10) (34) 

★ - 1.96(0)(6)(10)(12) 
- 1.9(0)(1)(1)(1) 



[77, 111] A 



4.77(15) 
4.79(07)(12) 
4.65(15)(32) 
4.53(1)(8)(23)(12) 
4.6(0) (2) (2) (1) 



0.448(06) (29) 
0.4818(96)(860) 
0.432(1)(9)(0)(39) 
0.42(0) (1)(0) (4) 



1.7(0) (1) (2) (2) 3.9(0) (1) (4) (2) 0.43(0) (1)(0) (8) 



RBC 07* [73] A ■ ■ ★ ★ - 3.02(27)(19) 5.49(20)(34) 0.550(31) 



* Values obtained by combining the HPQCD 10 result for m s with the MILC 09 results for m a /m u d and 
m u /m d - 

+ The fermion action used is tree-level improved. 

* The calculation includes quenched e.m. effects. 

a The masses are renormalized and run non-perturbatively up to a scale of 4 GeV in the Nf = 3 RI/MOM 
scheme. In this scheme, non-perturbative and N 3 LO running for the quark masses are shown to agree 
from 6 GeV down to 3 GeV to better than 1% [105]. 

Table 4: Lattice results for m u , and m u /md, together with the colour coding of the 
calculation used to obtain them. The upper part of the table lists results obtained with 
Nf = 2 + 1, while the lowest row represents a calculation with Nf = 2. 



Again, the quoted error exclusively accounts for the errors attached to the estimates (12) - 
contributions of non-leading order are ignored. The uncertainty in the leading order prediction 
is dominated by the one in the coefficient e, which specifies the difference between the mass 
splittings generated by the e.m. interaction in the kaon and pion multiplets. Since the input 
e = 1.2(5) used in MILC 09A [59] differs from the results e = 0.13(4) and e = 0.5(1) obtained 
by RBC 07 [73] and Blum 10 [74], respectively, it does not come as a surprise that the values 
for m u /md are different. Evaluating the relation (26) for these values of e and neglecting the 
self-energies of the neutral particles, we obtain mu/md = 0.46(4), 0.547(3) and 0.52(1) for 
MILC 09A, RBC 07 and Blum 10, respectively. Within errors, these numbers agree with the 
values for m u /md quoted in the last column of Table 4. This indicates that, as in the case of 
m s /m u d, the corrections to the leading order prediction are small. 

The uncertainties in the lattice results discussed so far for the ratio m u /md are dominated 
by those in the e.m. contributions to the meson masses. Given the exploratory nature of the 
RBC 07 calculation, we do not draw a conclusion concerning the values of m u and obtained 
with Nf = 2. As discussed in section 3.3.2, the results of Blum 10 do not pass our selection 
criteria either. We therefore resort to the phenomenological estimates of the electromagnetic 
self-energies discussed in section 3.1 and use the MILC data for Nf = 2 + 1 to assess the 
size of the corrections from higher orders of the chiral expansion. With this input, we can 
calculate the mass ratio of the two lightest quarks, with the result m u /md = 0.47(4). The 
two individual masses can then be worked out from the estimate (20) for their mean. The 
result reads 
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N f = 2 + 1: m u = 2.19(15) MeV , m d = 4.67(20) MeV , — = 0.47(4) . (28) 
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These estimates have uncertainties of order 7%, 4% and 9%, respectively. 

In contrast to the situation with the ratio m s /m u d, the error in the result for m u /rrid is 
totally dominated by the uncertainties in the input used for the electromagnetic corrections. 
In fact, the comparison of equations (27) and (28) indicates that more than half of the 
difference between the prediction m u /md = 0.558 obtained from Weinberg's mass formulae 
[128] and the result for m u jm < j L obtained on the lattice stems from electromagnetism - the 
higher orders in the chiral perturbation series only generate a small correction. 

To determine m u /md, BMW 10A, 10B [65, 105] follow a slightly different strategy. They 
obtain this ratio from their result for m s /m u d combined with a phenomenological determina- 
tion of the isospin breaking quark mass ratio Q = 22.3(8) from 77 — > 3ir decays [70]. In fact, 
as discussed in section 3.5, the central value of the e.m. parameter e in (15) is taken from the 
same source. Their results for m u /rrid, m u and md are thus correlated with our estimates in 
(28) and are compatible with them, but are more precise. 

In view of the fact that a massless up quark would solve the strong CP-problem, many 
authors have considered this an attractive possibility, but the results presented above exclude 
this possibility: the value of m u in (28) differs from zero by 15 standard deviations. We 
conclude that nature solves the strong CP-problem differently. This conclusion relies on 
lattice calculations of kaon masses and on the phenomenological estimates of the e.m. self- 
energies discussed in section 3.1. The uncertainties therein currently represent the limiting 
factor in determinations of m u and rrid- As demonstrated in [72-74], lattice methods can be 
used to calculate the e.m. self-energies. Further progress on the determination of the light 
quark masses hinges on an improved understanding of the e.m. effects. 



3.5 Estimates for R and Q 

The quark mass ratios 

m d -m u m d~ m i 

compare SU(3) breaking with isospin breaking. The quantity Q is of particular interest 
because of a low energy theorem [132], which relates it to a ratio of meson masses, 

Ql ^%- ^ K ~^ ' M^\{Ml + + Ml ), Ml^M^+Mlo). (30) 

M n M K° K+ 

Chiral symmetry implies that the expansion of Q 2 M in powers of the quark masses (i) starts 
with Q 2 and (ii) does not receive any contributions at NLO: 

NLO 

Qm = Q- (31) 

Inserting the estimates for the mass ratios m s /m u d, and m u /rrid given in equations (24) 
and (28), respectively, we obtain 

# = 36.6(3.8), Q = 22.8(1.2). (32) 

As is the case with the ratio m u /md, the errors are dominated by the uncertainties in the 
electromagnetic corrections. To investigate the sensitivity to the latter, we use to ex- 
press the quark mass ratios in terms of the pion and kaon masses in QCD and then again 
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use equations (12)-(14) to relate the QCD masses to the physical ones. Linearizing in the 
corrections, this leads to 



R = 43.9- 10.8 e + 0.2 e^o -0.2e x o - 10.6 e. 



(33) 



NLO 



Q 



24.3 - 3.0 e + 0.9 e^o - 0.1 e K o + 2.1 e 



m 



While the first relation only holds to LO of the chiral perturbation series, the second remains 
valid at NLO, on account of the low energy theorem mentioned above. The first terms on 
the right hand side represent the values of R and Q obtained with the Weinberg leading 
order formulae for the quark mass ratios [128]. Inserting the estimates (15), we find that the 

LO NLO 

e.m. corrections lower the Weinberg values to R = 36.7(7.6) and Q = 22.3(2.1), respectively. 
The comparison with the full results quoted above shows that there is little room for the higher 
order terms in the chiral expansion. 

As mentioned in section 3.1, there is a phenomenological determination of Q based on 
the decay 77 — > 37r [133, 134]. The key point is that the transition rj — > 3tt violates isospin 
conservation. The dominating contribution to the transition amplitude stems from the mass 
difference m u — rrid- At NLO of xPT, the QCD part of the amplitude can be expressed in a 
parameter free manner in terms of Q. It is well-known that the electromagnetic contributions 
to the transition amplitude are suppressed (a thorough recent analysis is given in [135]). This 
implies that the result for Q is less sensitive to the electromagnetic uncertainties than the 
value obtained from the masses of the Nambu-Goldstone bosons. For a recent update of this 
determination and for further references to the literature, we refer to [136]. Using dispersion 
theory to pin down the momentum dependence of the amplitude, the observed decay rate 
implies Q = 22.3(8) (since the uncertainty quoted in [136] does not include an estimate for 
all sources of error, we have retained the error estimate given in [129], which is twice as 
large). The formulae for the corrections of NNLO are available also in this case [137] - the 
poor knowledge of the effective coupling constants, particularly of those that are relevant 
for the dependence on the quark masses, is currently the limiting factor encountered in the 
application of these formulae. 

As was to be expected, the central value of Q obtained from 77-decay agrees exactly with 
the central value obtained from the low energy theorem: we have used that theorem to 
estimate the coefficient e, which dominates the e.m. corrections. Using the numbers for e m , 
e n o and e K o in (15) and adding the corresponding uncertainties in quadrature to those in the 
phenomenological result for Q, we obtain 



The estimate (15) for the size of the coefficient e is taken from here, except that the error 
bar is stretched, to account for the fact that the above relation is valid only up to NNLO 
contributions. 

The MILC data offer a welcome check on the size of the neglected higher orders: with the 
input e = 1.2(5) and e n o = e K o = e m = used in the MILC analysis, equation (34) implies 
Qm = 21.3(1.0). On the other hand, inserting the MILC results for the quark mass ratios 
msl m ud and m u /rrid (see Tables 3 and 4) in the definition of Q, equation (29), and ignoring 
a possible correlation between the results for the two ratios, we get Q = 21.7(1.1). For MILC 
09 A, the central values of Qm and Q thus agree to an accuracy of 2%, indicating that the 
corrections of NNLO or higher are smaller than the uncertainty of the result. 



NLO 

e = 



0.70(28) . 



(34) 
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In contrast to this, the picture which follows from the numbers obtained with the Nf = 2 
simulation of RBC 07 [73] is difficult to understand. Repeating the above calculation with 
the estimate e = 0.13(4) that corresponds to the kaon self-energy difference obtained in this 
reference and again neglecting the other terms in equation (34), we get Qm = 23.9(1), while 
the value of Q that follows from the quark mass ratios m s /m ud and m u /rrid given in the 
same reference, implies Q = 26.1(1.2). A difference of this size calls for large contributions 
of NNLO, in conflict with the data analysis, which assumes that it is meaningful to truncate 
the chiral series at NLO. 

The recent extension of this calculation to Nf = 2 + 1 dynamical flavours described in [74] 
appears to solve the puzzle. The values obtained for the kaon and pion self-energies obtained 
in this reference imply e = 0.5(1), e m = 0.04(5). Inserting this in equation (34) and neglecting 
the self-energies of the neutral particles, we obtain Qm = 22.8(3), to be compared with the 
value Q = 23.9(3.4) that follows from the quark mass ratios given in the same reference. This 
implies that the NNLO corrections to the low energy theorem (31) are too small to be seen 
at the accuracy of this calculation: Qm/Q = 0.96(12). 

The determination of R and Q in BMW 10A, 10B [65, 105] relies on the phenomenological 
value of Q discussed above [70] - their results do not bring new information on Q. The 
value of R which follows from their result for the ratio m s /m uc i and the input used for Q is 
R = 34.9(2.5), slightly different from our estimate in (32), but perfectly compatible with it. 7 

Our final results for the masses m u , m d , m s and the mass ratios m u /m d , m s /m u d, R, Q 
are collected in the table below. 

Nf m u m d m s m ud m u /m d m s /m ud R Q 

2+1 2.19(15) 4.67(20) 94(3) 3.43(11) 0.47(4) 27.4(4) 36.6(3.8) 22.8(1.2) 
2 95(6) 3.6(2) 27.3(9) 



Table 5: Our estimates for the masses of the three lightest quarks, and related ratios (the 
masses refer to the MS scheme at running scale = 2 GeV for Nf = 3, the numerical values 
are given in MeV units). 



7 Repeating the calculation with the same input for Q, but replacing the BMW result for m B /m u d by our 
estimate in equation (24), we obtain R = 35.0(2.6). 
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4 \V ud \ and \V US \ 

4.1 Experimental information concerning \V u d\, \V US \, /+(0) and fx/fw 

The following review relies on the fact that precision data on kaon decays very accurately 
determine the product |V^ s |/+(0) and the ratio \V us fK\/\V U( ifn\ [138]: 

\V us \f + (0) = 0.2163(5), 

V u d and V us are elements of the Cabibbo-Kobayashi-Maskawa matrix and f+(t) represents 
one of the form factors relevant for the semileptonic decay K° — > -k~£v, which depends on 
the momentum transfer t between the two mesons. What matters here is the value at t = 0: 
/+(0) = f+ n {t)\ t _^ - The pion and kaon decay constants are defined by 8 

(0|d7 M 75u|vr + (p)) = ip^U , (0| 57^75 u\K + (p)} =ip^K- 

In this normalization, f n ~ 130 MeV, fx — 155 MeV. 

The measurement of \V u d\ based on superallowed nuclear f3 transitions has now become 
remarkably precise. The result of the recent update of Hardy and Towner [50], which is based 
on 20 different superallowed transitions, reads 9 

\Vud\ = 0.97425(22) . (36) 

The matrix element \V US \ can be measured in r decays [148-151]. Separating the inclusive 
decay r —> hadrons + v into non-strange and strange final states, Gamiz et al. [152] obtain 

\V US \ = 0.2165(26^(5)^. (37) 

Maltman et al. [153] arrive at very similar values. 

As recently pointed out by Maltman [150], the theoretical uncertainties of the analysis 
which underlies the result (37) can be reduced by invoking the experimental information about 
the spectral function of the electromagnetic current, but the experimental uncertainties then 
play a more important role. Applying this method, the outcome for \V US \ reads [154] 

\V US \ = 0.2208(39) . (38) 

In principle, r decay offers a clean measurement of \V US \, but a number of open issues 
yet remain to be clarified. In particular, the measured exclusive decay rates for r — > ttv and 
t — > Ku are below the Standard Model predictions, which determine these rates in terms 

8 The pion decay constant represents a QCD-matrix-element - in the full Standard Model, the one-pion 
state is not a meaningful notion: the correlation function of the charged axial current does not have a pole at 
p 2 — M 2 + , but a branch cut extending from M 2 + to 00. The analytic properties of the correlation function and 
the problems encountered in the determination of f n are thoroughly discussed in [139]. The "experimental" 
value of fn depends on the convention used when splitting the sum £qcd + £qed into two parts (compare 
section 3.1). The lattice determinations of do not yet reach the accuracy where this is of significance, but 
at the precision claimed by the Particle Data Group [140], the numerical value does depend on the convention 
used [67-69, 139]. A recent analysis of the electromagnetic and strong isospin breaking effects in the leptonic 
decays of the pseudoscalar mesons is given in [141]. 

It is not a trivial matter to perform the data analysis at this precision. In particular, isospin-breaking 
effects need to be properly accounted for [142-146]. For a review of recent work on this issue, we refer to [147]. 



= 0.2758(5). (35) 

VudJn 
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RBC 06 


[163] 
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★ 


0.968(9)(6) 


JLQCD 05 


[164] 


2 


C 


■ 


■ 


★ 


0.967(6), 0.952(6) 



Table 6: Colour code for the data on /+(0). 



of the same matrix elements [V^/vi-l and |T4s/x| that govern the leptonic decays tt — >• \iv 
and K — )• [iv (violation of r/fi universality). On the other hand, the value obtained for 
|Kis/i<r|/|^ud/7r| from the ratio of the tau decay rates is perfectly consistent with the one 
from the leptonic n and A'-decays quoted in (35), but this does not shed any light on the 
value of \V US \. It is important to pursue the measurement of modes that have previously been 
studied only with low statistics (especially the large mode K 7r°7r _ ), as well as those of higher 
multiplicity (final states with more than two pions, for instance). The recent developments 
on the theoretical side [155-157] also need to be pursued. The most interesting possibility 
is that t decay involves new physics, but more work is required before r decay becomes a 
competitive source of information about the CKM matrix elements. 

4.2 Lattice results for /+(0) and f K j '/„. 

The traditional way of determining \ V US \ relies on using theory for the value of /+(0), invoking 
the Ademollo-Gatto theorem [165]. Since this theorem only holds to leading order of the 
expansion in powers of m u , and m s , theoretical models are used to estimate the corrections. 
Lattice methods have now reached the stage where quantities like /+(0) or fx/ f-w can be 
determined to good accuracy. As a consequence, the uncertainties inherent in the theoretical 
estimates for the higher order effects in the value of /+(0) do not represent a limiting factor 
any more and we shall therefore not invoke those estimates. Also, we will use the experimental 
results based on nuclear f3 decay and r decay exclusively for comparison - the main aim of 
the present review is to assess the information gathered with lattice methods and to use it for 
testing the consistency of the SM and its potential to provide constraints for its extensions. 

The database underlying the present review of /+ (0) and fx / fn is listed in Tables 6 and 
7. The properties of the lattice data play a crucial role for the conclusions to be drawn from 
these results: range of M n , size of LM n , continuum extrapolation, extrapolation in the quark 
masses, finite size effects, etc. The key features of the various data sets are characterized by 
means of the colour code specified in section 2. 
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The quantity /+(0) represents a matrix element of a strangeness changing null plane 
charge, /+(0) = (K\Q us \ir). The vector charges obey the commutation relations of the Lie 
algebra of SU(3), in particular [Q US ,Q SU ] = Q uu ~ ss . This relation implies the sum rule 

\(K\Q us \n)\ 2 — ^\ \(K\Q su \n)\ 2 = 1. Since the contribution from the one-pion intermedi- 
ate state to the first sum is given by /+(0) 2 , the relation amounts to an exact representation 
for this quantity [166]: 

/ + (0) 2 = 1-1] \{K\Q us \n)\ 2 + £ \(K\Q su \n)\ 2 . (39) 

n=£ir n 

While the first sum on the right extends over non-strange intermediate states, the second runs 
over exotic states with strangeness ±2 and is expected to be small compared to the first. 

The expansion of /+(0) in powers of m u , and m s starts with /+(0) = 1 + fa + fa + . . . 
[91]. Since all of the low energy constants occurring in fa can be expressed in terms of M w , 
Mk-, M?) and fa [167], the correction of NLO is known. In the language of the sum rule (39), 
fa stems from non-strange intermediate states with three mesons. Like all other non-exotic 
intermediate states, it lowers the value of /+(0): fa = —0.023. The corresponding expressions 
have also been derived in quenched or partially quenched chiral perturbation theory [168]. 
At the same order in the SU(2) expansion [169], /+(0) is parametrized in terms of M n and 
two a priori unknown parameters. The latter can be determined from the dependence of the 
lattice results on the masses of the quarks. Note that any calculation that relies on the xPT 
formula for fa is subject to the uncertainties inherent in NLO results: instead of using the 
physical value of the pion decay constant fa, one may, for instance, work with the constant 
fa that occurs in the effective Lagrangian and represents the value of fa in the chiral limit. 
Although trading fa for fa in the expression for the NLO term affects the result only at 
NNLO, it may make a significant numerical difference in calculations where the latter are 
not explicitly accounted for (the lattice results concerning the value of the ratio fa/ fa are 
reviewed in section 5.2). 

The lattice results shown in the left panel of Figure 4 indicate that the higher order 
contributions A/ = /+(0) — 1 — fa are negative and thus amplify the effect generated by fa. 
This confirms the expectation that the exotic contributions are small. The entries in the lower 
part of the left panel represent various model estimates for fa. In [179] the symmetry breaking 
effects are estimated in the framework of the quark model. The more recent calculations are 
more sophisticated, as they make use of the known explicit expression for the form 
factors to NNLO of xPT [180, 181]. The corresponding formula for fa accounts for the chiral 
logarithms occurring at NNLO and is not subject to the ambiguity mentioned above. The 
numerical result, however, depends on the model used to estimate the low energy constants 
occurring in fa [181-184]. The figure indicates that the most recent numbers obtained in this 
way correspond to a positive rather than a negative value for A/. 

A lattice determination of the low energy constants that occur in fa would be very useful 
to identify the origin of the problem. Lattice results concerning the dependence of the scalar 
form factor on the momentum transfer would also very useful, in particular they might help 
sorting out the discrepancies in the experimental results for the slope of this form factor 
[185-188]. Note that when analyzing lattice data on the basis of the xPT formulae, the loop 
integrals occurring therein need to be evaluated for the masses at which the Nambu-Goldstone 
bosons propagate in a given lattice simulation. 10 

10 Fortran programs for the numerical evaluation of the form factor representation in [181] are available on 
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Table 7: Colour code for the data on fx/ fir- 
4.3 Direct determination of /+(0) and fx] fir 

Figure 4 shows that the lattice results for /+(0) and fx/fir obtained by the various collab- 
orations which do a comprehensive error analysis are consistent with each other. We now 
proceed to form the corresponding averages, separately for the data with Nf = 2 + 1 and 
Nf = 2 dynamical flavours. As will be discussed in detail in section 4.5, CKM unitarity and 
experiment correlate /+(0) with fx/In- Indeed, the lattice data are consistent also with this 
correlation. First, however, we analyze the lattice information available separately for /+(0) 
and for fx/fn- We refer to results obtained in this way as "direct" determinations. 

The colour code in Table 5 shows that for /+(0), presently only the results of the ETM 
collaboration with Nf = 2 dynamical flavours of fermions are without a red tag. In the 
following, we rely on the published work ETM 09A [161]. The table contains two data 
sets with Nf = 2 + 1, but since one represents an update of the other, we only discuss 
RBC/UKQCD 10. These data carry a red tag, because a systematic study of cut-off effects 
is missing. These two computations of /+(0) are the most advanced ones and we quote their 
results: 



/+(0) = 0.9599(34)(lj£)(14), 
/+(0) = 0.9560(57)(62), 



(direct, Nf = 2 + 1) 
(direct, Nf = 2). 



(40) 



The first error in the result quoted for Nf = 2 + 1 is statistical, the second is due to the 
uncertainties in the chiral extrapolation of the lattice data and the third is an estimate of 
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* Estimates obtained from an analysis of the lattice data within the Standard Model, see section 4.5. 

Figure 4: Comparison of lattice results (red squares) for /+(0) and fx/f-n with various model 
estimates based on (blue triangles). Full and empty squares represent simulations with 
Nf = 2 + 1 and Nf = 2, respectively. The vertical bands indicate our estimates.* 



potential discretization effects. Flavour symmetry implies that if m u d is set equal to m s , the 
lattice data yield /+(0) = 1, irrespective of the lattice spacing or the size of the box and for 
any value of m s . Cut-off effects can therefore only affect the difference 1 — /+(0), which turns 
out to be about 0.04. Indeed, the estimate provided by RBC/UKQCD 10 for the uncertainties 
due to discretization effects shows that these are sub-dominant: inflating the corresponding 
error by a factor of up to 2 barely affects the net systematic uncertainty. 

In the result quoted for Nf = 2, the brackets indicate the statistical and systematic errors, 
respectively. The ETM collaboration provides a more comprehensive study of the systematics 
by presenting results for three lattice spacings [189] and simulating at lighter pion masses 
(down to M n = 260 MeV). This allows to better constrain the chiral extrapolation, using both 
SU(3) [167] and SU(2) [169] chiral perturbation theory. Moreover, a rough estimate for the 
size of the effects due to quenching the strange quark is given, based on the comparison of the 
result for Nf = 2 dynamical quark flavours [177] with the one in the quenched approximation, 
obtained earlier by the SPQcdR collaboration [190]. 

The quality criteria laid out in section 2 require a systematic study of lattice artifacts. 
As indicated by the colour code in Table 6, the errors due to the continuum extrapolation 
yet need to be investigated in more detail for the data with Nf = 2 + 1, while for Nf = 2, 
where the quoted uncertainties are larger, these errors are under somewhat better control. 
The value 

/+(0) = 0.956(8) (our estimate, direct) (41) 

covers both results in equation (40). In our opinion, it represents a conservative estimate for 
the range permitted by the presently available direct determinations of /+(0) in lattice QCD, 
not only for Nf = 2, but also for Nf = 2 + 1. 

For fx/ fn, Table 7 contains several simulations with Nf = 2 + 1 dynamical quark flavours. 
The latest update of the MILC program is reported in MILC 10 [171]. We use the results 
quoted there when forming averages. Three further data sets meet the criteria formulated in 
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Figure 5: The plot compares the information for \V u d\, \V US \ obtained on the lattice with 
the experimental result extracted from nuclear /3 transitions. The dotted arc indicates the 
correlation between \V u d\ and \V US \ that follows if the three- flavour CKM- matrix is unitary. 

the introduction: BMW 10 [41] and HPQCD/UKQCD 07 [175] with Nf = 2 + 1 and ETM 
09 [177] with Nf = 2 dynamical flavours. We ignore possible correlations due to the fact that 
MILC 10 and HPQCD/UKQCD 07 have partly used the same set of gauge configurations and 
apply the procedure outlined in section 2.2 to the three sets with Nf = 2 + 1. The resulting fit 
is of good quality, with fx/ fw = 1-193(4) and x 2 = 0-4 f° r 3 data points and 1 free parameter. 
The systematic errors of the individual data sets are larger: 0.005, 0.007 and 0.006 for MILC 
10, HPQCD/UKQCD 07 and BMW 10, respectively. Following the prescription of section 
2.2, we replace the error by the smallest one of these numbers. Together with the ETM 09 
result for Nf = 2, our estimates thus read 

fx/ U = 1-193(5) , (direct, Nf = 2 + 1), (42) 

f K /U = 1-210(6) (17) , (direct, N f = 2). 

It is instructive to convert the above results for /+(0) and fx/fir into a corresponding range 
for the CKM matrix elements V ud and V us , using the relations (35). Consider first the results 
for Nf = 2 + 1. The range for /+(0) in (40) is mapped into the interval V us = 0.2255(14), 
depicted as a horizontal gray band in Figure 5, while the one for fx/fn in (42) is converted 
into V us /V ud = 0.2312(11), shown as a green band. The red curve is the intersection of these 
two bands. More precisely, it represents the 68% likelihood contour, obtained by treating 
the above two results as independent measurements. A Gaussian in /+(0) corresponds to a 
Gaussian in the variable 1/V US . Since the width is small, the distribution in the variable V us is 
also approximately Gaussian. The corresponding likelihood function is given by x\ = (Vus ~ 
0.2255) 2 /0.0014 2 . Likewise, a Gaussian in fx/fn is mapped into an approximately Gaussian 
distribution of the variable V us /V ud , with x\ = (V us /V ud - 0.2312) 2 /0.0011 2 . Expressed 
in terms of the CKM matrix elements, the Nf = 2 + 1 results for /+(0) and fx/ fir are 
thus characterized by the likelihood function x 2 = x\ + Xh- r ^ ne minimum occurs at the 
intersection of the centers of the two bands, where x 2 vanishes. The contour shown is the line 
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Collaboration 


Ref. 


N f 


from 


\Vus\ 


MILC 10 


[171] 

L J 


2 + 1 


fic/fn 


0.2245(5)(+l 2 ) 


RBC/UKQCD 10A 


[106] 


2 + 1 


Ik/U 


0.2233(13) (44) 


RBC/UKQCD 10 


[158] 


2 + 1 


/+(0) 


0.2253(10) (^ 2 ) 


BMW 10 


[41] 


2 + 1 


fic/fn 


0.2254(13)(11) 


HPQCD/UKQCD 07 


[175] 


2 + 1 


Ik/U 


0.2260(5)(13) 


ETM 09 


[177] 


2 


Ik/U 


0.2222(11)(31) 


ETM 09A 


[161] 


2 


/+(0) 


0.2263(14)(15) 



Table 8: Values of \V US \ obtained from lattice determinations of /+(0) or fx/ fn with CKM 
unitarity. The first (second) number in brackets represents the statistical (systematic) error. 



where x 2 differs from the minimum by unity. Values of V us , V u d in the region enclosed by this 
contour are consistent with the lattice data for Nf = 2 + 1, within one standard deviation. In 
particular, the plot shows that the nuclear j3 decay result for V u d is perfectly consistent with 
these data. 

Repeating the exercise with the Nf = 2 results leads to the dashed ellipse. The figure 
thus indicates that the data are consistent within errors. We conclude that, at the accuracy 
reached up to now for the quantities /+(0), fx/ fir, Vus and V u d, the distortions generated by 
the quenching of the strange quark are too small to be visible. 

4.4 Testing the Standard Model 

In the Standard Model, the CKM matrix is unitary. In particular, the elements of the first 
row obey 

|K| 2 = \V ud \ 2 + \V US \ 2 + \V ub \ 2 = 1 . (43) 

The tiny contribution from \V u b\ is known much better than needed in the present context: 
\V u b\ = 3.89(44) • 10~ 3 [191]. In the following, we first discuss the evidence for the validity of 
the relation (43) and only then use it to analyze the lattice data within the Standard Model. 

In Figure 5, the correlation between V u d and V us imposed by the unitarity of the CKM 
matrix is indicated by a dotted arc (more precisely, in view of the uncertainty in V u b, the 
correlation corresponds to a band of finite width, but the effect is too small to be seen here). 
The plot shows that the data for Nf = 2 + 1 are perfectly consistent with this constraint. 
Numerically, the outcome for the sum of the squares of the first row of the CKM matrix reads 
|^u| 2 = 1.002(15). The Standard Model thus passes a nontrivial test that exclusively involves 
lattice data and well-established kaon decay branching ratios. Combining the lattice results 
for /+(0) and Jk/ fn in (40) and (42) with the /3-decay value of \V u d\ quoted in (36), the test 
sharpens considerably: the lattice result for /+(0) leads to |V^| 2 = 1.0000(7), while the one 
for fx/ fn implies \ V U \ 2 = 0.9999(6), thus confirming CKM unitarity at the permille level. 

Repeating the analysis for Nf = 2, we find |^u| 2 = 1.037(36) with the lattice data 
alone. The number is somewhat larger than 1, in accordance with the fact that the dot- 
ted curve passes just outside the blue contour, but the difference barely exceeds one stan- 
dard deviation. Moreover, it only concerns the comparison of the Nf = 2 results for /+(0) 
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Estimates obtained from an analysis of the lattice data within the Standard Model, see section 4.5. 



Figure 6: Results for \V US \ and \V u d\ that follow from the lattice data for /+(0) (red triangles) 
and fx/ fn (red squares), on the basis of the assumption that the CKM matrix is unitary. 
The black squares and the bands represent our estimates, obtained by combining these two 
different ways of measuring \V US \ and \V u d\ on a lattice. For comparison, the figure also 
indicates the results obtained if the data on nuclear (3 decay and r decay are analyzed within 
the Standard Model. 



with those for fx/In- Taken by themselves, these results are perfectly consistent with the 
value of \V u d\ found in nuclear /3 decay: combining this value with the data on /+(0) yields 
\V U \ 2 = 1.0004(10), combining it with the data on f K /U S ives IK| 2 = 0.9985(16). 

Note that the above tests also offer a check of the basic hypothesis that underlies our 
analysis: we are assuming that the weak interaction between the quarks and the leptons 
is governed by the same Fermi constant as the one that determines the strength of the 
weak interaction among the leptons and determines the lifetime of the muon. In certain 
modifications of the Standard Model, this is not the case. In those models it need not be true 
that the rates of the decays rr — > £v, K — > Iv and K — > ttIv can be used to determine the 
matrix elements 1^^/^!, and |V^ S /+(0)|, respectively and that \V u d\ can be measured 

in nuclear f3 decay. The fact that the lattice data are consistent with unitarity and with 
the value of \V u d\ found in nuclear (3 decay indirectly also checks the equality of the Fermi 
constants. 
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4.5 Analysis within the Standard Model 

The Standard Model implies that the CKM matrix is unitary. The precise experimental 
constraints quoted in (35) and the unitarity condition (43) then reduce the four quantities 
I Kid |) \ Vus\i /+(0)j Ik/ fn to a single unknown: any one of these determines the other three 
within narrow uncertainties. 

Figure 6 shows that the results obtained for |V^ S | and \ V u d\ from the data on fx/fn (red 
squares) are quite consistent with the determinations via /+(0) (red triangles). In order to 
calculate the corresponding average values, we restrict ourselves to those determinations that 
we have considered best in section 4.3. The corresponding results for \ V US \ are listed in Table 
8 (the noise in the experimental numbers used to convert the values of /+(0) and fx/fn into 
values for \V US \ is included in the statistical error). 

We consider the fact that the results from the four Nt = 2 + 1 data sets RBC/UKQCD 
10 [158], BMW 10 [41], MILC 10 [171] and HPQCD/UKQCD 07 [175] are consistent with 
each other to be an important reliability test of the lattice work. Treating the four sets as 
independent measurements, and applying the standard averaging procedure, we obtain a fit 
of good quality, with \V US \ = 0.2253(6) and x 2 = 0.6 for 4 data points and 1 free parameter. 
The results listed in the upper half of Table 8 show, however, that the standard procedure 
underestimates the systematic uncertainties also in this case. Applying the prescription of 
section 2.2, we arrive at a somewhat larger error: |V^ S | = 0.2253(9). This result is indicated 
on the left hand side of Fig. 6 by the narrow vertical band. The broader band shows the 
corresponding value for Nf = 2 (standard error analysis, \V US \ = 0.2253(17), with X 2 = 1-2 
for 2 data points and 1 free parameter, S = 1.09). The figure shows that the result obtained 
for the data with Nf = 2 is perfectly consistent with the one found for Nf = 2 + 1. 

Alternatively, we can solve the relations for \V u d\ instead of \V US \. Again, the result 
\V u d\ = 0.97428(21) which follows from the lattice data with Nf = 2 + 1 is perfectly consistent 
with the value \V u d\ = 0.97433(42) obtained from those with Nf = 2. The reduction of the 
uncertainties in the result for \V UC [\ due to CKM unitarity is to be expected from Figure 5: 
the unitarity condition reduces the region allowed by the lattice results to a nearly vertical 
interval. 

Next, we determine the value of /+(0) that follows from the lattice data within the Stan- 
dard Model. Using CKM unitarity to convert the lattice determinations of fx/fn into corre- 
sponding values for /+(0) and then combining these with the direct determinations of /+(0), 
we find /+(0) = 0.9599(38) from the data with A/ = 2 + 1 and /+(0) = 0.9604(75) for Nf = 2. 
The results are shown in the left panel of Fig. 4. 

Finally, we work out the analogous Standard Model fits for fx/ fx, converting the direct 
determinations of /+(0) into corresponding values for fx/ fir an d combining the outcome 
with the direct determinations of that quantity. The results read fx/ fir = 1.1927(50) for 
Nf = 2 + 1 and fx/fn = 1.194(10) for Nf = 2, respectively; they are shown in the right panel 
of Fig. 4. 

The results obtained by analyzing the lattice data in the framework of the Standard Model 
are collected in the upper half of Table 9. In the lower half of this table, we list the analogous 
results, found by working out the consequences of CKM-unitarity for the experimental values 
of \V u d\ and \V US \ obtained from nuclear /3-decay and r-decay, respectively. The comparison 
shows that the lattice result for \V u d\ not only agrees very well with the totally independent 
determination based on nuclear (3 transitions, but is also remarkably precise. On the other 
hand, the values of \ V u d\, /+(0) and fx/fn which follow from the r-decay data if the Standard 
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Ref. 


\v us \ 


\v ud \ 


/+(0) 




N f = 2 + 1 




0.2253(9) 


0.97428(21) 


0.9599(38) 


1.1927(50) 


N f = 2 




0.2251(18) 


0.97433(42) 


0.9604(75) 


1.194(10) 


/3-decay 


[50] 


0.22544(95) 


0.97425(22) 


0.9595(46) 


1.1919(57) 


r-decay 


[152] 


0.2165(26) 


0.9763(6) 


0.999(12) 


1.244(16) 


r-decay 


[154] 


0.2208(39) 


0.9753(9) 


0.980(18) 


1.218(23) 



Table 9: The upper half of the table shows our final results for \V US \, \V u d\, /+(0) and fx/fir, 
which are obtained by analyzing the lattice data within the Standard Model. For comparison, 
the lower half lists the values that follow if the lattice results are replaced by the experimental 
results on nuclear j3 decay and r decay, respectively. 



Model is assumed to be valid, are not in good agreement with the lattice results for these 
quantities. The disagreement is reduced considerably if the analysis of the r data is supple- 
mented with experimental results on electroproduction [154]: the discrepancy then amounts 
to little more than one standard deviation. 



5 Low-energy constants 



In studying the quark mass dependence of QCD observables calculated on the lattice it is good 
practice to invoke xPT. For a given quantity this framework predicts the nonanalytic quark 
mass dependence and it provides symmetry relations among different observables. These 
symmetry relations are best expressed with the help of a set of universal and independent 
low-energy constants (LECs), which appear as coefficients of the polynomial terms in different 
observables. If one expands around the SU(2) chiral limit, in the Chiral Effective Lagrangian 
there are 2 LECs at order p 2 



F = F W 



B 



(0|«u|0) 



F 2 



(44) 



and seven at order p 4 , indicated by £i, % = 1, . . . , 7. In the analysis of the SU(3) chiral limit, 
at order p 2 there are also only two LECs 



Fn = Ftt 



Bn 



(0\uu\0) 



mu,md,m s — >0 



F 2 



(45) 



but ten at order p 4 indicated by the capital letter Li(/j,), i = 1, . . . , 10. These constants are 
independent of the quark masses, 11 but they become scale dependent after renormalization 
(sometimes a superscript r is added). The SU(2) constants are scale independent, since they 
are defined at [i = M % (as indicated by the bar). For the precise definition of these constants, 
their scale dependence etc, we refer the reader to [91, 93]. 



11 More precisely, they are independent of the 2 or 3 light quark masses which are explicitly considered in 
the respective framework. However, all low-energy constants depend on the masses of the remaining quarks 



s,c,b,t or c, b, t in the SU(2) and SU(3) framework, respectively. 
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First of all, lattice calculations can be used to test if chiral symmetry is indeed broken 
as SU(iV/)iXSU(iV/)_R — >S\J(Nf)L+R by measuring non-zero chiral condensates and by veri- 
fying, for instance, the validity of the GMOR relation M 2 oc m close to the chiral limit. If 
the chiral extrapolation of quantities calculated on the lattice is made with the help of xPT, 
not only does one get the observable at the physical value of the quark masses, but one also 
determines the relevant LECs. This is a very important by-product for two reasons: 

1. All the LECs up to order p A (with the exception of B and Bq, since only the product 
of these times the quark masses can be estimated from phenomenology) have been 
either determined by comparison to experiments or estimated theoretically. A lattice 
determination of the better known ones provides therefore a stringent test of the xPT 
approach. 

2. The less well known LECs are those which describe the quark mass dependence of 
observables - these cannot be determined from experiments, and therefore the lattice 
provides for these unique quantitative information. This information is essential for 
improving those xPT predictions in which these LECs play a role. 

We stress that this program is based on the non-obvious assumption that xPT is valid in the 
region of masses used in the lattice simulation under consideration. 

The fact that, at large volume, the finite-size effects, which occur if a system undergoes 
spontaneous symmetry breakdown, are controlled by the Nambu-Goldstone modes, was first 
noted in solid state physics, in connection with magnetic systems [192, 193]. As pointed out 
in [194], in the context of QCD, the thermal properties of such systems can be studied in 
a systematic and model independent manner by means of the corresponding effective field 
theory, provided only the temperature is low enough. While finite volumes are not of physical 
interest in particle physics, lattice simulations are necessarily carried out in a box. As shown 
in [195-197], the ensuing finite size effects can also be studied on the basis of the effective 
theory - xPT in the case of QCD - provided the simulation is close enough to the continuum 
limit, the volume is sufficiently large and the explicit breaking of chiral symmetry generated 
by the quark masses is sufficiently small. Indeed, xPT represents a useful tool also for the 
analysis of the finite size effects in lattice simulations. 

In the following two sections we will summarize the lattice results for the SU(2) and SU(3) 
LECs at order p 2 and p 4 . The 0(p 2 ) constants are determined from the chiral extrapolation 
of masses and decay constants or, alternatively, from a finite-size study of correlators in the 
e-regime or from spectral observables. At order p some LECs influence two-point functions 
while others appear only in three- or four-point functions; the latter need to be determined 
via a chiral extrapolation of form factors or scattering amplitudes. The xPT analysis of the 
phenomenology is nowadays based on 0(p e ) formulae. 12 At this level the number of LECs 
explodes and we will not discuss any of these here. We will, however, discuss how comparing 
different orders and different expansions (x versus £, see below) can help to reliably assess 
the theoretical uncertainties of the LECs determined on the lattice. 

12 Some of the 0(p 6 ) formulae presented below have been derived in an unpublished note by three of us (GC, 
SD and HL) and Jiirg Gasser. We thank him for allowing us to publish them here. 
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5.1 SU(2) Low-Energy Constants 

5.1.1 Quark-mass dependence of pseudoscalar masses and decay constants 

The expansions of M 2 and F n in powers of the quark mass are known to next-to-next-to- 
leading order in the SU(2) chiral effective theory. 13 In the isospin limit, m u = = m, the 
explicit expressions may be written in the form [198] 

Ml = M 2 |l-lxln^ + ^x 2 (ln^|) 2 + ^ M + 0(x 3 )|, (46) 
( A 2 5 „ / A 2x 2 



F n = Fjl + xln^-^ln^) + ,-*.,+(,(.,■• 

The expansion parameter is given by 

M 2 . o „ 2mS 



--(4nFr ^ = 2Bm=^. (47) 

The logarithmic scales A3, A4 are related to the effective coupling constants £3, £4 of the chiral 
Lagrangian at running scale M v : 

in = ln W' n = l,...,7. (48) 

Note that in Eq. (46), the logarithms are evaluated at M 2 , and not at M 2 . The coupling 
constants k{ in Eq. (46) are mass-independent. The scales of the quadratic logarithms can be 
expressed in terms of the coupling constants: 

A 2 m 1 ( A? A 2 A 2 



ln # = 5Tl 281n M^ + 321n A^- 91n ^ + 49 J' <*> 

A 2 If A 2 A 2 A 2 A 2 

InT^ = — |141n— ir + 161n-Jr + 61n-Jr-61n-^ + 23 

M 2 30 V M 2 M 2 M 2 M 2 

By analyzing the quark mass dependence of M^ and F n with Eq. (46), possibly truncated 
at NLO, one can therefore determine 14 the 0(p 2 ) LECs B and F, as well as the 0(p 4 ) LECs 
£3 and £4. The quark condensate in the chiral limit is determined by S = F 2 B. With precise 
enough data at several low pion masses, one could in principle also determine Am, A^ and 
ku-, kp. This is not yet the case. The results for the LO and NLO constants will be presented 
in Sec. 5.1.6. 

Alternatively, one can invert Eq. (46) and express M and F as an expansion in 

M 2 

f=5£$I- (50) 



13 Here and in the following, we stick to the notation used in the papers where the x?T formulae were 
established, i.e. we work with = fn/V2 = 92.2(1) MeV and F K = /k/a/2- The occurrence of different 
normalization conventions is not convenient, but avoiding it by reformulating the formulae in terms of f n , fx 
is not a good way out. Since we are using different symbols, confusion cannot arise. 

14 Notice that one could analyze the quark mass dependence entirely in terms of the parameter M defined 
in Eq. (47) and determine equally well all other LECs. Using the determination of the quark masses described 
in Sec. 3 one can then extract B or E. 
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The corresponding expressions read 



M 2 



F 



M -{ 1+ 5 fl "|-^ 2 ( ln i) 2+ ^' +OK3) 



(51) 



The scales of the 



quadratic logarithms are determined by Ai, . . . , A4: 




(52) 



5.1.2 Two-point correlation functions in the e-regime 

The finite-size effects encountered in lattice calculations can be used to determine some of 
the Low-Energy Constants of QCD. In order to illustrate the method, we focus on the two 
lightest quarks, take the isospin limit m u = rrid = m and consider a box of size L s in the three 
space directions and size Lt in the time direction. If m is sent to zero at fixed box size, chiral 
symmetry is restored. The behaviour of the various observables in the symmetry restoration 
region is controlled by the parameter fj, = mE7, where V = l? a Lt is the volume of the box. 
Up to a sign and a factor of two, the parameter fi represents the minimum of the classical 
action that belongs to the leading order effective Lagrangian of QCD. 

For /1 > 1, the system behaves qualitatively as if the box was infinitely large. In that 
region, the p-expansion, which counts 1/L S , 1/Lt and M as quantities of the same order, is 
adequate. In view of 11 = \F 2 M 2 V, this region includes configurations with ML > 1 , where 
the finite-size effects due to pion loop diagrams are suppressed by the factor e~ ML . 

If [i is comparable to or smaller than 1, however, the chiral perturbation series must be 
reordered. The e-expansion achieves this by counting 1/L S , 1/Lt as quantities of 0(e), while 
the quark mass m is booked as a term of 0(e 4 ). This ensures that the symmetry restoration 
parameter \x represents a term of order 0(e ), so that the manner in which chiral symmetry 
is restored can be worked out. 

As an example, we consider the correlator of the axial charge carried by the two lightest 
quarks, q(x) = {u(x) , d(x)} . The axial current and the pseudoscalar density are given by 



%PT yields explicit finite-size scaling formulae for these quantities [197, 199, 200]. In the 





(54) 
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e-regime, the expansion starts with 



c AA {t) 



F 2 Ll 

2 r 6 



C PP {t) = Z Z L 



ap+ ^ h {£) +0(e4) 



(55) 



where the coefficients a a, &a> ap, bp stand for quantities of 0(e ). They can be expressed 
in terms of the variables L s , Lt and m and exclusively involve the two leading low-energy 
constants F and S. In fact, at leading order, only the combination /i = m £ L^Lt matters, the 
correlators are independent of t and the dependence on \i is fully determined by the structure 
of the groups involved in the spontaneous symmetry breakdown. In the case of SU(2)xSU(2) 
—> SU(2), relevant for QCD in the symmetry restoration region of the two lightest quarks, the 
coefficients can be expressed in terms of Bessel functions. The t-dependence of the correlators 
starts showing up at 0(e 2 ), in the form of a parabola: h\{r) = ~ (r — — ^ . Explicit 

expressions for a a, bA, ap, bp can be found in [197, 199, 200], where some of the correlation 
functions are worked out to NNLO. By matching the finite-size scaling of correlators computed 
on the lattice with these predictions one can extract F and E. 

The fact that the representation of the correlators to NLO is not "contaminated" by higher 
order unknown LECs, makes the e-regime potentially convenient for a clean extraction of the 
LO couplings. The determination of the LECs is then affected by different systematic un- 
certainties with respect to the standard case; simulations in this regime yield complementary 
information which can serve as a valuable cross-check to get a comprehensive picture of the 
low-energy properties of QCD. 

The effective theory can also be used to study the distribution of the topological charge in 
QCD [201] and the various quantities of interest may be defined for a fixed value of this charge. 
The expectation values and correlation functions then not only depend on the symmetry 
restoration parameter //, but also on the topological charge v. The dependence on these two 
variables can explicitly be calculated. It turns out that the two-point correlation functions 
considered above retain the form (55), but the coefficients cia, bA, ap, bp now depend on 
the topological charge as well as on the symmetry restoration parameter (see [202-204] for 
explicit expressions). 

A specific issue with e-regime calculations is the scale setting. Ideally one would perform 
a p-regime study with the same bare parameters to measure a hadronic scale (e.g. the proton 
mass). In the literature, sometimes a gluonic scale (e.g. ro) is used instead; this introduces 
an extra uncertainty of the order of 5%. 

It is important to stress that in the e-expansion higher order finite-volume corrections 
might be significant, and the physical box size (in fm) should still be large in order to keep 
them under control. The criteria for the chiral extrapolation and finite volume effects are 
obviously different with respect to the p-regime. For these reasons we have to adjust the 
colour coding defined in Sect. 2 (see 5.1.6 for more details). 

Recently, the chiral effective theory has been extended also to the so-called mixed regime, 
where some quarks are in the p-regime and others in the e-regime [205] . In [206] a technique 
is proposed to smoothly connect p- and e-regimes. These theoretical advances can be adopted 
for future extraction of the Low Energy Couplings from lattice simulations. 
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5.1.3 Energy levels of the QCD Hamiltonian in a box, 5-regime 

At low temperature, the properties of the partition function are governed by the lowest 
eigenvalues of the Hamiltonian. In the case of QCD, the lowest levels are due to the Nambu- 
Goldstone bosons and can be worked out with xPT [207]. In the chiral limit, the level 
pattern follows the one of a quantum-mechanical rotator: Eg = £(£ + 2)/(2 0), £ = 0, 1,2, . . . 
For a cubic spatial box and to leading order in the expansion in inverse powers of the box 
size L s , the moment of inertia is fixed by the value of the pion decay constant in the chiral 
limit: = F 2 L^. In order to analyze the dependence of the levels on the quark masses 
and on the parameters that specify the size of the box, a reordering of the chiral series 
is required, the so-called <5-expansion; the region where the properties of the system are 
controlled by this expansion is referred to as the 5-regime. Evaluating the chiral perturbation 
series in this regime, one finds that the expansion of the partition function goes in even 
inverse powers of FL S , that the rotator formula for the energy levels holds up to NNLO and 
the expression for the moment of inertia is now also known up to and including order (FL S )~ 
[208-210]. Since the level spectrum is governed by the value of the pion decay constant in 
the chiral limit, an evaluation of this spectrum on the lattice can be used to measure F. 
More generally, the evaluation of various observables in the <5-regime offers an alternative 
method for a determination of some of the low energy constants occurring in the effective 
Lagrangian. At present, however, the numerical results obtained in this way [211, 212] are 
not yet competitive with those found in the p- or e-regimes. 

5.1.4 Other methods for the extraction of the Low-Energy Constants 

An observable that can be used to extract the LECs is the topological susceptibility, which 
is defined as 

X t= ld A x{uj{x)uj{U)), (56) 



where ui(x) is the topological charge density, 



^ ( '''''-Tr /•;,,,!,•!/•, ,i.r! . (57) 



At infinite volume, the expansion of the topological susceptibility in powers of the quark 
masses starts with [213] 

Xt = rnT,{l + 0(m)} , m = ( — + — + — + .. . (58) 

\m u m d m s J 

The quark condensate £ can thus be extracted from the properties of the topological sus- 
ceptibility close to the chiral limit. The behaviour at finite volume, in particular also in the 
region where the symmetry is restored, are discussed in [200]. The dependence on the vacuum 
angle 6 and the projection on sectors of given topological charge have also been investigated 
[201]. For a discussion of the finite-size effects to NLO, including the dependence on 9, we 
refer to [204, 214]. 

Another method to compute directly the quark condensate has been proposed in [215], 
where it is shown that starting from the Banks-Casher relation [216], it is possible to extract 
the condensate from suitable (renormalizable) spectral observables, for instance the number 
of Dirac operator modes contained in a given interval. For those spectral observables higher 
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order corrections can be systematically computed in terms of the chiral effective theory. 

An alternative strategy is based on the fact that at LO in the e-expansion, the partition 
function in a given topological sector v is equivalent to the one of a chiral Random Matrix 
Theory (RMT) [217-220]. In RMT it is possible to extract the probability distributions 
of individual eigenvalues [221-223] in terms of two dimensionless variables £ = A£V and 
H = rnSF, where A represents the eigenvalue of the massless Dirac operator and m is the sea 
quark mass. Hence it is possible to match the QCD low-lying spectrum of the Dirac operator 
with the RMT predictions in order to extract 15 the chiral condensate X. The main issue 
concerning this method is that for the distributions of individual eigenvalues higher order 
corrections are still not known by means of the chiral effective theory, and this may introduce 
systematic effects which are not under control. 16 Another open question is that, while it is 
clear how the spectral density is renormalized [227], this is not the case for the individual 
eigenvalues, and one has to rely on assumptions. There have been many lattice studies [228- 
232] where the matching of the low-lying Dirac spectrum with RMT has been investigated. 
In this review we don't include the results of the LECs obtained in this way. 17 

5.1.5 Pion form factors 

The scalar and vector form factors of the pion are defined by the matrix elements 

{^(p2)\qqW{pi)) = 6 ij F£(t), (59) 
(Tr^l^rVgK'CPi)) = ie^^+p^Fm, 

where the operators only contain the lightest two quark flavours, r , r 2 , r 3 are the Pauli 
matrices and t = (pi — P2) 2 denotes the momentum transfer. 

The vector form factor has been measured by several experiments for timelike as well 
as for spacelike values of t. The scalar form factor is not measurable but can be evaluated 
theoretically from data on the tttt and ttK phase shifts [233] on the basis of analyticity and 
unitarity, i.e. in a model-independent way. Lattice calculations can be compared with data 
or model-independent theoretical evaluations at different values of t. At present, however, 
most lattice calculations concentrate on the region close to t = and on the evaluation of the 
slope and curvature, defined as: 

F$(t) = l + \(r 2 y v t + c v t 2 + ... , (60) 
m) = F^0)[l + l(r 2 y s t + c s t 2 + ...] . 

The slopes are related to the mean quadratic vector and scalar radii which are the quantities 
on which most experiments and lattice calculations concentrate. 

In chiral perturbation theory, the form factors are known at NNLO [234]. The corre- 
sponding formulae are available in fully analytical form and are compact enough that they 
can be used for the chiral extrapolation of the data (as done, for example in [235, 236]). The 

15 By introducing an imaginary isospin chemical potential, this framework can be extended such that the low 
spectrum of the Dirac operator is sensitive also to the pseudoscalar decay constant F at LO [224]. 

16 Higher order systematic effects in the matching with RMT have been recently investigated in [225, 226]. 

17 The results obtained for E and F lie in the same range as the determinations reported in Tables 10 and 
11. 
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expressions for the scalar and vector radii and for the c<j v coefficients at two-loop level reads 
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and k<pg y 



and k csv are four constants independent of the quark masses. Their expression in 



terms of the t{ and of the 0(p ) constants Cj is known but will not be reproduced here. 



5.1.6 Results 

In this section we summarize the lattice results for the SU(2) couplings in a set of tables 
(10-13) and figures (7-9). The tables present our usual colour coding which summarizes the 
main aspects related to the treatment of the systematic errors of the various calculations. 

A delicate issue in the lattice determination of chiral LECs (in particular at NLO) which 
cannot be reflected by our colour coding is a reliable assessment of the theoretical error. We 
add a few remarks on this point: 

1. Using both the x and the £ expansion is a good way to test how the ambiguity of the 
chiral expansion (at a given order) affects the numerical values of the LECs that are 
determined from a particular set of data. For instance, to determine £4 (or A4) from 
lattice data for F n as a function of the quark mass, one may compare the fits based 
on the parametrization F w = F{1 + xln(A 2 /M 2 )} [see Eq. (46)] with those obtained 
from = F/{1 — £ ln(A 2 /M 2 )} [see Eq. (51)]. The difference between the two results 
provides an estimate of the uncertainty due to the truncation of the chiral series. 18 

18 Notice however that this difference could be accidentally small and then lead to an underestimate of the 
true systematic error. 
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Which central value one chooses is in principle arbitrary, but we find it advisable to use 
the one obtained with the £ expansion, 19 in particular because it makes the comparison 
with phenomenological determinations (where it is an established practice to use the £ 
expansion) more meaningful. The fact that there is an apparent discrepancy in Tab. 13 
between the lattice and the phenomenological determination of £\— £2, but not for (r 2 )g 
(on which it is based), calls for a re-assessment of the theoretical error which includes 
an x versus £ expansion analysis. 

2. As an alternative one could try to estimate the influence of higher chiral orders by 
reshuffling irrelevant higher-order terms. For instance, in the example mentioned above 
one might use F w = F/{1 — a; ln(A|/M 2 )} as a different functional form at NLO. Another 
way to make this estimate is through introducing by hand "analytic" higher-order terms 
(e.g. "analytic NNLO"), as done, in the past, by MILC [6]. In principle it would be 
preferable to include all NNLO terms or none, such that the structure of the chiral 
expansion is preserved at any order; this is what ETM [237] and JLQCD/TWQCD [95] 
are doing now for SU(2) xPT and MILC for SU(3) [59]. In any case it is not 
advisable to include terms to which the data are not sensitive. The use of priors in fits 
has not yet found general consensus in the lattice community. 

3. Another issue concerns the s-quark mass dependence of the LECs or Aj in the SU(2) 
framework. This has been studied analytically in a series of papers [91, 238, 239]. An 
analysis of the difference of the LECs determinations done in Nf = 2 or Nf = 2 + 1 
lattice simulations provides a unique means to analyze quantitatively this dependence. 

4. Finally, the determination of the LECs is in general affected by discretization effects, 
and it is important that these are corrected for with the help of the appropriate effective 
chiral Lagrangian (as done, e.g. by MILC). For actions which respect chiral symmetry 
at finite lattice spacing (like the overlap formulation adopted by JLQCD/TWQCD or 
domain wall fermions used by RBC/UKQCD in the limit of large lattice size in the fifth 
dimension), this step is unnecessary. 

In the tables and figures we summarize the results of various lattice collaborations for the 
SU(2) LECs at leading order (F or F/F n , B or S) and at NLO (I x - £ 2 , £3, h, k, h)- The 
tables group the results into those which stem from Nf = 2 + 1 calculations and those which 
stem from Nf = 2 calculations. Furthermore, we make a distinction whether the results are 
obtained from simulations in the p-regime or whether alternative methods (e-regime, spectral 
quantities, topological susceptibility, etc.) have been used. For comparison we add, in each 
case, a few phenomenological determinations with high standing. 

There is only one set of Nf = 2 lattice results for the SU(2) LECs which does not have any 
red tag, the one by the ETM collaboration [235, 237]. These results are at present the lattice 
determinations of these quantities with the best control of systematic effects (for Nf = 2 and 
in some cases even overall) and the only ones which qualify for making an average according 
to our criteria. We therefore offer no Nf = 2 averages, and the dashed bands in the plots 
coincide with the lattice points of the ETM collaboration. 

19 There are theoretical arguments suggesting that the £ expansion is preferable to the x expansion, based 
on the observation that the coefficients in front of the squared logs in (46) are somewhat larger than in (51). 
This can be traced to the fact that a part of every formula in the x expansion is concerned with locating the 
position of the pion pole (at the previous order) while in the £ expansion the knowledge of this position is built 
in exactly. Numerical evidence supporting this view is presented in [95]. 
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Table 10: Lattice results for the quark condensate S, in the MS scheme at scale \x = 2 GeV. 
We separate Nf = 2 + 1 from Nf = 2 results and, in the latter case, results obtained in 
the p-regime (middle) from those obtained in the e-regime (bottom). For MILC 09A and 
PACS-CS 08, converted values from SU(3) fits and values from direct SU(2) fits are included. 



The determination of the SU(2) LECs in Nf = 2 + 1 lattice calculations requires more 
discussion: MILC has published values for E, F and F n /F in their summary paper [6], but 
not for the £i, which were instead extracted from a SU(2) fit only in conference proceedings 
[59, 103]. As discussed in the previous sections, these results have no red tag and those 
presented in a regular article fully qualify for entering our averages. The RBC/UKQCD 
collaboration has calculated all the SU(2) LECs discussed here. After a first study with a 
single lattice spacing [108], a more recent publication with 2 lattice spacings appeared [106]. 
Those results have no red tag and are included in our averages. PACS-CS [63] has also 
performed a comprehensive analysis of SU(2) LECs, and at remarkably small quark masses - 
their calculation has two red tags related to the continuum extrapolation and finite-size effects 
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Figure 7: Summary of lattice results for the quark condensate £ (in the MS scheme at scale 
fx = 2 GeV). The filled symbols represent lattice simulations with Nf = 2 + 1 dynamical 
flavours, while the empty ones refer to Nt = 2 computations. Squares and left triangles 
indicate determinations from correlators in the p- and e-regimes, respectively. Up triangles 
refer to extractions from the topological susceptibility, diamonds to determinations from the 
pion form factor. The gray (dashed) band indicates our estimate for Nf = 2 + l(Nf = 2). 



and should also not be considered in averages. Finally JLQCD [240, 241] has performed a 
calculation of £ (but no other SU(2) LECs) with Nf = 2 + 1 overlap dynamical quarks, but 
at a single lattice spacings and has also a red tag. 

The ETM collaboration has recently started performing calculations with Nf = 2 + 1 + 1 
dynamical quarks and a first paper [60] on the SU(2) LECs has appeared (a more recent 
paper [249] gives a progress report but not yet an update of the LECs). Assigning a color 
code to this calculation is not straightforward because of the even larger splitting between 
the neutral and charged pion than observed in Nf = 2 calculations, but as explained in the 
footnote in Table 11 we did not assign it any red tag. We also refer to the careful analysis 
by Bar [40] on the effect of this splitting on the LEC determination. All this considered, we 
compare and take into account their numbers when drawing our conclusions on Nf = 2 + 1 
calculations of the SU(2) LECs. 

A generic comment applies to the issue of the scale setting. Since none of the lattice 
studies involves simulations in the p-regime at the physical value of m uc i, the setting of the 
scale a -1 via an experimentally measurable quantity involves a chiral extrapolation. As a 
result, dimensionful quantities are particularly sensitive to this scale-setting ambiguity, while 
in dimensionless ratios such as F n /F, F/Fq, B/Bq, £/£o> the problem is much reduced, and 
often finite lattice-to-continuum renormalization factors drop out. In those cases where the 
collaboration offers separate results for numerator and denominator but not for the ratio, 
we have calculated the quotient ourselves. Since the systematic errors are dominated by 
those due to the renormalization of the scale and these cancel in F n /F, we have dropped the 
systematic errors altogether - a lattice determination of the ratio that properly accounts for 
the correlations would of course be preferable. In the tables, results obtained in this way are 
shown in slanted fonts. Since the difference between F v and F manifests itself only at NNLO 
of the e-expansion, the data taken in the e-regime do not by themselves allow a determination 
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Figure 8: Summary of lattice results for the ratio F n /F, which compares the physical value of 
the pion decay constant with the value in the limit m u , md — > 0. The meaning of the symbols 
is the same as in Fig. 7. The gray band here represents our estimate, whereas the dashed one 
the current Nf = 2 value. 



of the ratio F n /F. The corresponding entries in Tab. 11 and Fig. 8 are obtained by dividing 
the physical value of F w (cf. footnote 13) through the lattice result for F. 

It is worth repeating here that the standard colour coding scheme of our tables is nec- 
essarily schematic and cannot do justice to every calculation. For instance, in the e-regime 
we routinely assign 20 a green star in the "chiral extrapolation" column (since the pertinent 
expressions directly involve F and B, even if the calculation is done at finite quark mass), 
while the "infinite volume" assessment is exclusively based on the number of different volumes 
with L > 1.5 fm have been used (the M^L criterion does not make sense here). Similarly, in 
the calculation of form factors and charge radii the tables do not reflect whether an interpola- 
tion to the desired q 2 has been performed or whether the evaluation has taken place directly 
at the relevant q 2 , or very close to it, by means of "partially twisted boundary conditions" 
[250]. Nevertheless, we feel that these tables give an adequate overview of the qualities of the 
various calculations. 

For most quantities our tables and figures show an overall reasonable consistency among 
the various lattice determinations. One notable exception is the case of S. The PACS-CS 08 
determination [63] of this quantity is in flat disagreement with all other calculations. A 
glimpse at Tab. 10 reveals that this calculation is the only one which received three red boxes. 
In particular, as was pointed out by the authors, for this quantity perturbative renormalization 
may lead to theoretical errors which are hard to quantify. But even among the two Nf = 2 + 1 
calculations without red tags, MILC 10A and RBC/UKQCD 10A there is some tension, at 
the level of 2.6 sigmas. As pointed out above, the origin of the problem may lie in the scale 
setting issue: the ratio £/£o would be a much more stringent test. In any case, since it is not 
meaningful to average these two calculations, the estimate shown in Figure 7, £ = 269(18) is 
constructed so as to cover both of them. 

As far as F n /F is concerned, there is some tension among the MILC and the RBC/UKQCD 
results, with the latest ETM Nf = 2 + 1 + 1 agreeing with the latter. The tension is at the 

20 Also for [241] and for [240] (for Nf = 2 + 1,3) the colour-coding criteria for the e-regime have been applied. 
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1.02(5) 


HHS 08 


[247] 


2 


A 


* 


■ 


• 


* 


90(4) 


1.02(5) 


JLQCD/TWQCD 07 


[248] 


2 


A 


* 


■ 


■ 


★ 


87.3(5.6) 


1.06(7) 


CD 03 


[251] 














86.2(5) 


1.0719(52) 



' The values of M^+L correspond to a green tag, while those of M n oL imply a red one 

play a role in finite-volume effects, we opt for orange. 

§ Result for roF converted into a value for F with ro = 0.49 fm. 



since both masses 



Table 11: Results for the leading order SU(2) low energy constant F and for the ratio F n /F. 
Numbers in slanted fonts have been calculated by us (see text for details). Horizontal lines 
establish the same grouping as in Tab. 10. The table shows that the precision reached on the 
lattice is now comparable to the accuracy of the prediction, which is indicated in the last row. 



level of about two sigmas, but it is worth stressing that the quoted lattice uncertainties are 
all below one percent. We believe it is fair and useful to quote an estimate which covers all 
three lattice determinations and has an uncertainty of about 1.5%, 



Fn 

F 



1.073(15) our estimate. 



(63) 



The ETM result [237] shows that the value obtained with Nf = 2 dynamical flavours is the 
same within errors. 21 The crucial point here is that the lattice results confirm the theoretical 
prediction [251]: the breaking of chiral symmetry generated by the small masses of the two 
lightest quarks is significantly larger than what dimensional reasoning would indicate, because 
it is amplified by a chiral logarithm with a large coefficient. The amplification is generated by 
the propagation of Nambu-Goldstone boson pairs that necessarily accompany the spontaneous 
breakdown of the symmetry. 

At NLO, the ratio F^/F is determined by the effective coupling £4. As we are dealing with 
an expansion in powers of the two lightest quark masses here, the corrections of higher order 



21 On the other hand, the result of JLQCD/TWQCD [95] shows a 2-sigma deviation from the other results 
obtained with Nf = 2. Since the quoted value for F^/F has been obtained by dividing the physical value of F„ 
by the lattice result, we can't exclude that systematic errors are properly taken into account in this estimate. 
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are expected to be very small (not necessarily at the quark masses used in the simulations, 
but at their physical values). Unfortunately, however, some of the collaborations appear to 
underestimate the uncertainties in the determination of this coupling: in view of the small 
errors quoted, the results are not coherent (see the middle panel in Fig. 9). The dust needs 
to settle before a meaningful estimate can be given for £4. 

The leading term in the chiral expansion of the scalar radius (r 2 )g is also determined by 
the coupling £4. The bottom panel of Table 13 shows that the result for the scalar radius 
obtained from a lattice calculation of the scalar form factor by JLQCD/TWQCD 09 [252] is 
in good agreement with the result of CGL 01 [198], which is based on dispersion theory, but 
the red tags indicate that some of the systematic uncertainties in that calculation yet need 
to be studied. 

The current situation with £3 is better than the one with £4: the N* = 2+1 determinations 
agree well with one another and the two ETM calculations for Nf = 2 and Nj = 2 + 1 + 1 
corroborate the result further. There is some discrepancy between two among the most recent 
determinations which have rather small error bars, namely ETM 10 [60] and RBC/UKQCD 
10A [106]. We stress, however, that the value: 

£ 3 = 3.2(8) our estimate , (64) 

comfortably covers all available calculations and represents a substantial progress with respect 
to earlier estimates. The uncertainty is three times smaller than the one attached to the old 
estimate in [93], which was based on analysis of the mass formulae for the pseudoscalar octet. 

Table 13 collects the results obtained for the electromagnetic form factor of the pion: 
charge radius, curvature and the coupling constant £q, that determines the charge radius to 
leading order of the chiral expansion. The experimental information concerning the charge 
radius is excellent and the curvature is also known very precisely, on the basis of e + e~ data 
and dispersion theory. The form factor calculations thus present an excellent testing ground 
for the lattice calculations. The table shows that most of the available lattice results pass 
the test. Concerning the value of £q, the situation is worse than for the charge radius also 
because the extraction of this constant is done differently by different groups (in particular 
some use the x rather than the ^ expansion and NLO rather than NNLO formulae): here the 
dust yet needs to settle. 

Perhaps the most important physics result of this section is that the lattice simulations 
confirm the approximate validity of the Gell-Mann-Oakes-Renner formula and show that the 
square of the pion mass indeed grows in proportion to m uc i- The formula represents the 
leading term of the chiral perturbation series and necessarily receives corrections from higher 
orders. At first non-leading order, the correction is determined by the effective coupling 
constant £3. The results collected in Tab. 12 and in the top panel of Fig. 9 show that £3 is now 
known quite well. They corroborate the conclusion drawn in the pioneering work of reference 
[253]: the lattice confirms the early estimate of £3 derived in [93]. In the graph of M% versus 
m u d, the values found on the lattice for £3 correspond to remarkably little curvature: the 
Gell-Mann-Oakes-Renner formula represents a crude first approximation out to values of m u d 
that exceed the physical value by an order of magnitude. 

As emphasized by Stern and collaborators [255-257], the analysis in the framework of 
xPT is coherent only if (i) the leading term in the chiral expansion of M% dominates over 
the remainder and (ii) the ratio m s /m u d is close to the value 25.6 that follows from Wein- 
berg's leading order formulae. In order to investigate the possibility that one or both of 
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Collaboration 



Ref. 



o 



ETM 10 



[60] 



2+1+1 A 



3.70(7) (26) 



4.67(3)(10) 



MILC 10A 


[103] 


2+1 


C 




★ 


★ 


2-85(81) (±«) 


3.98(32) 


MILC 10 


[171] 


2+1 


C 


• 


★ 


★ 


3.18(50)(89) 


4.29(21)(82) 


RBC/UKQCD 10A 


[106] 


2+1 


P 


• 


• 


★ 


2.57(18) 


3.83(9) 


MILC 09A 


[59] 


2+1 


c 


• 


★ 


★ 


3.32(64)(45) 


4.03(16) (17) 


MILC 09A 


[59] 


2+1 


c 


• 


★ 


* 


3-0(6) (±J) 


3.9(2)(3) 


PACS-CS 08 


[63] 


2+1 


A 


★ 


■ 


■ 


3.47(11) 


4.21(11) 


PACS-CS 08 


[63] 


2+1 


A 


★ 


■ 


■ 


3.14(23) 


4.04(19) 


RBC/UKQCD 08 


[108] 


2+1 


A 


• 


■ 


★ 


3.13(33) (24) 


4.43(14) (77) 


ETM 09C 


[237] 


2 


A 


• 


★ 


• 


3-50(9) (_+ 3 9 ) 


4.66(4) (_ + 3 4 3 ) 


JLQCD/TWQCD 09 


[252] 


2 


A 


• 


■ 


■ 




4.09(50)(52) 


ETM 08 


[235] 


2 


A 


• 


• 


• 


3.2(8)(2) 


4.4(2)(1) 


JLQCD/TWQCD 08 A 


[95] 


2 


A 


• 


■ 


■ 


3.38(40) (24) 


4.12(35) (30) (+_«) 


CERN-TOV 06 


[254] 


2 


A 


• 


• 


■ 


3-0(5)(l) 





CGL 01 


[198] 




4.4(2) 


GL 84 


[93] 


2.9(2.4) 


4.3(9) 



Table 12: Results for the coupling constants £3 and £4 of the effective SU(2) Lagrangian. For 
MILC 09A and PACS-CS 08, converted values from SU(3) fits and values from direct SU(2) 
fits are included. The MILC 10 results are obtained by converting the SU(3) LECs, while the 
MILC 10A results are obtained with a direct SU(2) fit. For comparison, the last two lines 
show the results obtained from a phenomenological analysis. 
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of* 

$ i 



o 



Collaboration 



Ref. N, 



f 



(r 2 Vv [fm 2 



c v (GeV" 



RBC/UKQCD 08A [250] 


2+1 


A 




■ 


★ 


0.418(31) 




12.2(9) 


LHP 04 [259] 


2+1 


A 




■ 




0.310(46) 






JLQCD/TWQCD 09 [252] 


2 


A 


• 


■ 


■ 


0.409(23) (37) 


3.22(17)(36) 


11.9(0.7) (1.0) 


ETM 08 [235] 


2 


A 


• 


• 


• 


0.456(30) (24) 


3.37(31) (27) 


14.9(1.2) (0.7) 


QCDSF/UKQCD 06A[260] 


2 


A 


★ 


• 


• 


0.441(19) (56) (29) 






BCT 98 [234] 












0.437(16) 


3.85(60) 


16.0(0.5) (0.7) 


NA7 86 [261] 












0.439(8) 






GL 84 [93] 
















16.5(1.1) 



Collaboration 



Ref. 



{r 2 )l 



fm 2 



JLQCD/TWQCD 09 


[252] 2 


A 


• ■ 


■ 0.617(79)(66) -2.9(0.9) (1.3) 


CGL 01 


[198] 






0.61(4) -4.7(6) 


Table 13: Top panel: 


vector form 


factor 


of the pion. 


Lattice results for the charge radius 



(r 2 )y, the curvature cy and the effective coupling constant 1% are compared with the ex- 
perimental value obtained by NA7 and some phenomenological estimates (£q determines the 
leading contribution in the chiral expansion of the charge radius) . Bottom panel: scalar form 
factor of the pion. Lattice results for the scalar radius (r 2 ) J and the combination i\ — £2 of 
effective coupling constants are compared with the outcome of a dispersive calculation of these 
quantities [198] (the leading term in the chiral expansion of the scalar radius is determined 
by the coupling constant £4, for which the available information is collected in Tab. 12; the 
lattice estimate for £\ — £2 stems from an analysis of the momentum dependence of the vector 
and scalar form factors, based on the two-loop formulae of ;\TT [234]). 



these conditions might fail, the authors proposed a more general framework, referred to as 
"Generalized xPT", which includes xPT as a special case. The results found on the lattice 
demonstrate that QCD does satisfy both of the above conditions - in the context of QCD, 
the proposed generalization of the effective theory does not appear to be needed. There is a 
modified version, however, referred to as "Resummed xPT" [258], which is motivated by the 
possibility that the Zweig rule violating couplings L4 and Lq might be larger than expected. 
The available lattice data do not exclude this possibility (see section 5.2.4). 
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3 4 5 



ETM 10 

MILC 10A [SU(2) fit] 
MILC10 [SU(3) fit] 
RBC/UKQCD 10A 
MILC 09A [SU(2) fit] 
MILC 09A [SU(3) fit] 
PACS-CS 08 [SU{2) fit] 
PACS-CS 08 [SU<3) fit] 
RBC/UKQCD 08 



ETM 09C 
ETM 08 

JLQCD/TWQCD 08A 
CERN-TOV 06 

our estimate 



GL84 



3 4 



2.5 



3 3.5 4 4.5 5 




ETM 10 



MILC 10A [SU(2) fit] 
MILC 10 [SU(3) fit] 
RBC/UKQCD 10A 
MILC09A [SU(3) fit] 
MILC 09A [SU(2) fit] 
PACS-CS 08 [SU(3) fit] 
PACS-CS 08 [SU(2) fit] 
RBC/UKQCD 08 



H <^ h- 

h □— — i 



ETM 09C 

JLQCD/TWQCD 09 
ETM 08 

JLQCD/TWQCD 08A 



J i I i L 



CGL01 
GL84 



2.5 



3 3.5 4 4.5 5 
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~\ 1 r 
-♦ — i 



RBC/UKQCD 08A 



JLQCD/TWQCD 09 
ETM 08 




Figure 9: Summary of lattice results for the effective coupling constants ^3,^4 and £q. The 
filled symbols refer to lattice calculations with Nf = 2 + 1, while the empty ones represent 
simulations with Nf = 2. Squares indicate determinations from correlators in the p-regime, 
diamonds refer to determinations from the pion form factor. The gray band represents our 
estimate of £3, whereas the dashed ones indicate the current Nf = 2 values of £3 and £4. 
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5.2 SU(3) Low-Energy Constants 

5.2.1 Quark-mass dependence of pseudoscalar masses and decay constants 

In the isospin limit, the relevant SU(3) formulae take the form [91] 

Ml = 2B m ud {l + fin - -fir, + 16m ud (2L 8 -L 5 ) + 16(m s +2m ud )(2L e -L 4 ) |, 
M K = B Q (m s +m ud )\l+'^fi ri + ^ 8(m s +m ud )(2L 8 - L 5 ) + l6(m s +2m ud )(2L 6 - L 4 ) ) , 

nlo r -^o r 1 1 

F-k = Fq^I - - hk + -jjq [ 8m udL5 + 8(m s +2m ud )L^ j , (65) 
Fr = F Q ^l- ^/x^ - ^fi K - + 4:(m s +m ud )L 5 + 8(m s + 2m ud )L 4 | , 



where Bq = Y^q/Fq and Fq denote the condensate parameter and the pseudoscalar decay 
constant in the SU(3) chiral limit, respectively, and 

Ml , [Ml\ , . 

"^^fH/)- < 66 > 

At the accuracy of these formulae, the quantities fix, fJt] can equally well be evaluated 
with the leading order expressions for the masses, 

, LO LO x r, LO r, s . 

Ml = 2B m ud , M K = B (m s +m ud ) , M< = | J B (2m s +m ud ) . (67) 

Throughout, Lj denotes the renormalized low-energy constant/coupling (LEC) at scale \x. 
The normalization used for the decay constants is specified in footnote 13. 

5.2.2 Charge radius 

The SU(3) formula for the slope of the pion vector form factor reads [167] 

, 9v LO If , Ml, , 12L 9 , . 

{r = -3Sw{ 3 + 2 ln( 7#» + ln< 7# ) } + V • (68) 

while the expression (r 2 )s iOC t for the octet part of the scalar radius does not contain any NLO 
low-energy constant at the 1-loop order [167]. 

5.2.3 Partially quenched formulae 

The term "partially quenched QCD" is used in two ways. For heavy quarks (c, b and sometimes 
s) it usually means that these flavors are included in the valence sector, but not to the 
functional determinant. For the light quarks (u,d and sometimes s) it means that they are 
present in both the valence and the sea sector of the theory, but with different masses (e.g. a 
series of valence quark masses is evaluated on an ensemble with a fixed sea quark mass). 

The program of extending the standard (unitary) SU(3) theory to the (second version 
of) "partially quenched QCD" has been completed at the 2-loop (NNLO) level for masses 
and decay constants [262]. These formulae tend to be complicated, with the consequence 
that a state-of-the-art analysis with 0(2000) bootstrap samples on O(20) ensembles with 
0(5) masses each [and hence O(200'000) different fits] will require significant computational 
resources for the global fits. For an up-to-date summary of recent developments in Chiral 
Perturbation Theory relevant to lattice QCD we refer to [263]. 
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o 







Ref. N, 



f 



7 $ 



J 



F [MeV] 



F/F 



B/B 



JLQCD/TWQCD 10[240] 3 A ■ ■ ■ ★ 71 (3) (8) 



MILC 10 [171] 2+1 C • * * • 80.3(2. 5)(5.4) 

MILC 09A 22 [59] 2+1 C • ★ ★ • 78.3(1.4) (2.9) 1.104(3)(41) 1.21(4)(+ 5 6 ) 

MILC 09 22 [6] 2+1 A • ★ ★ • 1.15(5)(±* 3 ) 1.15(16)(tf 3 ) 

PACS-CS 08 [63] 2+1 A ★ ■ ■ ■ 83.8(6.4) 1.078(44) 1.089(15) 

RBC/UKQCD 08 [108] 2+1 A ■ ★ ★ 66.1(5.2) 1.229(59) 1.03(05) 



Ref. 



N f 



,1/3 



[GeV] 



S/E 



JLQCD/TWQCD 10 


[240] 


3 


A 


■ 


■ 


■ 


★ 


0.214(6) (24) 


1.31(13)(52) 


MILC 09A 22 
MILC 09 22 
PACS-CS 08 
RBC/UKQCD 08 


[59] 

[6] 
[63] 
[108] 


2+1 
2+1 
2+1 
2+1 


C 
A 
A 
A 


• 
• 

★ 
• 


★ 
★ 
■ 
■ 


★ 
★ 
■ 

★ 


• 
• 
■ 

* 


0.245(5) (4) (4) 

0.242(9)(+ 5 7 )(4) 

0.290(15) 


1.48(9) (8) (10) 
1.52(17) (±f 6 ) 
1.245(10) 
1.55(21) 



Table 14: Summary of lattice results for the low-energy constants Fq, Bq and Eo = F^Bq, 
which specify the effective SU(3) Lagrangian at leading order. The ratios F/Fq, B/Bq, E/Sq, 
which compare these with their SU(2) counterparts, indicate the strength of the Zweig-rule 
violations in these quantities: in the large--/V c limit, they tend to unity. Numbers in slanted 
fonts are calculated by us, from the information given in the quoted references. 



5.2.4 Recent lattice determinations 

To date, there are three comprehensive papers with results based on lattice QCD with Nf = 2+1 
dynamical flavors [6, 63, 108]. It is an open issue whether the data collected at m s ~ms hys 
allow for an unambiguous determination of SU(3) low-energy constants (cf. the discussion in 
[108]). So far only MILC has some data at considerably smaller m s [6]. Furthermore, we are 
aware of a few papers with a result on one SU(3) low-energy constant each [176, 250, 264]. 
Some particulars of the comprehensive papers are shown in Tab. 14. 

Results for the SU(3) low-energy constants of leading order are found in Tab. 14 and 
analogous results for some of the effective coupling constants that enter the chiral SU(3) 
Lagrangian at NLO are collected in Tab. 15. From PACS-CS [63] only those results are 
quoted which have been corrected for finite-size effects (misleadingly labeled "w/FSE" in 
their tables). For staggered data our "chiral" colour coding criterion is slightly ambiguous; 
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Ref. N-, 



o 



4 



10 3 L 4 



10 3 L 6 



10 3 (2L 6 -L 4 ) 



JLQCD/TWQCD 10 [240] 



A 



0.03(7)(17) 



MILC 10 


[171] 


2+1 


C 


• 


★ 


★ 


-0.08(22)(li) 


-0.02(16)(±ll) 


0.03(24) (lf 7 ) 


MILC 09A 22 


[59] 


2+1 


c 


• 


★ 


★ 


0.04(13) (4) 


0.07(10)(3) 


0.10(12)(2) 
0.3(1)(1 2 ) 


MILC 09 22 


[6] 


2+1 


A 


• 


★ 


★ 


0.1(3)(t?) 


0.2(2)(± 2 ) 


PACS-CS 08 


[63] 


2+1 


A 


★ 


■ 


■ 


-0.06(10)(-) 


0.02(5)(-) 


0.10(2)(-) 


RBC/UKQCD 08 


[108] 


2+1 


A 


• 


■ 


★ 


0.14(8)(-) 


0.07(6)(-) 


0.00(4)(-) 



Bijnens 09 


[263] 


0.0(-)(5) 


0.0(-)(l) 


0.0(-)(5) 


GL85 


[91] 


-0.3(5) 


-0.2(3) 


-0.1(8) 



Ref. N, 



f 



o 



ft? 



10 3 L 8 



10 3 (2L 8 



MILC 10 


[171] 


2+1 


C 


• 


★ 


★ 


0.98(16)(t2) 


0.42(10)(tf 3 ) 


-0.15(11)(±«) 


MILC 09A 22 


[59] 


2+1 


C 


• 


★ 


★ 


0.84(12)(36) 


0.36(5) (7) 


-0.12(8)(21) 


MILC 09 22 


[6] 


2+1 


A 


• 


★ 


★ 


1.4(2)(± 2 ) 


0.8(1)(1) 


0.3(1)(1) 


PACS-CS 08 


[63] 


2+1 


A 


★ 


■ 


■ 


1.45(7)(-) 


0.62(4)(-) 


-0.21(3)(-) 


RBC/UKQCD 


3 [108] 


2+1 


A 


• 


■ 


★ 


0.87(10)(-) 


0.56(4)(-) 


0.24(4) (-) 



Bijnens 09 


[263] 


0.97(11)(-) 


0.60(18)(-) 


0.23(38) (-) 


GL85 


[91] 


1.4(5) 


0.9(3) 


0.4(8) 



Ref. 



N f 



J 



o 



io-% 



10 3 L! 



JLQCD 08A 
RBC/UKQCD 08A 
NPLQCD 06 


[264] 
[250] 
[176] 


2 

2+1 
2+1 


A 
A 
A 


■ ■ 

■ ★ 

■ ■ 


1-42(2)0 


3.08(23)(51) 


-5.2(2)(t«) 


Bijnens 09 
GL85 


[263] 
[91] 








0.97(11)(-) 
1.4(5) 


5.93(43) (-) 
6.9(7) 


-5.5(7) 



Table 15: Summary of lattice results for some of the coupling constants that enter the effective 
SU(3) Lagrangian at NLO (running scale /i = 770 MeV - the values in [6, 59, 91, 171] are 
evolved accordingly). The PACS-CS entry for Lg is obtained from their results for 1L§ — L/± 
and L4 (and similarly for other entries in slanted fonts). 
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Figure 10: Summary of various lattice determinations of SU(3) low-energy constants. 

we solve this issue pragmatically. 22 

A graphical summary of the lattice results for the coupling constants L4, L5, Lq and 



The "chiral extrapolation" rating of [6] is based on the information on the RMS masses given in [59]. 
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L$, which determine the masses and the decay constants of the pions and kaons at NLO 
of the chiral SU(3) expansion, is displayed in Fig. 10, along with the two phenomenological 
determinations quoted in the above tables. In general, there is a rather convincing overall 
consistency, although some of the lattice results possibly underestimate the theoretical error. 
Since only the calculation by the MILC collaboration is free of red tags, the gray bands in 
Fig. 10 mark their results. 

In spite of this apparent consistency, there is a point which needs to be clarified as soon 
as possible. Some collaborations (RBC/UKQCD and PACS-CS) find that they are having 
difficulties in fitting their partially quenched data to the respective formulas for pion masses 
above ~ 400 MeV. Evidently, this indicates that the data are stretching the regime of validity 
of these formulas. To date it is, however, not clear which subset of the data causes the 
troubles, whether it is the unitary part extending to too large values of the quark masses or 
whether it is due to m val /m sea differing too much from one. In fact, little is known, in the 
framework of partially quenched xPT, about the shape of the region of applicability in the 

m val versus m sea pl ane f or fixed Nf. 

In the large-iV c limit, the Zweig-rule becomes exact, but the quarks have iV c = 3. The 
work done on the lattice is ideally suited to disprove or confirm the approximate validity 
of this rule. Two of the coupling constants entering the effective SU(3) Lagrangian at NLO 
disappear when N c is sent to infinity: L4 and Lq. The upper part of Tab. 15 and the left panels 
of Fig. 10 show that the lattice results for these are quite coherent. At the scale /i = M p , 
L4 and Lq are consistent with zero, indicating that these constants do approximately obey 
the Zweig-rule. As mentioned above, the ratios F/Fq, B/Bq and E/So also test the validity 
of this rule. Their expansion in powers of m s starts with unity and the contributions of first 
order in m s are determined by the constants L4 and Lq, but they also contain terms of higher 
order. Quite apart from measuring the Zweig-rule violations, an accurate determination will 
thus also allow us to determine the range of m s where the first few terms of the expansion 
represent an adequate approximation. Unfortunately, at present, the uncertainties in the 
lattice data on the ratios are too large to draw conclusions, both concerning the relative size 
of the subsequent terms in the chiral perturbation series and concerning the magnitude of 
the Zweig-rule violations. The data do appear to confirm the paramagnetic inequalities [257], 
which require F/Fq > 1, X/So > 1 and it appears that the ratio B/Bq is also larger than 
unity, but the numerical results need to be improved before further conclusions can be drawn. 

In principle, the matching formulae in [91] can be used to calculate the SU(2) couplings 
li from the SU(3) couplings Lj. 23 This procedure, however, yields less accurate results than 
a direct determination within SU(2), as it relies on the expansion in powers of m s , where 
the omitted higher order contributions generate comparatively large uncertainties. We plead 
with every collaboration performing Nt = 2 + 1 simulations to directly analyze their data in 
the SU(2) framework. In practice, lattice simulations are performed at values of m s close to 
the physical value and the results are then corrected for the difference of m s from its physical 
value. If simulations with more than one value of m s have been performed, then this can be 
done by interpolation. Alternatively one can use the technique of reweighting (for a review 
see [5]), which applies an a posteriori correction to shift m s into its physical value. 

23 For instance, for the MILC data this yields T 3 = 3.32(64) (45) and l 4 = 4.03(16)(17) [59]. 
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6 Kaon 5-parameter Bk 



6.1 Indirect CP-violation and e K 

The mixing of neutral pseudoscalar mesons plays an important role in the understanding of 
the physics of CP-violation. In this section we will only focus on K° — K° oscillations, which 
probe the physics of indirect CP-violation. We collect and comment the basic formulae; for 
extended reviews on the subject see, among others, refs. [265-267]. Indirect CP-violation 
arises in Kl — > irir transitions through the decay of the CP = +1 component of Kl into two 
pions (which are also in a CP = +1 state). Its measure is defined as 

A[K L -» (^r) J=0 ] 
eK A[K S -+ (irir) I=0 ] ' [ ™ } 

with the final state having total isospin zero. The parameter ex may also be expressed in 
terms of K° — K° oscillations. In particular, to lowest order in the electroweak theory, the 
contribution to these oscillations arises from so-called box diagrams, in which two W^-bosons 
and two "up-type" quarks (i.e. up, charm, top) are exchanged between the constituent down 
and strange quarks of the X-mesons. The loop integration of the box diagrams can be 
performed exactly. In the limit of vanishing external momenta and external quark masses, 
the result can be identified with an effective four-fermion interaction, expressed in terms of 
the "effective Hamiltonian" 

H AS=2 = 9^Mk^ Q AS=2 + h c (?0) 

In this expression, Gp is the Fermi coupling, Jl% the VF-boson mass, and 

Q^=2 = [s 7m (1 - 75 )rf] [s 7 ^(l - J5 )d\ = OyV+AA " OVA+AV , (71) 

is a dimension-six, four-fermion operator. The function J 70 is given by 

T° = X 2 c S (x c ) + X 2 S (x t ) + 2X c X t S (x c ,x t ) , (72) 

where A a = V* s V a d, and a = c,t denotes a flavour index. The quantities So(x c ), So(x t ) and 
So(x c ,xt) with x c = rric/Myy, xt = m^/M^ are the Inami-Lim functions [268], which express 
the basic electroweak loop contributions without QCD corrections. The contribution of the 
up quark, which is taken to be massless in this approach, has been taken into account by 
imposing the unitarity constraint X u + A c + Xt = 0. 

When strong interactions are included, AS = 2 transitions can no longer be discussed 
at the quark level. Instead, the effective Hamiltonian must be considered between mesonic 
initial and final states. Since the strong coupling constant is large at typical hadronic scales, 
the resulting weak matrix element cannot be calculated in perturbation theory. The operator 
product expansion (OPE) does, however, factorize long- and short- distance effects. For 
energy scales below the charm threshold, the K° — K° transition amplitude of the effective 
Hamiltonian can be expressed as 

fi2 M 2 r 

(^°|^ e A /= 2 |K°) = X 2 c S (x c ) m + X 2 t S (x t ) m + 2X c XtS (x c ,x t ) m 
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Air J { (Air 



2/?o 2 



(K°\Q& S=2 (ri\K») + h.c. , (73) 
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where g{jJi) and Q A5=2 (//) are the renormalized gauge coupling and four-fermion operator in 
some renormalization scheme. The factors 771,772 and 773 depend on the renormalized coupling 
g, evaluated at the various flavour thresholds nit, nib, 77i c and M\y, as required by the OPE and 
RG-running procedure that separates high- and low-energy contributions. Explicit expressions 
can be found in [266] and references therein. We follow the same conventions for the RG- 
equations as in ref. [266]. Thus the Callan-Symanzik function and the anomalous dimension 
7(5) of Q^ s=2 are defined by 



dg 
dln/i 

with perturbative expansions 
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We stress that Po,Pi and 70 are universal, i.e. scheme-independent. K° — K° mixing is 
usually considered in the naive dimensional regularization (NDR) scheme of MS, and below 
we specify the perturbative coefficient 71 in that scheme: 
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Note that for QCD the above expressions must be evaluated for N = 3 colours, while Nt 
denotes the number of active quark flavours. As already stated, eq. (73) is valid at scales 
below the charm threshold, after all heavier flavours have been integrated out, i.e. Nf = 3. 

In eq. (73), the terms proportional to 771, 772 and 773, multiplied by the contributions 
containing g(fj-) 2 , correspond to the Wilson coefficient of the OPE, estimated at NLO in 
perturbation theory. Its dependence on the renormalization scheme and scale ji is cancelled 
by that of the weak matrix element (K°\Q^ s=2 (fi)\K°) . The latter corresponds to the long- 
distance effects of the effective Hamiltonian and must be computed non-perturbatively. For 
historical, as well as technical reasons, it is convenient to express it in terms of the S-parameter 
Bk, defined as 

B M = ' 8 R f2 2 • (77) 

3 J K m K 

The four-quark operator Q A5=2 (/j) is renormalized at scale fj, in some regularization scheme, 
usually taken to be NDR. The renormalization group independent (RGI) S-parameter is 
related to Bk(h) by the exact formula 
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\ 4vr 
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At NLO in perturbation theory, the above reduces to 
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which, to this order, is the scale-independent product of all /j-dependent quantities in eq. (73). 

The implementation of non-perturbative renormalization in lattice calculations is based on 
the numerical evaluation of a non-perturbative matching condition between the bare operator 
matrix element and its counterpart in some intermediate renormalization scheme, where the 
latter has a controlled perturbative relation to continuum schemes such as MS. Examples of 
intermediate schemes are the RI/MOM scheme [115] (also dubbed the "Rome-Southampton 
method") and the Schrodinger functional (SF) scheme [114]. Typical scales for the transition 
between lattice regularization and intermediate scheme are /i = 0(1 GeV). When applied 
to the case at hand, this implies that -Bk(/J ~ 1 GeV) in the intermediate scheme must 
then be matched to the RGI l?-parameter via eq. (79). Obviously, due to asymptotic free- 
dom, the NLO relation provides a more reliable link if it is applied at scales far greater than 
jj, = 0(1 GeV). Therefore, in order to remove any doubts about the use of perturbation theory 
at NLO when determining B^, it is advantageous to run non-perturbatively to larger values 
of /z, where eq. (79) provides an accurate matching relation. Indeed, neglecting higher orders 
in perturbation theory has been identified by the authors of [269] as one of the dominant sys- 
tematic uncertainties in their result for B^. We note that the short-distance QCD corrections 
in eq. (73) are known in perturbation theory, to NLO for 771 , 772 and to NNLO for 773 [270]. 
This implies that an uncertainty of 0(a s (m c )) pertains to the K° — K° transition amplitude 
and ex, even if the computation of the kaon l?-parameter is essentially free from perturbative 
effects. However, since 771 and 772 are currently being determined beyond NLO as well [271], 
the issue of eliminating systematic uncertainties arising from the use of perturbation theory 
in the matching procedure becomes more acute. 

The "master formula" , for ex , which connects the experimentally observable quantity ex 
to the matrix element of n^ s s=2 , is [267, 272, 273] 

e K = exp(^) sm(0 £ ) [ — + j , (80) 

for A n real and positive; the phase of ex is given by 

AM K . . 

(f) € = arctan — — — . (81) 
ZYL K /Z 

The quantities AMk and AI^ are the mass- and decay width-differences between long- and 
short-lived neutral Kaons, while Aq is the amplitude of the Kaon decay into an isospin-0 two 
pion state. Experimentally known values of the above quantities are: 

\e K \ = 2.280(13) x 10~ 3 , 

<p e = 43.51(5)° , (82) 

AM K = 3.491(9) x 10~ 12 MeV , 

AT K = 7.335(4) x 10~ 15 GeV , 

while the last term in the square brackets of eq. (80) has been estimated in ref. [274]. Long 
distance contributions, omitted from the first term in the square brackets of eq. (80) have 
been recently discussed in ref. [275]. 

6.2 Lattice computation of Bk 

The lattice calculation of Bk is affected by the same systematic effects, discussed in pre- 
vious sections. However, the issue of renormalization merits special attention. The reason 
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is that the multiplicative renormalizability of the relevant operator Q AS=2 is lost, once the 
regularized QCD action ceases to be invariant under chiral transformations. With Wilson 
fermions, Q^ s=2 mixes with four additional dimension-six operators which belong to differ- 
ent representations of the chiral group, with mixing coefficients that are finite functions of the 
gauge coupling. This complicated renormalization pattern is the main source of systematic 
error in computations of Bk with Wilson quarks. It can be bypassed via the implementation 
of specifically designed methods, which are either based on Ward identities [276] or on a 
modification of the Wilson quark action, known as twisted mass QCD [52, 277]. An advan- 
tage of staggered fermions is the presence of a remnant U(l) chiral symmetry. However, at 
non-vanishing lattice spacing, the symmetry among the extra unphysical degrees of freedom 
(tastes) is broken. As a result, mixing with other dimension-six operators cannot be avoided 
in the staggered formulation, which complicates the determination of the l?-parameter. 

Fermionic lattice actions based on the Ginsparg- Wilson relation [278] are invariant under 
the chiral group, and hence four-quark operators such as Q^ s=2 renormalize multiplicatively. 
However, depending on the particular formulation of Ginsparg- Wilson fermions, residual chi- 
ral symmetry breaking effects may be present in actual calculations. For instance, in the case 
of domain wall fermions, the finiteness of the extra 5th dimension implies that the decou- 
pling of modes with different chirality is not exact, which produces a residual non-zero quark 
mass in the chiral limit. Whether or not a significant mixing with dimension-six operators is 
induced as well must be investigated on a case-by-case basis. 

In this section we focus on recent results for , obtained for Nf = 2 and 2 + 1 flavours of 
dynamical quarks. A compilation of results is shown in Table 16 and Fig. 11. An overview of 
the quality of systematic error studies is represented by the colour coded entries in Table 16. 
In Appendix B.4 we gather the simulation details and results from different collaborations, the 
values of the most relevant lattice parameters, and comparative tables on the various estimates 
of systematic errors. Note that some references do not quote results for both £>k and Bk- In 
this case we have performed the conversion ourselves by evaluating the proportionality factor 
in eq. (79) for [i = 2GeV. This requires fixing the value of a s (2GeV) in the two- and three- 
flavour cases. For Nf = 2 we have inserted the non-perturbative result for the A-parameter 
by the ALPHA Collaboration [279], i.e. A( 2 ) = 245(16) (16) MeV, into the NLO expression for 
a s , neglecting the quoted error. This gives Bk/Bk = 1.412, in accordance with ref. [280]. 
For Nf = 3 we use the value a s {Mz) = 0.1184 [122] and run it across the quark thresholds 
down to \i = 2GeV, using the four-loop RG /3-function. We then obtain Bk/Bk = 1.368 in 
the three-flavour theory. 

At present four collaborations have reported estimates of Bk for Nf = 2 + 1. The 
RBC/UKQCD collaboration uses domain wall fermions, while the results by HPQCD/UKQCD 
are based on the configurations generated by MILC, using Asqtad improved, rooted staggered 
fermions. In Aubin 09 a mixed-action approach is adopted in which staggered sea quarks 24 
are combined with HYP-smeared [289] domain wall fermions. The SWME result reported 
in [281] is also based on the Asqtad MILC ensembles, but uses HYP-smeared staggered quarks 
in the valence sector. 

In view of our quality criteria we note that one main shortcoming of the results by 
HPQCD/UKQCD 06 [285] is that they are based on data obtained at a single lattice spac- 
ing. This is also true for the results of RBC/UKQCD 07A, 08 [108, 284], which, however, 
have recently been updated by including data at a second value of a (see RBC/UKQCD 10B 

24 Aubin 09 also use configurations generated by MILC. 
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# 



.$ / S £ $ J° 

i f g / / # 

Collaboration Ref. N f <f o° c? <f B K B f 



SWME 11 


[281] 


2+1 


P 


★ 


• 


• 


■ 


- 0.523(7)(26) 


0.716(10)(35) 


RBC/UKQCD 10B 


[282] 


2+1 


P 


• 


• 


★ 


★ 


a 0.549(5)(26) 


0.749(7) (26) 


SWME 10 


[283] 


2+1 


A 


★ 


• 


• 


■ 


- 0.529(9)(32) 


0.724(12) (43) 


Aubin 09 


[269] 


2+1 


A 


• 




• 


★ 


- 0.527(6)(21) 


0.724(8)(29) 


RBC/UKQCD 07A, 08 


[108, 284] 


2+1 


A 


■ 


• 


★ 


★ 


- 0.524(10)(28) 


0.720(13)(37) 


HPQCD/UKQCD 06 


[285] 


2+1 


A 


■ 


* 


★ 


■ 


- 0.618(18) (135) 


0.83(18) 



ETM 10A [286] 2 A ★ • • ★ b 0.516(18)(12) 0.729(25)(17) 

JLQCD 08 [280] 2 A ■ • ■ * - 0.537(4)(40) 0.758(6)(71) 

RBC 04 [287] 2 A ■ ■ ■+ * - 0.495(18) 0.699(25) 

UKQCD 04 [288] 2 A ■ ■ B f ■ - 0.49(13) 0.69(18) 



* This result has been obtained with only two "light" sea quark masses. 

' These results have been obtained at (M^L) m i n > 4 in a lattice box with a spatial extension L < 2 fin. 
D In this mixed action computation, the lightest valence pion weighs ~ 230 MeV, while the lightest sea 

taste-pseudoscalar, used in the chiral fits, weighs ~ 370 MeV. 
a Bk is renormalized non-perturbatively at a scale of 2 GeV in a couple of Nf = 3 RI/SMOM schemes. A 

careful study of perturbative matching uncertainties has been performed by comparing results in the two 

schemes in the region of 2 GeV to 3 GeV [282]. 
b Bk is renormalized non-perturbatively at scales 1/a ~ 2 + 3 GeV in the Nf =2 RI/MOM scheme. In this 

scheme, non-perturbative and NLO running for the are shown to agree from 4 GeV down 2 GeV to better 

than 3% [102, 286]. 

Table 16: Results for the kaon S-parameter together with a summary of systematic errors. 
If information about non-perturbative running is available, this is indicated in the column 
"running", with details given at the bottom of the table. 



[282]). The good agreement between the earlier results and the latest update indicates that 
discretization effects for Bk, computed using domain wall fermions appear to be small at the 
present level of accuracy. Aubin 09 have also data at two values of the lattice spacing, simi- 
lar to those of RBC/UKQCD 10B, while the SWME 10 [283] and SWME 11 [281] results are 
obtained at three and four lattice spacings, respectively. Having two smaller lattice spacings 
at their disposal compared to other collaborations, SWME 10, 11 have performed the most 
complete study of discretization effects for this quantity so far. 

The determination of the renormalization factors merits special attention, owing to the 
fact that the renormalization and mixing patterns of four-quark operators such as Q AS=2 
are typically much more complicated than for quark bilinears. The renormalization factors 
used by HPQCD/UKQCD 06 are based on perturbation theory at one loop, which is by far 
the biggest source of systematic uncertainty quoted by these authors. The same is true for 
the SWME 10, 11 results [281, 283]: Here the estimated uncertainty in their latest estimate 
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Figure 11: Recent unquenched lattice results for Bk (MS-scheme, scale = 2GeV) and for 
the corresponding RGI parameter Bk- The gray band indicates our estimate. 



for Bk [281] arising from neglecting higher orders in perturbation theory amounts to 4.4%, 
which accounts for about 95% of the entire systematic error reported by this collaboration. 
The two other groups, RBC/UKQCD 07A,08, 10B and Aubin 09, have both implemented 
a non-perturbative renormalization procedure based on the RI/MOM scheme. Nonetheless, 
Aubin 09 state that their biggest single systematic uncertainty of 3.2% in Bk (or 85% of 
the total systematic error) is still associated with the uncertainty in the renormalization 
factor which links the i?-parameter of the bare operator to that in the MS-scheme at 2 GeV. 
RBC/UKQCD 10B [282] have investigated several different RI/MOM schemes, including so- 
called non-exceptional momenta. The resulting spread of results in different schemes is then 
included in the uncertainty for the renormalization factor. As observed in [282], this error can 
be significantly reduced if the conversion to the MS-scheme is performed at 3 GeV instead 
of the conventional choice, /i = 2 GeV. For this reason, the RBC/UKQCD collaboration 
prefers to quote B^ s (3 GeV), but also gives the value obtained at the conventional scale, for 
comparison with other results. All of the numbers listed in the first column of Table 16 refer 
to the conventional scale. In ref. [290] the RBC/UKQCD collaboration also investigated the 
effects of residual chiral symmetry breaking induced by the finite extent of the 5th dimension 
in the domain wall fermion formulation and found that the mixing of Q AS=2 with operators 
of opposite chirality was negligibly small. 

The rules of section 2 stipulate that our averages only involve results which are free of red 
tags and are published in a refereed journal. Papers which, at the time of writing, are still 
unpublished but are obvious updates of earlier published results should also be taken into 
account. In view of the above, we combine the results of Aubin 09 and RBC/UKQCD 10B, 
by applying the following procedure: in a first step statistical and systematic errors are added 
in quadrature, in order to have an overall error for the result of a particular collaboration. 
In a second step the results from the two groups are combined in a weighted average. This 
yields 

N f = 2 + 1: B K m (2GeV) = 0.536(17) B K = 0.738(20) . (83) 

Here we have used the RBC/UKQCD 10B result at the scale of 2 GeV, in spite of the bigger 
uncertainty reported for the associated renormalization factor. 
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Due to the use of one-loop perturbation theory for the renormalization factor, the SWME 
results are assigned a red tag and are thus not included in the average. However it must be 
stressed that they are in excellent agreement with eq. (83). Also the earlier HPQCD/UKQCD 
estimate, albeit with bigger statistical and systematic errors, agrees with our average. The 
observation that independent calculations produce compatible results - despite the fact that 
according to the colour coding in Table 16 some systematic effects are not fully controlled - 
lends further confidence to the result quoted in eq. (83). 

We now pass over to the results computed for Nf = 2 flavours of dynamical quarks. The 
UKQCD 04 result [288], computed with 0(a) improved Wilson fermions, is afflicted by many 
systematic errors and has a rather large overall error. Results by the JLQCD collaboration 
[280] are computed using overlap fermions in both the sea and valence quark sectors. Besides 
the usual systematic uncertainties listed in Table 16, they quote a 5% systematic error due 
to the uncertainty in the physical value of tq, which sets the scale in their simulation. They 
also quote a 1.4% systematic error, owing to the fact that they determine Bk on gauge 
configurations with a fixed value of the topological charge. In view of the high numerical 
cost of simulations with dynamical overlap quarks, it is not surprising that the lattice sizes of 
ref. [280] are smaller than those in present simulations with Nf = 2 + 1 flavours of dynamical 
staggered or domain wall fermions. Consequently, finite-volume effects are less well controlled 
in ref. [280]. The RBC 04 result [287] has been obtained using domain wall fermions at rather 
heavy pseudoscalar (pion) masses; otherwise it has similar systematics as the JLQCD 08 run. 
The result by ETM 09D [291] is based on a set of ensembles which allows for an extensive 
investigation of systematic uncertainties. In particular, it is so far the only Nf = 2 calculation 
involving data computed at three values of the lattice spacing. Being the only result without 
red tags, it is quoted as the currently best overall estimate for Nf = 2: 

N f = 2: B^(2GeV) = 0.516(18) (12) B K = 0.729(25)(17) . (84) 

It is clear that with the present level of accuracy, any dependence of Bx on the number of 
dynamical quark flavours, Nf, is not visible. 

It is interesting to compare these results to the best available quenched data, computed 
for large (valence) quark masses and with degenerate down and strange quarks. Avoiding a 
detailed discussion about the systematics, we show an (incomplete) compilation of indicative 
results in Tab. 17. We see that they agree with the ones obtained with Nf = 2 and Nf = 2 + 1 
at the level of two standard deviations. 



Collaboration 


Ref. 


N } 


B K 


B K 


ALPHA 09 


[292] 





0.532(25) 


0.73(3) 


CP-PACS 08 


[293] 





0.565(6) 


0.782(9) 


ALPHA 07 


[294] 





0.534(52) 


0.74(7) 


JLQCD 97 


[295] 





0.628(42) 


0.86(6) 



Table 17: Quenched results for the S-parameter B^ from various collaborations. Errors have 
been combined in quadrature. 
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A Glossary 

A.l Lattice actions 

In this appendix we give brief descriptions of the lattice actions used in the simulations and 
summarize their main features. 

A. 1.1 Gauge actions 

The simplest and most widely used discretization of the Yang-Mills part of the QCD action 
is the Wilson plaquette action [296] : 

^G^EEt 1 "^ TvW^\x)), (85) 

X fl<U 

where the plaquette, W^ 1 (x), is the product of link variables around an elementary square 
of the lattice, i.e. 

Ku\ x ) = U^(x)U u (x + afi)U^x + a£)- x U v {x)- x . (86) 

This expression reproduces the Euclidean Yang-Mills action in the continuum up to cor- 
rections of order a 2 . There is a general formalism, known as the "Symanzik improvement 
programme" [18, 19], which is designed to cancel the leading lattice artifacts, such that ob- 
servables have an accelerated rate of convergence to the continuum limit. The improvement 
programme is implemented by adding higher-dimensional operators, whose coefficients must 
be tuned appropriately in order to cancel the leading lattice artifacts. The effectiveness of 
this procedure depends largely on the method with which the coefficients are determined. 
The most widely applied methods (in ascending order of effectiveness) include perturbation 
theory, tadpole-improved (partially resummed) perturbation theory, renormalization group 
methods, and the non-perturbative evaluation of improvement conditions. 

In the case of Yang-Mills theory, the simplest version of an improved lattice action is 
obtained by adding rectangular 1x2 loops to the plaquette action, i.e. 

= 0£ U £ ( l ~ \ Re Tr W™{xj) + d £ (l - Tr W^\x)) | , (87) 

X y [1<V fl,V ) 
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where the coefficients cq , c\ satisfy the normalization condition cq + 8c\ = 1 . The Symanzik- 
improved [297], Iwasaki [298], and DBW2 [299, 300] actions are all defined through eq. (87) 
via particular choices for cq, c\. Details are listed in Table 18 together with the abbreviations 
used in the summary tables. 



Abbrev. 


Cl 


Description 


Wilson 





Wilson plaquette action 


tlSym 


-1/12 


tree-level Symanzik-improved gauge action 


tadSym 


variable 


tadpole Symanzik-improved gauge action 


Iwasaki 


-0.331 


Renormalization group improved ("Iwasaki") action 


DBW2 


-1.4088 


Renormalization group improved ("DBW2") action 



Table 18: Summary of lattice gauge actions. The leading lattice artifacts are 0(a 2 ) or better 
for all discretizations. 
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A. 1.2 Quark actions 

If one attempts to discretize the quark action, one is faced with the fermion doubling problem: 
the naive lattice transcription produces a 16-fold degeneracy of the fermion spectrum. 

Wilson fermions 

Wilson's solution to the doubling problem is based on adding a dimension-5 operator which 
removes the doublers from the low-energy spectrum. The Wilson-Dirac operator for the 
massless case reads [296] 

£>w = ^(V M + V*) + oV*V MJ (88) 

where V M , V* denote lattice versions of the covariant derivative. Adding the Wilson term, 
aV* V M , results in an explicit breaking of chiral symmetry even in the massless theory. Fur- 
thermore, the leading order lattice artifacts are of order a. With the help of the Symanzik im- 
provement programme, the leading artifacts can be cancelled by adding the so-called "Clover" 
or Sheikholeslami-Wohlert (SW) term. The resulting expression in the massless case reads 

^ sw — D w -\ — C sw <7 fiu-Ffxi/ , (89) 

where = an d F^ u is a lattice transcription of the gluon field strength tensor 

F^ u . Provided that the coefficient c sw is suitably tuned, observables computed using D sw 
will approach the continuum limit with a rate proportional to a?. Chiral symmetry remains 
broken, though. The coefficient c sw can be determined perturbatively at tree-level (tree-level 
impr., c sw = 1 or tlSW in short), via a mean field approach [301] (mean-field impr. or mfSW) 
or via a non-perturbative approach [302] (non-perturbativley impr. or npSW). 

Finally, we mention "twisted mass QCD" as a method which was originally designed to ad- 
dress another problem of Wilson's discretization: the Wilson-Dirac operator is not protected 
against the occurrence of unphysical zero modes, which manifest themselves as "exceptional" 
configurations. They occur with a certain frequency in numerical simulations with Wilson 
quarks and can lead to strong statistical fluctuations. The problem can be cured by intro- 
ducing a so-called "chirally twisted" mass term, after which the fermionic part of the QCD 
action in the continuum assumes the form [52] 

5 tm;cont = j fajtyfr^ + m + i^ T ^ x ). ( 90 ) 

Here, ^ q is the twisted mass parameter, and r 3 is a Pauli matrix. The standard action in 
the continuum can be recovered via a global chiral field rotation. The lattice action of twisted 
mass QCD (tmWil) for Nf = 2 flavours is defined as 

S^[U,^ V] = a 4 Yl V^XAv + ™o + ^ q 75T 3 )V(z)- (91) 

xeA B 

Although this formulation breaks physical parity and flavour symmetries, is has a number of 
advantages over standard Wilson fermions. In particular, the presence of the twisted mass pa- 
rameter fiq protects the discretized theory against unphysical zero modes. Another attractive 
feature of twisted mass lattice QCD is the fact that the leading lattice artifacts are of order 
a 2 without the need to add the Sheikholeslami-Wohlert term [53]. Although the problem of 
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explicit chiral symmetry breaking remains, the twisted formulation is particularly useful to 
circumvent some of the problems that are encountered in connection with the renormalization 
of local operators on the lattice, such as those required to determine £?k- 

Staggered fermions 

An alternative procedure to deal with the doubling problem is based on so-called "staggered" 
or Kogut-Susskind (KS) fermions [303]. Here the degeneracy is only lifted partially, from 
16 down to 4. It has become customary to refer to these residual doublers as "tastes" in 
order to distinguish them from physical flavours. At order a 2 different tastes can interact via 
gluon exchange, thereby generating large lattice artifacts. The improvement programme can 
be used to suppress taste-changing interactions, leading to "improved staggered fermions", 
with the so-called "Asqtad" and "HISQ" actions as the most widely used versions [304, 305]. 
The standard procedure to remove the residual doubling of staggered quarks ("four tastes per 
flavour" ) is to take fractional powers of the quark determinant in the QCD functional integral. 
This is usually referred to as the "fourth root trick" . The validity of this procedure has not 
been rigorously proven so far. In fact, it has been questioned by several authors, and the is- 
sue is still hotly debated (for both sides of the argument see the reviews in refs. [55-58, 306]). 

Ginsparg- Wilson fermions 

Fermionic lattice actions, which do not suffer from the doubling problem whilst preserving 
chiral symmetry go under the name of "Ginsparg- Wilson fermions". In the continuum the 
massless Dirac operator anti-commutes with 75. At non-zero lattice spacing chiral symmetry 
can be realized even if this condition is relaxed according to [307, 308] 

{D, 75 } = aD 75 D, (92) 

which is now known as the Ginsparg- Wilson relation [278]. A lattice Dirac operator which 
satisfies eq. (92) can be constructed in several ways. The "domain wall" construction proceeds 
by introducing a fifth dimension of length N§ and coupling the fermions to a mass defect 
(i.e. a negative mass term) [309]. The five-dimensional action can be constructed such that 
modes of opposite chirality are trapped at the four-dimensional boundaries in the limit of 
an infinite extent of the extra dimension [310]. In any real simulation, though, one has to 
work with a finite value of N$, so that the decoupling of chiral modes is not exact. This 
leads to a residual breaking of chiral symmetry, which, however, is exponentially suppressed. 
A doubler-free, (approximately) chirally symmetric quark action can thus be realized at the 
expense of simulating a five-dimensional theory. 

The so-called "overlap" or Neuberger-Dirac operator can be derived from the domain wall 
formulation [311]. It acts in four space-time dimensions and is, in its simplest form, defined 
by 

Dn = i(l--£=), A = l + s-aD w , a=-^—, (93) 

where D w is the massless Wilson-Dirac operator, and \s\ < 1 is a tunable parameter. The 
overlap operator removes all doublers from the spectrum, and can easily be shown to 
satisfy the Ginsparg- Wilson relation. The occurrence of an inverse square root in Z?n ren- 
ders the application of Z?n in a computer program potentially very costly, since it must be 
implemented using, for instance, a polynomial approximation. 

The third example of an operator which satisfies the Ginsparg- Wilson relation is the so- 
called fixed-point action [312, 313]. This construction proceeds via a renormalization group 
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approach. A related formalism are the so-called "chirally improved" fermions [314]. 
Smearing 

A simple modification which can help improve the action as well as the computational perfor- 
mance is the use of smeared gauge fields in the covariant derivatives of the fermionic action. 
Any smearing procedure is acceptable as long as it consists of only adding irrelevant (local) 
operators. Moreover, it can be combined with any discretization of the quark action. The 
"Asqtad" staggered quark action mentioned above [304] is an example which makes use of 
so-called "Asqtad" smeared (or "fat") links. Another example is the use of n-HYP smeared 
[289, 315], n-stout smeared [316, 317] or n-HEX (hypercubic stout) smeared [318] gauge links 
in the tree-level clover improved discretization of the quark action, denoted by "n-HYP tlSW" , 
"n-stout tlSW" and "n-HEX tlSW" in the following. 

In Table 19 we summarize the most widely used discretizations of the quark action and their 
main properties together with the abbreviations used in the summary tables. Note that 
in order to maintain the leading lattice artifacts of the actions as given in the table in non- 
spectral observables (like operator matrix elements) the corresponding non-spectral operators 
need to be improved as well. 



A. 2 Setting the scale 

In simulations of lattice QCD quantities such as hadron masses and decay constants are 
obtained in "lattice units" i.e. as dimensionless numbers. In order to convert them into 
physical units they must be expressed in terms of some experimentally known, dimensionful 
reference quantity Q. This procedure is called "setting the scale". It amounts to computing 
the non-perturbative relation between the bare gauge coupling go (which is an input parameter 
in any lattice simulation) and the lattice spacing a expressed in physical units. To this end 
one chooses a value for go an d computes the value of the reference quantity in a simulation: 
This yields the dimensionless combination, (aQ)\ go , at the chosen value of go- The calibration 
of the lattice spacing is then achieved via 

- M MeV]=«^a. (94) 

where Qlexp denotes the experimentally known value of the reference quantity. Common 
choices for Q are the mass of the nucleon, the $7 baryon or the decay constants of the pion 
and the kaon. Vector mesons, such as the p or K*-meson, are unstable and therefore their 
masses are not very well suited for setting the scale, despite the fact that they have been used 
over many years for that purpose. 

Another widely used quantity to set the scale is the hadronic radius ro, which can be 
determined from the force between static quarks via the relation [319] 

F{ro)rl = 1.65. (95) 

If the force is derived from potential models describing heavy quarkonia, the above relation 
determines the value of ro as ro ~ 0.5 fm. A variant of this procedure is obtained [320] by 
using the definition F{ri)r\ = 1.00, which yields r\ ~ 0.32 fm. It is important to realize that 
both ro and r\ are not directly accessible in experiment, so that their values derived from phe- 
nomenological potentials are necessarily model-dependent. Inspite of the inherent ambiguity 
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Abbrev. 


Discretization 


Leading lattice 
artifacts 


Chiral symmetry 


Remarks 


Wilson 


Wilson 


O(a) 


broken 




tmWil 


Twisted Mass Wilson 


0(a 2 ) at 
maximal twist 


broken 


flavour symmetry breaking: 
(M° s ) 2 - (M p ±) 2 ~ 0(a 2 ) 


tlSW 


Sheikholeslami-Wohlert 


0(g 2 a) 


broken 


tree- level impr., c sw = 1 


n-HYP 
tlSW 


Sheikholeslami-Wohlert 


0(g 2 a) 


broken 


tree-level impr., c sw = 1, 
n-HYP smeared gauge links 


n-stout 
tlSW 


Sheikholeslami-Wohlert 


0(g 2 a) 


broken 


tree-level impr., c sw = 1, 
n-stout smeared gauge links 


mfSW 


Sheikholeslami-Wohlert 


0(g 2 a) 


broken 


mean- field impr. 


npSW 


Sheikholeslami-Wohlert 


0(a 2 ) 


broken 


non-perturbatively impr. 


KS 


Staggered 


0(a 2 ) 


U(l) <g)U(l) subgr. 
unbroken 


rooting for Nf < 4 


Asqtad 


Staggered 


0(a 2 ) 


U(l) <g)U(l) subgr. 
unbroken 


Asqtad smeared gauge links, 
rooting for Nf < 4 


HISQ 


Staggered 


0(a 2 ) 


U(l) <g>U(l) subgr. 
unbroken 


HISQ smeared gauge links, 
rooting for Nf < 4 


DW 


Domain Wall 


asymptotically 
0(a 2 ) 


remnant breaking 
exponentially suppr. 


exact chiral symmetry and 
0(a) impr. only in the limit 
L s — > oo 


overlap 


Neuberger 


0(a 2 ) 


exact 





Table 19: The most widely used discretizations of the quark action and some of their proper- 
ties. Note that in order to maintain the leading lattice artifacts of the action in non-spectral 
observables (like operator matrix elements) the corresponding non-spectral operators need to 
be improved as well. 
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whenever hadronic radii are used to calibrate the lattice spacing, they are very useful quanti- 
ties for performing scaling tests and continuum extrapolations of lattice data. Furthermore, 
they can be easily computed with good statistical accuracy in lattice simulations. 

A. 3 Matching and running 

The lattice formulation of QCD amounts to introducing a particular regularization scheme. 
Thus, in order to be useful for phenomenology, hadronic matrix elements computed in lattice 
simulations must be related to some continuum reference scheme, such as the MS-scheme of 
dimensional regularization. The matching to the continuum scheme usually involves running 
to some reference scale using the renormalization group. 

In principle, the matching factors which relate lattice matrix elements to the MS-scheme, 
can be computed in perturbation theory formulated in terms of the bare coupling. It has been 
known for a long time, though, that the perturbative expansion is not under good control. 
Several techniques have been developed which allow for a non-perturbative matching between 
lattice regularization and continuum schemes, and are briefly introduced here. 

Regularization-independent Momentum Subtraction 

In the Regularization-independent Momentum Subtraction ("RI/MOM" or "RI") scheme 
[115] a non-perturbative renormalization condition is formulated in terms of Green functions 
involving quark states in a fixed gauge (usually Landau gauge) at non-zero virtuality. In this 
way one relates operators in lattice regularization non-perturbatively to the RI scheme. In a 
second step one matches the operator in the RI scheme to its counterpart in the MS-scheme. 
The advantage of this procedure is that the latter relation involves perturbation theory for- 
mulated in the continuum theory. The uncontrolled use of lattice perturbation theory can 
thus be avoided. A technical complication is associated with the accessible momentum scales 
(i.e. virtualities) , which must be large enough (typically several GeV) in order for the per- 
turbative relation to MS to be reliable. The momentum scales in simulations must stay well 
below the cutoff scale (i.e. 2ir over the lattice spacing), since otherwise large lattice artifacts 
are incurred. Thus, the applicability of the RI scheme traditionally relies on the existence of 
a "window" of momentum scales, which satisfy 

Aqcd < P < 27TG- 1 . (96) 

However, solutions for mitigating this limitation, which involve continuum limit, non-perturbative 
running to higher scales in the RI/MOM scheme, have recently been proposed and imple- 
mented [65, 105, 282, 321]. 
Schrddinger functional 

Another example of a non-perturbative matching procedure is provided by the Schrodinger 
functional (SF) scheme [114]. It is based on the formulation of QCD in a finite volume. If 
all quark masses are set to zero the box length remains the only scale in the theory, such 
that observables like the coupling constant run with the box size L. The great advantage is 
that the RG running of scale-dependent quantities can be computed non-perturbatively using 
recursive finite-size scaling techniques. It is thus possible to run non-perturbatively up to 
scales of, say, 100 GeV, where one is sure that the perturbative relation between the SF and 
MS-schemes is controlled. 
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Perturbation theory 

The third matching procedure is based on perturbation theory in which higher order are effec- 
tively resummed [301]. Although this procedure is easier to implement, it is hard to estimate 
the uncertainty associated with it. 



In Table 20 we list the abbreviations used in the compilation of results together with a short 
description. 



Abbrev. 


Description 


RJ 


Regularization-independent momentum subtraction scheme 


SF 


Schrodinger functional scheme 


PTU 


matching/running computed in perturbation theory at one loop 


PT2£ 


matching/running computed in perturbation theory at two loop 



Table 20: The most widely used matching and running techniques. 



A. 4 Chiral extrapolation 

As mentioned in section 1.1, Symanzik's framework can be combined with Chiral Pertur- 
bation Theory. The well-known terms occurring in the chiral effective Lagrangian are then 
supplemented by contributions proportional to powers of the lattice spacing a. The additional 
terms are constrained by the symmetries of the lattice action and therefore depend on the 
specific choice of the discretization. The resulting effective theory can be used to analyze 
the a-dependence of the various quantities of interest - provided the quark masses and the 
momenta considered are in the range where the truncated chiral perturbation series yields an 
adequate approximation. Understanding the dependence on the lattice spacing is of central 
importance for a controlled extrapolation to the continuum limit. 

For staggered fermions, this program has first been carried out for a single staggered flavor 
(a single staggered field) [24] at 0(a 2 ). In the following, this effective theory is denoted by 
SxPT. It was later generalized to an arbitrary number of flavours [28, 30], and to next-to- 
leading order [33]. The corresponding theory is commonly called Rooted Staggered chiral 
perturbation theory and is denoted by RSxPT. 

For Wilson fermions, the effective theory has been developed in [23, 26, 31] and is called 
WxPT, while the theory for Wilson twisted mass fermions [32, 34, 40] is termed tmWxPT. 

Another important approach is to consider theories in which the valence and sea quark 
masses are chosen to be different. These theories are called partially quenched. The acronym 
for the corresponding chiral effective theory is PQxPT [20-22, 25]. 

Finally, one can also consider theories where the fermion discretizations used for the sea 
and the valence quarks are different. The effective chiral theories for these "mixed action" 
theories are referred to as MAxPT [27, 29, 35-39]. 
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A. 5 Summary of simulated lattice actions 

In the following two tables we summarize the gauge and quark actions used in the various 
calculations. Abbreviations are explained in section A. 1.1 and A. 1.2, and summarized in 
tables 18 and 19. 
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Collab. 


Ref. 


Nf 


gauge 
action 


quark 
action 


ALPHA 05 


[98] 


2 


Wilson 


npSW 


CERN-TOV 06 


[254] 


2 


Wilson 


Wilson/npSW 


CERN 08 


[215] 


2 


Wilson 


npSW 


Bernardoni 10 


[244] 


2 


Wilson 


npSW t 


CP-PACS 01 


[101] 


2 


Iwasaki 


mfSW 


ETM 07, 07A, 08, 09, 09A-D, 10B, D 


[49, 54, 94, 160, 
161, 177, 235, 
237, 246, 291] 


2 


tlSym 


tmWil 


HHS 08 


[247] 


2 


tadSym 


n-HYP tlSW 


JLQCD 08 


[280] 


2 


Iwasaki 


overlap 


JLQCD 02, 05 


[100, 164] 


2 


Wilson 


npSW 


JLQCD/TWQCD 07, 08A, 10 


[95, 240, 248] 


2 


Iwasaki 


overlap 


QCDSF 07 


[162] 


2 


Wilson 


npSW 


QCDSF/UKQCD 04, 06, 06A, 07 


[96, 99, 178, 260] 


2 


Wilson 


npSW 


RBC 07 


[73] 


2 


DBW2 


DW 


RBC 04, 06 


[163, 287] 


2 


DBW2 


DW 


SPQcdR 05 


[97] 


2 


Wilson 


Wilson 


UKQCD 04, 07 


[159, 288] 


2 


Wilson 


npSW 



The calculation uses overlap fermions in the valence quark sector. 



Table 21: Summary of simulated lattice actions with Nf = 2 quark flavours. 



Collab. 


Ref. 


N f 


gauge 
action 


quark 
action 


Aubin 08, 09 


[173, 269] 


2 + 1 


tadSym 


Asqtad ' 


SWME 10, 11 


[281, 283] 


2 + 1 


tadSym 


Asqtad + 


Blum 10 


[74] 


2 + 1 


Iwasaki 


DW 


BMW 10A-C 


[65, 75, 105] 


2 + 1 


tlSym 


2-HEX tlSW 


BMW 10 


[41] 


2 + 1 


tlSym 


6-stout tlSW 


CP-PACS/JLQCD 07 


[109] 


2 + 1 


Iwasaki 


npSW 


ETM 10, 10E 


[60, 170] 


2 + 1 + 1 


Iwasaki 


tmWil 


HPQCD 05 


[110] 


2 + 1 


tadSym 


Asqtad 


HPQCD/UKQCD 06 


[285] 


2 + 1 


tadSym 


Asqtad 


HPQCD/UKQCD 07 


[175] 


2 + 1 


tadSym 


Asqtad * 


HPQCD/MILC/UKQCD 04 


[111] 


2 + 1 


tadSym 


Asqtad 


JLQCD 09 


[241] 


2 + 1 


Iwasaki 


overlap 


JLQCD/TWQCD 08B, 09A 


[172, 243] 


2 + 1 


Iwasaki 


overlap 


JLQCD/TWQCD 10 


[240] 


2 + 1,3 


Iwasaki 


overlap 


LHP 04 


[259] 


2 + 1 


tadSym 


Asqtad ' 


MILC 04, 07, 09, 09A, 10, 10A 


[6, 77, 103, 111, 171, 322] 


2 + 1 


tadSym 


Asqtad 


NPLQCD 06 


[176] 


2 + 1 


tadSym 


Asqtad ' 


PACS-CS 08, 08A, 09, 10 


[42, 63, 64, 174] 


2 + 1 


Iwasaki 


npSW 


RBC/UKQCD 07, 08, 08A, 10, 10A-B[106, 108, 158, 250, 282, 284] 


2 + 1 


Iwasaki 


DW 


TWQCD 08 


[242] 


2 + 1 


Iwasaki 


DW 



' The calculation uses domain wall fermions in the valence quark sector. 

* The calculation uses HISQ staggered fermions in the valence quark sector. 

+ The calculation uses HYP smeared improved staggered fermions in the valence quark sector. 

Table 22: Summary of simulated lattice actions with Nf = 2 + 1 or Nf = 2 + 1 + 1 quark 
flavours. 
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B Notes 

B.l Notes to section 3 on quark masses 



Collab. 


Ref. 


N f 


a [fm] 


Description 


PACS-CS 10 


[64] 


2+1 


0.09 


cf. PACS-CS 08 


M1LO 1UA 


[103] 


2+1 




CI. 1V11LU 09, 09A 


BMW 10A, 10B 


[65, 105] 


2+1 


0.116,0.093,0.077, 

U.U0o,U.Uo4 


Scale setting via M n , Mk, Mn. 


RBC/UKQCD 10A 


[106] 


2+1 


0.114, 0.087 


Scale set through Mq. 


Blum 10 


[74] 


2+1 


0.11 


Relies on RBC/UKQCD 08 scale settir 


PACS-CS 09 


[42] 


2+1 


0.09 


Scale setting via Mq. 


HPQCD 09, 10 


[104, 107] 


2+1 






MILC 09A, 09 


[6, 59] 


2+1 


0.045, 0.06, 0.09 


Scale set through ri and T and contin- 
uum extrapolation based on RSxPT. 


PACS-CS 08 


[63] 


2+1 


0.09 


Scale set through Mn. Non- 
perturbatively 0(a)-improved. 


RBC/UKQCD 08 


[108] 


2+1 


0.11 


Scale set through Mq. Automatic 
0(a)-improvement due to appoximate 
chiral symmetry. (Aqcdij) 2 ~ 4% sys- 
tematic error due to lattice artifacts 
added. 


CP-PACS/JLQCD 07 


[109] 


2+1 


0.07,0.10,0.12 


Scale set through Mk or M$. Non- 
perturbatively 0(a)-improved. 


HPQCD 05 


[110] 


2+1 


0.09,0.12 


Scale set through the T — T' mass dif- 
ference. 


HPQCD/MILC/UKQCD 04, 
MILC 04 


[77, 111] 


2+1 


0.09,0.12 


Scale set through ri and T and contin- 
uum extrapolation based on RSxPT. 



Table 23: Continuum extrapolations/estimation of lattice artifacts in determinations of m U( i 
and m s with Nf = 2 + 1 quark flavours. 
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Collab. 


Ref. 


Nf 


a [fm] 


Description 


ETM 10B 


[94] 


2 


0.098, 0.085, 
0.067, 0.054 


Scale set through Fw 


JLQCD/TWQCD 08 A 


[95] 


2 


0.12 


Scale set through ro. 


RBC 07 


[73] 


2 


0.12 


Scale set through M p . 


ETM 07 


[49] 


2 


0.09 


Scale set through F w . 


QCDSF/UKQCD 06 


[96] 


2 


0.065-0.09 


Scale set through ro. 


SPQcdR 05 


[97] 


2 


0.06,0.08 


Scale set through Mk* ■ 


ALPHA 05 


[98] 


2 


0.07-0.12 


Scale set through ro. 


QCDSF/UKQCD 04 


[99] 


2 


0.07-0.12 


Scale set through ro. 


JLQCD 02 


[100] 


2 


0.09 


Scale set through M p . 


CP-PACS 01 


[101] 


2 


0.11,0.16,0.22 


Scale set through M p . 



Table 24: Continuum extrapolations/estimation of lattice artifacts in determinations of m u d 
and m s with Nf = 2 quark flavours. 
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Collab. 


Ref. 


Nf 


M., min [MeV] 


Description 


PACS-CS 10 


[64] 


2+1 




cf. PACS-CS 08 


MILC 10A 


[103] 


2+1 




NLO SU(2) SxPT. Cf. also MILC 09A,09. 


BMW 10A, 10B 


[65, 105] 


2+1 


135 


Interpolation to the physical point. 


RBC/UKQCD 10A 


[1061 

L J 


2+1 


290 




Blum 10 


[74, 108] 


2+1 


242 (valence), 
330 (sea) 


Extrapolation done on the basis of PQxPT 
formulae with virtual photons. 


PACS-CS 09 


[421 

L J 


2+1 


135 


Physical point reached by reweighting 
technique, no chiral extrapolation needed. 


HPQCD 09, 10 


[104, 107] 


2+1 






MILC 09A, 09 


[6, 59] 


2+1 


177, 240 


NLO SU(3) RSxPT, continuum xPT 
at NNLO and NNNLO and NNNNLO 
analytic terms. The lightest Nambu- 
Goldstone mass is 177 MeV (09A) and 224 
MeV (09) (at a =0.09fm) and the lightest 
RMS mass is 258MeV (at a =0.06fm). 


PACS-CS 08 


[63] 


2+1 


156 


NLO SU(2) X PT and SU(3) (Wilson)xPT. 


RBC/UKQCD 08 


[108] 


2+1 


242 (valence), 
330 (sea) 


SU(3) PQxPT and heavy kaon NLO SU(2) 
PQxPT fits. 


CP-PACS/JLQCD 07 


[109] 


2+1 


620 


NLO Wilson x?T fits to meson masses. 


HPQCD 05 


[110] 


2+1 


240 


PQ RSxPT fits. 


HPQCD/MILC/UKQCD 04, 
MILC 04 


[77, 111] 


2+1 


240 


PQ RSxPT fits. 



Table 25: Chiral extrapolation/minimum pion mass in determinations of m u d and m s with 
Nf = 2 + 1 quark flavours. 



84 



Collab. 


Ref. 


N f 


A^, min [MeV] 


Description 


ETM 10B 


[94] 


2 


270 


Fits based on NLO x?T and Symanzik ex-, 
pansion up to 0(a 2 ) 


JLQCD/TWQCD 08 A 


[95] 


2 


290 


NLO xPT fits. 


RBC 07 


[73] 


2 


440 


NLO fit including O(a) effects. 


ETM 07 


[49] 


2 


300 


Polynomial and PQxPT fits. 


QCDSF/UKQCD 06 


[96] 


2 


520 (valence), 
620 (sea) 


NLO (PQ)xPT fits. 


SPQcdR 05 


[97] 


2 


600 


Polynomial fit. 


ALPHA 05 


[98] 


2 


560 


LO xPT fit. 


QCDSF/UKQCD 04 


[99] 


2 


520 (valence), 
620 (sea) 


NLO (PQ)xPT fits. 


JLQCD 02 


[100] 


2 


560 


Polynomial and x?T fits. 


CP-PACS 01 


[101] 


2 


430 


Polynomial fits. 



Table 26: Chiral extrapolation/minimum pion mass in determinations of m u d and m s with 
Nf = 2 quark flavours. 
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Collab. 


Ref. 


N f 


L [fm] 




Description 


PACS-CS 10 


[64] 


2+1 






cf. PACS-CS 08 


MILC 10A 


[103] 


2+1 






cf. MILC 09A,09 


BMW 10A, 10B 


[65, 105] 


2+1 


> 5.0 


> 4.0 


FS corrections below 5 per mil 
on the largest lattices. 


RBC/UKQCD 10A 


[106] 


2+1 


2.7 


£ 4.0 




Blum 10 


[74] 


2+1 


1.8, 2.7 




Simulations done with dy- 
namical photons; large fi- 
nite volume effects analyti- 
cally corrected for, but not re- 
lated to M n L. 


PACS-CS 09 


[42] 


2+1 


2.9 


2.0 


Only one volume. 


HPQCD 09, 10 


[104, 107] 


2+1 








MILC 09A, 09 


[6 591 


2+1 


2.5, 2.9, 3.4, 

3.6, 3.8, 5.8 


4.1, 3.8 




PACS-CS 08 


[63] 


2+1 


2.9 


2.3 


Correction for FSE from X PT 
using [323]. 


RBC/UKQCD 08 


[108] 


2+1 


1.8, 2.7 


4.6 


Various volumes for compar- 
ison and correction for FSE 
from xPT [194, 195, 323]. 


CP-PACS/JLQCD 07 


[109] 


2+1 


2.0 


6.0 


Estimate based on the com- 
parison to a L = 1.6 fm vol- 
ume assuming powerlike de- 
pendence on L. 


HPQCD 05 


[110] 


2+1 


2.4, 2.9 


3.5 




HPQCD/MILC/UKQCD 04, 
MILC 04 


[77, 111] 


2+1 


2.4, 2.9 


3.5 


NLO SxPT. 



Table 27: Finite volume effects in determinations of m uc i and m s with Nf = 2 + 1 quark 
flavours. 
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Collab. 


Ref. 


N f 


L [fm] 




Description 


ETM 10B 


[94] 


2 


> 2.0 


3.5 


One volume L — 1.7 fm at 
rriTr — 495, a = 0.054fm. 


JLQCD/TWQCD 08 A 


[95] 


2 


1.9 


2.8 


Corrections for FSE based on 
NLO xPT. 


RBC 07 


[73] 


2 


1.9 


4.3 


Estimate of FSE based on a 
model. 


ETM 07 


[49] 


2 


2.1 


3.2 


NLO PQyPT 


QCDSF/UKQCD 06 


[96] 


2 


1.4-1.9 


4.7 




SPQcdR 05 


[97] 


2 


1.0-1.5 


4.3 


Comparison between 1.0 and 
1.5 fm. 


ALPHA 05 


[98] 


2 


2.6 


7.4 




QCDSF/UKQCD 04 


[99] 


2 


1.7-2.0 


4.7 




JLQCD 02 


[100] 


2 


1.8 


5.1 


Numerical study with three 
volumes. 


CP-PACS 01 


[101] 


2 


2.0-2.6 


5.7 





Table 28: Finite volume effects in determinations of m U( i and m s with Nf = 2 quark flavours. 
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Ooilab. 


Ker. 


N f 


Description 


PACS-CS 10 


[64] 


2+1 


Non-perturbative renormalization and running; 
Schrodinger functional method. 


MILC 10A 


[103] 


2+1 


cf. MILC 09A,09 


BMW 10A, 10B 


[65, 105] 


2+1 


Non-perturbative renormalization (tree-level im- 
proved RI-MOM), non-perturbative running. 


RBC/UKQCD 10A 


[106] 


2+1 


Non-perturbative renormalization (RI/SMOM). 


Blum 10 


[74] 


2+1 


Relies on non-perturbative renormalization fac- 
tors calculated by RBC/UKQCD 08; no QED 
renormalization. 


PACS-CS 09 


[42] 


2+1 


Non-perturbative renormalization; Schrodinger 
functional method. 


HPQCD 09, 10 


[104, 107] 


2+1 


Lattice calculation of m s /m c : m a derived from 
the perturbative determination of m c . 


MILC 09A, 09 


[6, 59] 


2+1 


2-loop perturbative renormalization. 


PACS-CS 08 


[63] 


2+1 


1-loop perturbative renormalization. 


RBC/UKQCD 08 


[108] 


2+1 


Non-perturbative renormalization, 3-loop pertur- 
bative matching. 


CP-PACS/JLQCD 07 


[109] 


2+1 


1-loop perturbative renormalization, tadpole im- 
proved. 


HPQCD 05 


[110] 


2+1 


2-loop perturbative renormalization. 


HPQCD/MILC/UKQCD 04, 
MILC 04 


[77, 111] 


2+1 


1-loop perturbative renormalization. 



Table 29: Renormalization in determinations of m u d and m s with Nf = 2 + 1 quark flavours. 
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Collab. 


Ref. 


Nf 


Description 


ETM 10B 


[94] 


2 


Non-perturbative renormalization. 


JLQCD/TWQCD 08 A 


[95] 


2 


Non-perturbative renormalization. 


RBC 07 


[73] 


2 


Non-perturbative renormalization. 


ETM 07 


[49] 


2 


Non-perturbative renormalization. 


QCDSF/UKQCD 06 


[96] 


2 


Non-perturbative renormalization. 


SPQcdR 05 


[97] 


2 


Non-perturbative renormalization. 


ALPHA 05 


[98] 


2 


Non-perturbative renormalization. 


QCDSF/UKQCD 04 


[99] 


2 


Non-perturbative renormalization. 


JLQCD 02 


[100] 


2 


1-loop perturbative renormalization. 


CP-PACS 01 


[101] 


2 


1-loop perturbative renormalization. 



Table 30: Renormalization in determinations of m U( i and m s with Nf = 2 quark flavours. 
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B.2 Notes to section 4 on \V u d\ and \V V 



Collab. 


Ref. 


N } 


a [fm] 


Description 


RBC/UKQCD 07,10[158, 159] 


2+1 


0.114(2) 


Scale fixed through Q baryon mass. Add 
(Aqcdo) 2 ~ 4% systematic error for lat- 
tice artifacts. Fifth dimension with exten- 
sion L s = 16, therefore small residual chi- 
ral symmetry breaking and approximate 
O(a)- imrovement . 


ETM 10D 


[160] 


2 


0.054, 0.068, 
0.086, 101 


Scale set through F % . Automatic O(a) 
impr., flavour symmetry breaking: 
(M° s ) 2 -(M± s ) 2 ~0(a 2 ). 


ETM 09A 


[161] 


2 


0.069, 0.883(6), 0.103 


Scale set through F n . Automatic 0(a) 
impr., flavour symmetry breaking: 
(M° PS f - (Mp s f ~ 0(a 2 ). 
Three lattice spacings only for pion mass 
470 MeV. 


QCDSF 07 


[162] 


2 


0.075 


Scale set with ro- Non-perturbatively 
0(a)-improved Wilson fermions, not clear 
whether currents improved. 


RBC 06 


[163] 


2 


0.12 


Scale set through M p . Automatic 0(a)- 
improvement due to approximate chiral 
symmetry of the action. 


JLQCD 05 


[164] 


2 


0.0887 


Scale set through M p . Non-perturbatively 
0(a)-improved Wilson fermions. 



Table 31: Continuum extrapolations/estimation of lattice artifacts in determinations of /+(0). 



90 



Collab. 


Ref. 


N f 


M^min [MeV] 


Description 


RBC/UKQCD 07,10 


[158, 159] 


2+1 


330 


NLO SU(3) x"T with phenomenological 
ansatz for higher orders. 


ETM 10D 


[160] 


2 


260 


NLO heavy kaon SU(2) x PT and NL 
SU(3) xPT and phenomenological ansatz 
for higher orders. Average of /+(0)-fit and 
joint /+(0)-/WA-nt. 


ETM 09A 


[161] 


2 


260 


NLO heavy kaon SU(2) xPT and NLO 
SU(3) xPT an d phenomenological ansatz 
for higher orders. 


QCDSF 07 


[162] 


2 


591 


Only one value for the pion mass. 


RBC 06 


[163] 


2 


490 


NLO SU(3) x?T and phenomenological 
ansatz for higher orders. 


JLQCD 05 


[164] 


2 


550 


NLO SU(3) x?T and phenomenological 
ansatz for higher orders. 



Table 32: Chiral extrapolation/minimum pion mass in determinations of /+(0). 



Collab. 


Ref. 


Nf 


L [fm] 




Description 


RBC/UKQCD 07,10 


[158, 159] 


2+1 


2.7 


4.7 


Two volumes for all but the lightest 
pion mass. 


ETM 10D 


[160] 


2 


1.7-2.8 


3.7 




ETM 09A 


[161] 


2 


2.1, 2.8 


3.2 


Two volumes at M n = 300 MeV 
and xPT-motivated estimate of the 
error due to FSE. 


QCDSF 07 


[162] 


2 


1.9 


5.4 




RBC 06 


[163] 


2 


1.9 


4.7 




JLQCD 05 


[164] 


2 


1.8 


4.9 





Table 33: Finite volume effects in determinations of /+(0). 
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Collab. 


Ref. 


N f 


a [fm] 


Description 


ETM 10E 


[170] 


2+1+1 


0.061, 0.078 


Scale set through f-ir/m^. Two lattice 
spacings but a-dependence ignored in all 
fits. Finer lattice spacing from [249]. 


RBC/UKQCD 10A 


[106] 


2+1 


0.114, 0.087 


Scale set through Mq. 


MILC 10 


[171] 


2+1 


0.045 - 0.09 


3 lattice spacings, continuum extrapola- 
tion by means of RS^PT. 


BMW 10 


[41] 


2+1 


0.065, 0.083, 
0.124 


Scale set through Mo =. Perturbative 
O ( a ) - improvement . 


JLQCD/TWQCD 08 A 


[95] 


2+1 


0.100(2) 


Scale set through F n . Automatic 0(a)- 
improvement due to chiral symmetry of ac- 
tion. 


MILC 09A 


[59] 


2+1 


0.045, 0.06, 0.09 


Scale set through n and T and continuum 
extrapolation based on RSxPT. 


MILC 09 


[6] 


2+1 


0.045, 0.06, 
0.09, 0.12 


Scale set through n and T and continuum 
extrapolation based on R»xP 1 • 


Aubin 08 


[173] 


2+1 


0.09, 0.12 


Scale set through n and T and continuum 
extrapolation based on MAxPT. 


PACS-CS 08, 08A 


[63, 174] 


2+1 


0.0907(13) 


Scale set through Mq. Non-perturbatively 
0(a)-improved. 


HPQCD/UKQCD 07 


[175] 


2+1 


0.06, 0.09, 0.15 


Scale set through n and T and contin- 
uum extrapolation on continuum-xPT mo- 
tivated ansatz. Taste breaking of sea 
quarks ignored. 


RBC/UKQCD 08 


[108] 


2+1 


0.114(2) 


Scale set through Afo. Automatic 0(a)- 
improvement due to appoximate chiral 
symmetry. (Aqcdoi) 2 ~ 4% systematic er- 
ror due to lattice artifacts added. 


NPLQCD 06 


[176] 


2+1 


0.125 


Scale set through ro and F n . Taste break- 
ing of sea quarks ignored. 


ETM 10D 


[160] 


2 


0.054, 0.068, 
0.086, 0.101 


Scale set through f+ Automatic 
0(a) impr., flavour symmetry breaking: 
(M° PS ) 2 - (M±„) 2 ~ 0(a 2 ). 


ETM 09 


[177] 


2 


0.07, 0.09, 0.10 


Scale set through F w . Automatic O(a) 
impr., flavour symmetry breaking: 
(M° s ) 2 - (M± s ) 2 ~ 0(a 2 ). 


QCDSF/UKQCD 07 


[178] 


2 


0.06, 0.07 


Scale set through F v . Non-perturbativc 
0(a)- improvement . 



Table 34: Continuum extrapolations/estimation of lattice artifacts in determinations of 

fK/f*. 
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Oollab. 


Ker. 


A/ 


MTr.inln [MeVJ 


Description 


ETM 10E 


[170] 


2+1+1 


265 


Pion mass on finer lattice from [324]. 


J.V/-L>V^/ V • I \ V, / V.. 1 J \_ UA 


[1 Ofil 

[1UUJ 


9-1- 1 


290 


Ti pci i It c aTP \~t!A tipri /~in n pnw nvi 1\TT <(l >T}( | 

J-VjCoLIILo CLlC UCloCLl Ull llCdv V JYCll_H± ll Uvy O \ iJ J 

PQ X PT. 


MILC 10 


[171] 


2+1 


258 


Lightest Nambu-Goldstone mass is 177 MeV (at 
0.09 fm) and lightest RMS mass is 258 MeV (at 
0.06 fm). NLO rSxPT and NNLO xPT. 


BMW 10 


[41] 


2+1 


190 


Comparison of various fit-ansatze: SU(3) xPT, 
heavy kaon SU(2) xPT, polynomial. 


JLQCD/TWQCD 08 A 


[95] 


2+1 


340 


NNLO SU(3) xPT. 


MILC 09A 


[59] 


2+1 


177 


NLO SU(3) RSxPT, continuum xPT at NNLO 
and NNNLO and NNNNLO analytic terms and 
heavy kaon SU(2) RSxPT with NNLO contin- 
uum chiral logs on a smaller sub-set of the lat- 
tices. The lightest Nambu-Goldstone mass is 
177 MeV (at a = 0.09 fm) and the lightest RMS 
mass is 258 MeV (at a = 0.06 fm). 


MILC 09 


[6] 


2+1 


240 


NLO SU(3) RSxPT with continuum x?T 
NNLO and NNNLO analytic terms added. 
According to [59] the lightest sea Nambu- 
Goldstone mass is 224 MeV and the lightest 
RMS mass is 258 MeV (at a = 0.06 fm). 


Aubin 08 


[173] 


2+1 


240 


NLO MAxPT. According to [59] the lightest sea 
Nambu-Goldstone mass is 246 MeV (at a = 0.09 
fm) and the lightest RMS mass is 329 MeV (at 
a = 0.09 fm). 


PACS-CS 08, 08A 


[63, 174] 


2+1 


156 


NLO SU(2) xPT and SU(3) (Wilson) X PT. 


HPCD/UKQCD 07 


[175] 


2+1 


240 


NLO SU(3) chiral perturbation theory with 
NNLO and NNNLO analytic terms. The sea 
RMS mass for the employed lattices is heavier 
than 250 MeV. 


RBC/UKQCD 08 


[108] 


2+1 


330 


While SU(3) PQxPT fits were studied, final 
results are based on heavy kaon NLO SU(2) 
PQ X PT. 


IN r LiV^OJJ UD 


[176J 




oUU 


1N.HJ ou^oj \r i ana some ininj^u terms, ine 
sea RMS mass for the employed lattices is heav- 
ier. 


ETM 10D 


[160] 


2 


260 


NLO SU(3) xPT arm phenomenological ansatz 
for higher orders. Joint /+(0)-/i<-// T -fit. 


ETM 09 


[177] 


2 


270 


NLO heavy meson SU(2) X PT and NLO SU(3) 
X PT. 


QCDSF/UKQCD 07 


[178] 


2 


300 


Linear extrapolation of lattice data. 



Table 35: Chiral extrapolation/minimum pion mass in determinations of fx/ fit. 
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Collab. 


Ref. 


N f 


L [fm] 




Description 


ETM 10E 


[170] 


2+1+1 


1.9 - 2.9 




Simulation parameters from [249, 
324]. 


RBC/UKQCD 10A 


[106] 


2+1 


2.7 


£ 4.0 




MILC 10 


[171] 


2+1 


2.52 


4.11 


L> 2.9 fm for the lighter masses. 


BMW 10 


[41] 


2+1 


2.0, 2.6, 2.7, 
3.0, 4.0 


4.0 


Various volumes for comparison 
and correction for FSE from x?T 

using [OZOJ. 


JLQCD/TWQCD 08 A 


[95] 


2+1 


1.6 


2.8 


Estimate of FSE using x?T [323, 
ozoj 


MILC 09A 


[59] 


2+1 


2.5, 2.9, 3.4, 

3.6, 3.8, 5.8 


4.1 




MILC 09 


[6] 


2+1 


2.4, 2.5, 2.9, 
3.0, 3.4, 3.6, 
o.o, o.o 


3.8 


Various volumes for comparison 
and correction for FSEs from 


Aubin 08 


[173] 


2+1 


2.4, 3.4 


4.0 


Correction for FSE from MA X PT. 


PACS-CS 08, 08A 


[63, 174] 


2+1 


2.9 


2.3 


Correction for FSE from X PT using 
[323]. 


HPCD/UKQCD 07 


[175] 


2+1 


2.3, 2.4, 2.5, 
2.9 


3.8 


Correction for FSE from x?T using 
[323]. 


RBC/UKQCD 08 


[108] 


2+1 


1.8, 2.7 


4.6 


Various volumes for comparison 
and correction for FSE from xPT 
[194, 195, 323]. 


NPLQCD 06 


[176] 


2+1 


2.5 


3.8 


Correction for FSE from SxPT [28, 
30] 


ETM 10D 


[160] 


2 


2.0-2.5 


3.7 




ETM 09 


[177] 


2 


2.0, 2.4, 2.7 


3.3 


Correction for FSE from xPT [194, 
195, 323]. 


QCDSF/UKQCD 07 


[178] 


2 


1.4,..., 2.6 


4.2 


Correction for FSE from X PT 



Table 36: Finite volume effects in determinations of fx/f^- 
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B.3 Notes to section 5 on Low-Energy Constants 



Collab. 


Ref. 


Nf 


a [fm] 


Description 


Bernardoni 10 


[244] 


2 


0.0784(10) 


Scale fixed through Mk- Non-perturbative 0(a) im- 
provement. No estimate of systematic error. 


JLQCD/TWQCD 10 


[240] 


2 


0.11 


One lattice spacing, scale fixed through ro. Estimate 
of 7% systematic error due to finite lattice spacing. 


JLQCD/TWQCD 09 


[252] 


2 


0.1184(21) 


Automatic 0(a) impr., exact chiral symmetry. Scale 
fixed through ro. Scale ambiguity added as system- 
atic error. 


ETM 09B 


[246] 


2 


0.063, 0.073 


Automatic O(a) impr. Estimate of lattice artifacts. 
(4% — 10% effect, not added in final result) by com- 
paring two lattice spacings. ro = 0.49 fm used to 
convert results in physical units. 


ETM 09C 


[2371 


2 


051 - 1 


Automatic 0(&) impr. Scale fixed through F^- 4 
lattice spacings, continuum extrapolation. 


ETM 08 


[235] 


2 


0.07-0.09 


Automatic O(tt) impr. Discretization effects esti- 
mated by comparing two lattice spacings. Scale fixed 
through Fk- 


JLQCD/TWQCD 08 A 
JLQCD 08A 


[95] 
[264] 


2 


0.1184(3)(21) 


Automatic O(a) impr., exact chiral symmetry. Scale 
fixed through ro. Scale ambiguity added as system- 
atic error. 


CERN 08 


[215] 


2 


0.0784(10) 


Scale fixed through Mk- Non-perturbative 0(a) im- 
provement. No estimate of systematic error. 


HHS 08 


[247] 


2 


0.1153(5) 


Tree level O(o) improvement. Scale fixed through ro. 
Estimate of lattice artifacts via W^PT in [326] . 


JLQCD/TWQCD 07 


[248] 


2 


0.1111(24) 


Automatic 0(a) impr., exact chiral symmetry. Scale 
fixed through ro. No estimate of lattice artifacts. 


JLQCD/TWQCD 07A 


[245] 


2 


~ 0.12 


Automatic 0(a) impr., exact chiral symmetry. Scale 
fixed through ro. No estimate of lattice artifacts. 


CERN-TOV 06 


[254] 


2 


0.0717(15), 

0.0521(7), 

0.0784(10) 


Scale fixed through Mk- The lattice with a — 
0.0784(10) is obtained with non-perturbative O(a) 
improvement. No estimate of systematic error. 


QCDSF/UKQCD 06A 


[260] 


2 


0.07-0.115 


5 lattice spacings. Non-perturbative 0(a) improve- 
ment. Scale fixed through ro. 



Table 37: Continuum extrapolations/estimation of lattice artifacts in Nf = 2 determinations 
of the Low-Energy Constants. 
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Collab. Ref. Nf a [fm] Description 



ETM 10 [60] 2+1+1 0.078, 0.086 Scale fixed through F v jM v . Systematic error 

taken into account by considering two lattice spac- 
ings. 



MILC 10, 10A [103, 171] 2+1 0.045 - 0.09 3 lattice spacings, continuum extrapolation by 

means of RS^PT. 



JLQCD/TWQCD 10 [240] 2+1, 3 0.11 One lattice spacing, scale fixed through m fl . 



RBC/UKQCD 10A [106] 2+1 0.1106(27), Two lattice spacings. Data combined in global 

0.0888(12) chiral-continuum fits. 



JLQCD 09 [241] 2+1 0.1075(7) Scale fixed through r . Systematic error associ- 

ated to lattice artifacts of order 7.4% estimated 
through mismatch of the lattice spacing obtained 
from different inputs. 

MILC 09, 09A [6, 59] 2+1 0.045 - 0.18 Total of 6 lattice spacings, continuum extrapola- 

tion by means of RSxPT. 



TWQCD 08 [242] 2+1 0.122(3) Scale fixed through m p , r . No estimate of sys- 

tematic effects. 



JLQCD/TWQCD 08B [243] 2+1 0.1075(7) Scale fixed through r . No estimate of systematic 

effects for the observable adopted in this work. 



PACS-CS 08 [63] 2+1 0.0907 Only one lattice spacing, no attempt to estimate 

size of cut-off effects. 



RBC/UKQCD 08 [108] 2+1 0.114 Only one lattice spacing, attempt to estimate size 

of cut-off effects via formal argument. 



RBC/UKQCD 08A [250] 2+1 0.114 Only one lattice spacing, attempt to estimate size 

of cut-off effects via formal argument. 



NPLQCD 06 [176] 2+1 0.125 Only one lattice spacing, continuum xPT used. 



LHP 04 [259] 2+1 ~ 0.12 Only one lattice spacing, mixed discretization ap- 

proach. 



Table 38: Continuum extrapolations/estimation of lattice artifacts in Nf = 2 + 1 and 2 + 1 + 1 
determinations of the Low-Energy Constants. 
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Collab. 


Ref. 


N f 




M., min [MeV] 


Description 


Bernardoni 10 


[244] 


2 




297, 377, 426 


NLO SU(2) fit of the topological sus- 
ceptibility. 


JLQCD/TWQCD 10 


[240] 


2 


mEV < 1 (e-reg.) 

m ~ m s /6 — m s (p-reg.) 


Quark masses both in the p and in the 
e-regime. NLO chiral fit of the spectral 
density interpolating the two regimes. 
( + i 8%) f° r ^ 1 ^ 3 associated to chiral fit. 


JLQCD/TWQCD 09 


[252] 


2 




290 


LECs extracted from NNLO chiral fit 
of vector and scalar radii (r 2 )y s . 


ETM 09B 


[246] 


2 


A* = 


= mEV ~ 0.11 - 0.35 


NLO SU(2) e-regime fit. 


ETM 09C 


[237] 


2 




280 


NNLO SU(2) fit. Systematic error esti- 
mated by performing a large set of dif- 
ferent fits. 


Hj 1 1V1 Uo 


Too c 1 
[ZoDJ 


Z 




oar* 
260 


LECs extracted from pion form factor 
using NNLO x?T and fixing the pion 
scalar radius to the experimental value. 


JLQCD/TWQCD 08 A 
JLK^vjL) UoA 


[95] 

ro£ a l 
[264] 


2 




290 


NNLO SU(2) fit up to At* = 750 MeV. 
Estimate of systematic error due to 
truncation of chiral series. 


CERN 08 


[215] 


2 




mq, m i n =13 MeV 


NLO SU(2) fit for the mode number of 
the Dirac operator. 


HHS 08 


[247] 


2 


/'i 


= mEVi ~ 0.7-2.9 
= mSV 2 =^ 2.1 - 5.0 


e-regime: higher values might be too 
large to be considered in the e-regime. 
NLO SU(2) e-regime fit. 


JLQCD/TWQCD 07 


[248] 


2 




= mSV ~ 0.556 


NLO SU(2) e-regime fit. 


JLQCD/TWQCD 07A 


[245] 


2 




?n u <i = m a /6 — m s 


Quark condensate extracted from topo- 
logical susceptibility, LO chiral fit. 


CERN-TOV 06 


[254] 


2 




403, 381, 377 


NLO SU(2) fit up to M w ~ 540 MeV. 


QCDSF/UKQCD 06A 


[260] 


2 




400 


Several fit functions to extrapolate the 
pion form factor. 



Table 39: Chiral extrapolation/minimum pion mass in Nf = 2 determinations of the Low- 
Energy Constants. 
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Collab. 


Ref. 


Nf 


M^, min [MeV] 


Description 


ETM 10 


[60] 


2+1+1 


270 


SU(2) NLO and NNLO fits. 


MILC 10, 10A 


103, 171] 


2+1 




Cf. MILC 09A. 


JLQCD/TWQCD 10 


[240] 


2+1,3 


m ud T,V < 1 (e-reg) 

In 

m u d = m s /b — m s 

(p-reg.), 

m s ~ phys. 


For the Nf = 2 + 1 runs, light quark 
masses both in the e- and in the p- 
regime. For the Nf = 3 runs, only 
masses in the p-regime. NLO chiral fit of 
the spectral density interpolating the two 
regimes. {_ 02 %) systematic error associ- 
ated to E 1/3 and ±8.7% to eJ /3 . 


RBC/UKQCD 10A 


[106] 


2+1 


290-420 


The range is 225-420 MeV for partially 
quenched pions. NLO SU(2) X PT fits are 
used. 


JLQCD 09 


[241] 


2+1 


m ud Y>V < 1 (e-reg) 
m u d = m s /6 — m s 
(p-reg.), 
m s ~ phys. 


Quark masses both in the e and in the p- 
regime. NLO chiral fit of the spectral den- 
sity interpolating the two regimes. ( + q ^) 
% systematic error estimated by varying 
the fit range and using formulae for SU(2) 
versus SU(3) xPT. 


MILC 09 


[6] 


2+1 


240 


Value taken from text of [6]; according 
to [59] lightest Nambu-Goldstone mass 
is 224 MeV and lightest RMS mass is 
258 MeV (at 0.06 fm). 


MTT.n 09A 




9-1-1 


258 


177 MeV (at 0.09 fm) and lightest RMS 
mass is 258 MeV (at 0.06 fm). 


TWQCD 08 


[242] 


2+1 


m ud = (m s — 
3ra s )/4, m s ~ phys. 


Quark condensate extracted from topolog- 
ical susceptibility, LO chiral fit. 


JLQCD/TWQCD 08B 


[243] 


2+1 


m u d — rn s /6 — m s , 
m s ~ phys. 


Quark condensate extracted from topolog- 
ical susceptibility, LO chiral fit. 


PACS-CS 08 


[63] 


2+1 


156 


To date, lightest published quark mass 
reached in a direct simulation. 


RBC/UKQCD 08 


[108] 


2+1 


330 


Other unitary masses are 419, 557, 672 
MeV, lightest non-unitary pion goes down 
to 242 MeV. 


RBC/UKQCD 08A 


[250] 


2+1 


330 


Pion electromagnetic form factor com- 
puted at one pion mass. Low values 

i* 2 i iii i i I'll 

of q readied thanks to partially twisted 
boundary conditions. LECs determined 
from NLO X PT fit in q 2 . 


NPLQCD 06 


[176] 


2+1 


460 


Value refers to lightest RMS mass at a = 
0.125 fm as quoted in [59]. 


LHP 04 


[259] 


2+1 


318 


Pion form factor extracted from vector 
meson dominance fit. 



Table 40: Chiral extrapolation/minimum pion mass in Nf = 2+1 and 2 + 1 + 1 determinations 
of the Low-Energy Constants. 
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Collab. 


Ref. 


N f 


L [fin] 




Description 


Bernardoni 10 


[244] 


2 


1.88 


2.8 


Finite-volume effects included in the 
NLO chiral fit. 


JLQCD/TWQCD 10 


[240] 


2 


1.8-1.9 




Finite volume effects estimated 
through differences between different 
topological sectors. Systematic error 
of 3.7% on E 1/3 . 


JLQCD/TWQCD 09 


[252] 


2 


1.89 


2.9 


Finite-volume corrections for form 
factor estimated through NLO x?T 
\oz i , ozoj . /vactitionai unite- volume 
effect for fixed topological charge 
[329]. 


ETM 09B 


[246] 


2 


1.3, 1.5 


e-regime 


Topology: not fixed. Comparison be- 
tween 2 volumes, estimate of 8% sys- 
tematic effect (not quoted in the final 
result). Volume might be too small 
for e-cxpansion. 


ETM 09C 


[237] 


2 


2.0-2.5 


4.4, 3.7, 3.2, 
3.5 


Several volumes. Finite- volume ef- 
tects estimated tnrougn [323]. 


ETM 08 


[235] 


2 


2.1, 2.8 


3.4, 3.7 


For pion form factor only data with 
M^L > 4 are considered in the final 
analysis. FSE for pion masses and de- 
cay constants are taken into account 
according to [323]. 


JLQCD/TWQCD 08 A 
JLQCD 08A 


[95] 
[264] 


2 


1.89 


2.9 


Finite-volume corrections estimates 
through [323]. Additional finite- 
volume effect for fixed topological 
charge [329]. 


CERN 08 


[215] 


2 


1.88, 2.51 


- 


For the specific observable stud- 
ied here the finite-volume effects oc 
e -M A V2 ) with Ma ^ Mn _ Finite- 
volume effects checked by comparing 
2 volumes. 


HHS 08 


[247] 


2 


1.84, 2.77 


e-regime 


Topology: not fixed. Comparison be- 
tween 2 volumes. 


JLQCD/TWQCD 07 


[248] 


2 


1.78 


e-regime 


Topology: fixed to v = 0. 


JLQCD/TWQCD 07A 


[245] 


2 


1.92 




Topology fixed to v — gives finite- 
volume effects estimated through 
[329]. 


CERN-TOV 06 


[254] 


2 


1.72, 1.67, 1.88 


3.5, 3.2, 3.6 


No estimate of finite- volume effects. 


QCDSF/UKQCD 06A 


[260] 


2 


1.4-2.0 


3.8 


Use NLO xPT [327] to investigate 
finite- volume effects. 



Table 41: Finite volume effects in Nt = 2 determinations of the Low-Energy Constants. 
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Collab. 


Ref. 


N f 


L [fm] 




Description 


ETM 10 


[60] 


2+1+1 


1.9-2.8 


3.0 


Finite-size effects studied us- 
ing [323]. LECs are extrac- 
ted from data which satisfy 
M„+L > 4, but M^oL ~ 2. 


A/TTT r in in A 


[lUo, 1 ( 1J 


z+1 


Z.oZ 


4.11 


L >_ z.y im tor tne lignter 
masses. 


JLQCD/TWQCD 10 


[240] 


2+1, 3 


1.9, 2.7 




Finite volume effects esti- 
mated by comparing with 
larger volume for a fixed quark 
mass. For F estimated to be 
3%, for S 1 / 3 and £* /3 (+°;^). 


RBC/UKQCD 10A 


[106] 


2+1 


2.7 


~ 4 


Finite-volume corrections in- 
cluded by means of yPT. 


JLQCD 09 


[241] 


2+1 


1.72 


e/p-regime 


Finite-volume effects esti- 
mated by comparing with 
larger volume (L — 2.58 fm) 
for a fixed m u d, m s- 


MILC 09, 09A 


[6 591 


2+1 


2.4/2.9 


3.5/4.11 


Lmin = 2.4 fm, but L > 2.9 fm 
for lighter masses; one coarse 
ensemble with M^L = 3.5, all 
fine ones have M 7r L>4.11. 


TWQCD 08 


[242] 


2+1 


1.95 


- 


No estimate of finite-volume 
effects. 


JLQCD/TWQCD 08B 


[243] 


2+1 


1.72 




Topology fixed to v = gives 
finite- volume effects estimated 
through [329]. 


PACS-CS 08 


[63] 


2+1 


2.9 


2.3 


Possible finite volume effects 
are the main concern of the 
authors. 


RBC/UKQCD 08 


[108] 


2+1 


2.74 


4.6 


Good setting to correct resid- 
ual finite-volume effects by 
means of x?T. 


RBC/UKQCD 08A 


[250] 


2+1 


2.74 


4.6 


Finite-size effects on the form 
factor and radius estimated to 
be < 1% using xPT 


NPLQCD 06 


[176] 


2+1 


2.5 


3.7 


Value refers to lightest valence 
pion mass. 


LHP 04 


[259] 


2+1 


~ 2.4 


3.97 


Value refers to domain-wall 
valence pion mass. 



Table 42: Finite volume effects in Nf = 2 + 1 and 2 + 1 + 1 determinations of the Low-Energy 
Constants. 
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Collab. 


Ref. 


N, 


Description 


ETM 10 


[60] 


2+1+1 


Non-perturbative 


JLQCD/TWQCD 10 


[240] 


2+1, 3 


Non-perturbative 


MILC 10, 10A 


[103, 171] 


2+1 


2 loop 


RBC/UKQCD 10A 


[106] 


2+1 


Non-perturbative 


JLQCD 09 


[241] 


2+1 


Non-perturbative 


MILC 09, 09A 


[6, 59] 


2+1 


2 loop 


TWQCD 08 


[242] 


2+1 


Non-perturbative 


JLQCD/TWQCD 08B 


[243] 


2+1 


Non-perturbative 


PACS-CS 08 


[63] 


2+1 


1 loop 


RBC/UKQCD 08, 08A 


[108, 250] 


2+1 


Non-perturbative 


NPLQCD 06 


[176] 


2+1 




LHP 04 


[259] 


2+1 




All collaborations 




2 


Non-perturbative 



Table 43: Renormalization in determinations of the Low-Energy Constants. 
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B.4 Notes to section 6 on Kaon £>-parameter B K 

In the following, we summarize the characteristics (lattice actions, pion masses, lattice spac- 
ings, etc.) of the recent Nf = 2 + 1 and Nf = 2 runs. We also provide brief descriptions of 
how systematic errors are estimated by the various authors. 



Collab. 


Ref. 


N f 


a [fm] 


Description 


SWME 11 


[281] 


2+1 


0.12, 0.09, 0.06, 0.04 


Continuum extrapolation of results ob- 
tained at four lattice spacings; resid- 
ual discretization error of 0.21% from 
difference to result at smallest lattice 
spacing. 


RBC/UKQCD 10B 


[282] 


2+1 


0.114, 0.087 


Two lattice spacings. Combined chiral 
and continuum fits. 


SWME 10 


[283] 


2+1 


0.12, 0.09, 0.06 


Continuum extrapolation of results ob- 
tained at four lattice spacings; resid- 
ual discretization error of 0.21% from 
difference to result at smallest lattice 
spacing. 


Aubin 09 


[269] 


2+1 


0.09,0.12 


Two lattice spacings; quote 0.3% dis- 
cretization error, estimated from vari- 
ous a 2 -terms in fit function 




[1 08 9841 




0.114(2) 


^i n o'l p 1 pp <st~i£i pino - ' nuritp 4% H 1 <3— 

cretization error, estimated from the 
difference between computed and ex- 
perimental values of f n . 


HPQCD/UKQCD 06 


[285] 


2+1 


0.12 


Single lattice spacing; 3% discretiza- 
tion error quoted without providing 
more details. 


ETM 10A 


[286] 


2 


0.07, 0.09, 0.1 


Three lattice spacings; 1.2% error 
quoted. 


JLQCD 08 


[280] 


2 


0.118(1) 


Single lattice spacing; no error quoted. 


RBC 04 


[287] 


2 


0.117(4) 


Single lattice spacing; no error quoted. 


UKQCD 04 


[288] 


2 


0.10 


Single lattice spacing; no error quoted. 



Table 44: Continuum extrapolations/estimation of lattice artifacts in determinations of Br. 
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Collab. 


Ref. 




M,, mij) [MeV] 


Description 


SWME 11 


[281] 


2+1 


309, 292, 275, 
334 


The M^min entries correspond to the four lattice 
spacings. Chiral extrapolations based on SU(2) 
staggered xPT at NLO, including some analytic 
NNLO terms. Combined 1.1% error from various 
different variations in the fit procedure. 


RBC/UKQCD 10B 


[282] 


2+1 


327, 289 


The M^min entries correspond to two lattice spac- 
ings. Combined chiral and continuum extrapola- 
tions. 


SWME 10 


[283] 


2+1 


309, 292, 275 


The M^min entries correspond to the three lattice 
spacings. Chiral extrapolations based on SU(2) 
staggered xPT at NLO, including some analytic 
NNLO terms. SU(3) staggered xPT as cross- 
check. Combined 1.1% error from various differ- 
ent variations in the fit procedure. 


Aubin 09 


[269] 


2+1 


240/370 


The M^min entries correspond to the smallest va- 
lence/sea quark masses. Chiral & continuum fits 
based on NLO mixed action x?T at NLO, includ- 
ing a subset of NNLO terms. Systematic error 
estimated from spread arising from variations in 
the fit function. Relevant pion mass in the sea 
quark sector is the singlet pion mass of 370 MeV. 


RBC/UKQCD 07A, 08[108, 284] 


2+1 


330 


Fits based on SU(2) PQxPT at NLO. Effect of 
neglecting higher orders estimated at 6% via dif- 
ference between fits based on LO and NLO ex- 










pressions. 


HPQCD/UKQCD 06 


[285] 


2+1 


360 


3% uncertainty from chiral extrapolation quoted, 
without giving further details. 


ETM 10A 


[286] 


2 


300, 270, 400 


Each M Wtin i n entry corresponds to a different lat- 



tice spacing. Simultaneous chiral & continuum 
extrapolations, based on x?T at NLO, are car- 
ried out. Systematic error from several sources, 
including lattice calibration, quark mass calibra- 
tion, chiral and continuum extrapolation etc., es- 
timated at 3.1%. 



JLQCD 08 [280] 2 290 Fits based on NLO PQxPT. Range of validity in- 

vestigated. Fit error included in statistical uncer- 
tainty. 

RBC 04 [287] 2 490 Fits based on NLO PQxPT. Fit error included in 

statistical uncertainty. 

UKQCD 04 [288] 2 780 Fits to continuum chiral behaviour at fixed sea 

quark mass. Separate extrapolation in sea quark 
mass. Fit error included in overall uncertainty. 



Table 45: Chiral extrapolation/minimum pion mass in determinations of Bk- 
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Collab. 



Ref. Nf L [fm] 



Af-^minL Description 



Each L-entry corresponds to a dif- 
ferent lattice spacing, with two 
volumes at the coarsest lattice. 
Central value obtained via finite- 
volume staggered xPT- Residual 
uncertainty of 0.6% from difference 
obtained using xPT in infinite vol- 
ume. 

RBC/UKQCD 10B [282] 2+1 2.7 > 4.0 Finite-volume xPT applied; resid- 

ual finite-volume error estimated 
from the difference to fit to xPT in 
infinite volume. 



Each L-entry corresponds to a dif- 
ferent lattice spacing, with two vol- 
umes at the coarsest lattice. Finite- 
volume error of 0.9% estimated 
from difference obtained on two vol- 
umes at the coarsest lattice spacing. 



Aubin 09 [269] 2+1 2.4, 3.4 3.5 Each L entry corresponds to a 

different lattice spacing. Keep 
m w L > 3.5; no comparison of re- 
sults from different volumes; 0.6% 
error estimated from mixed action 
xPT correction. 



RBC/UKQCD 07A, 08[108, 284] 2+1 1.83/2.74 4.60 Each L entry corresponds to a dif- 

ferent volume at the same lattice 
spacing; 1% error from difference in 
results on two volumes. 



HPQCD/UKQCD 06 [285] 2+1 2.46 4.49 Single volume; no error quoted. 



ETM 10A [286] 2 2.2, 2.2/2.9, 3.3, 3.3/4.3, Each L entry corresponds to a dif- 

2.1 5 ferent lattice spacing, with two vol- 

umes at the intermediate lattice 
spacing. Results from these two 
volumes at ~ 300 MeV are 

compatible. 



JLQCD 08 [280] 2 1.89 2.75 Single volume; data points with 

?n V ai < fisca excluded; 5% error 
quoted as upper bound of PQxPT 
estimate of the effect. 



RBC 04 


[287] 


2 


1.87 


4.64 


Single volume; no error quoted. 


UKQCD 04 


[288] 


2 


1.6 


6.51 


Single volume; no error quoted. 



Table 46: Finite volume effects in determinations of Bk- 



SWME 11 [281] 2+1 2.4/3.3, 2.4, > 3.5 

2.8, 2.8 



SWME 10 [283] 2+1 2.4/3.3, 2.4, > 3.5 

2.8 
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Collab. 


Ref. 


N f 


Ren. 


running 
match. 


Description 


sWMt 11 


[281] 


2+1 


PTU 


PTU 


Uncertainty from neglecting higher orders 
estimated at 4.4% by identifying the un- 
known 2-loop coefficient with result at the 
smallest lattice spacing. 


RBC/UKQCD 10B 


[282] 


2+1 


RI 


PTU 


Variety of different RI-MOM schemes in- 
cluding non-exceptional momenta. Resid- 
ual uncertainty of 2% uncertainty in run- 
iiing 06 iiiaLLiiiiig. 


SWME 10 


[283] 


2+1 


PTU 


PTU 


Uncertainty from neglecting higher orders 
estimated at 5.5% by identifying the un- 
known 2-loop coefficient with result at the 
smallest lattice spacing. 


Aubin 09 


[269] 


2+1 


RI 


PTU 


Total uncertainty in matching & running of 
3.3%, estimated from a number of sources, 
including chiral extrapolation fit ansatz 
for n.p. determination, strange sea quark 
mass dependence, residual chiral symme- 
try breaking, perturbative matching & 
running. 


RBC/UKQCD 07A, 08 [108, 284] 


2+1 


RI 


PTU 


Uncertainty from n.p. determination of 
ren. factor included in statistical error; 2% 
systematic error from perturbative match- 
ing to MS estimated via size of correction 
itself. 


HPQCD/UKQCD 06 


[285] 


2+1 


PTU 


PTU 


Uncertainty due to neglecting 2-loop order 
in perturbative matching and running es- 
timated by multiplying result by a . 


ETM 10A 


[286] 


2 


RI 


PTU 


Uncertainty from RI renormalization esti- 
mated at 2.5%. 


JLQCD 08 


[280] 


2 


RI 


PTU 


Uncertainty from n.p. determination of 
ren. factor included in statistical error; 
2.3% systematic error from perturbative 
matching to MS estimated via size of cor- 
rection itself. 


RBC 04 


[287] 


2 


RI 


PTU 


Uncertainty from n.p. determination of 
ren. factor included. 


UKQCD 04 


[288] 


2 


ptu 


PTU 


No error quoted. 



Table 47: Running and matching in determinations of Br. 
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